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i.o  iitroductior 


The  "MINI"  Electromagnetics  Cod* ,  or  MlRIREC,  is  a  method  of  moewnts  r.os- 
put*r  prograa  for  analysis  of  thin  wlr*  antennas  (reference  1  >.  A  Galerkm 
procedure  Is  applied  to  an  electric  field  integral  equation  to  solve  for  the 
wlr*  currents  following  an  approach  suggested  by  Wilton  (reference  2).  This 
formulation  results  in  an  unusually  short  computer  prograi  suitable  for 
implementation  on  a  microcomputer.  Hence,  MIRIRIC  is  written  in  •  BAS  I’' 
language  compatible  with  many  popular  microcomputers. 

MlRIREC  solves  for  impedance  and  currents  on  arbitrarily  oriented  wires, 
including  configuration*  with  multiple  wlr*  Junctions,  in  free  space  and  over 
a  perfectly  conducting  ground  plan*.  Options  include  lumped  parameter 
impedance  loading  of  wires  and  calculation  of  near  zone  and  far  zone  fields. 
Both  near  electric  fields  and  near  magnetic  fields  can  be  determined  for  free 
space  and  over  a  perfectly  conducting  ground.  The  far  zone  electric  fields 
and  radiation  pattern  (power  pattern)  can  also  be  determined  for  free  space 
and  perfectly  conducting  ground. 

Additional  radiation  pattern  options  include  a  Fresnel  reflection 
coefficient  correction  to  the  patterns,  for  finite  conducting  grounds  (real 
earth  surface  Impedance).  Up  to  five  changes  in  surface  impedance  due  to  real 
ground  are  allowed  in  a  linear  or  circular  "cl  iff"  model.  The  cliff  may  take 
on  any  elevation  (including  zero,  i.*.,  a  flat  surface),  however,  there  is  no 
correction  for  diffraction  from  cliff  edges.  In  the  case  of  a  circular  cliff 
model,  the  first  media  may  include  a  correction  for  the  surface  impedance  of  a 
densely  spaced,  burled,  radial  wlr*  ground  screen. 

The  first  version  of  WIRIREC  given  by  ROSC  TD  516  (reference  i),  calcu¬ 
lated  currents  and  radiation  patterns  for  wire  antennas  in  free  space  and  over 
a  perfectly  conducting  ground  plane.  Wires  attached  to  ground  were  required 
to  intersect  at  a  right  angle  and  could  not  be  impedance  loaded  at  the  connec¬ 
tion  point.  Subsequent  revisions  corrected  these  shortcomings  culminating  ir, 
version  2  or  MIRIREC(2) ,  given  by  Li,  et  al  .  (reference  3).  All  previous 
versions  of  MIRIREC  require  user  specification  of  wire  end  connections. 
However,  MIHIREC(3)  determines  connection  information  for  itself  from  user 
defined  wire  end  coordinate?.  MIWIWEC(3)  also  displays  the  currents  wire  t> 
wire,  and  at  all  wire  ends,  including  wire  Junctions.  MIRIREC(3>  features  ar. 
Improved,  faster  solution  routine  and  has  been  completely  restructured  using  a 
more  modular  programming  style,  including  the  use  of  helpful  comment 
statements. 

1 . 1  BACKGROURD 

The  Rumerical  Electromagnetics  Code  (NEC)  found  in  reference  ^  is  the 
most  advanced  computer  code  available  for  the  analysis  of  thin  wire  antennas. 
It  is  a  highly  user-oriented  computer  code  offering  a  comprehensive  capability 
for  analysis  of  the  interaction  of  electromagnetic  waves  with  conducting  struc¬ 
tures.  The  program  is  based  on  the  nunerical  solution  of  integral  equations 
for  the  currents  induced  on  the  structure  by  an  exciting  field. 

REC  combines  an  integral  equation  for  smooth  surfaces  with  one  for  wires 
to  provide  convenient  and  accurate  modeling  for  a  wide  range  of  applications. 
A  REC  model  may  include  nonradiating  networks  and  transmission  lines,  perfect 
and  imperfect  conductors,  limped  element  loading,  and  ground  planes.  The 
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ground  planes  say  b«  perfectly  or  imperfectly  conducting.  Eicitatior,  may  be 
vie  en  applied  voltage  ecurce  or  incident  plane  wave  .  The  output  aay  int.ude 
induced  currents  and  charges,  near  ur  far  zone  eiectru  ur  magnetic  fields, 
and  uspedance  or  admittance.  Many  other  cceaaonly  used  parameters  sue  r  as  gair. 
and  directivity,  power  be  Iget ,  and  antenna  to  antenna  coupling  are  iiic 

av  a  1 1  at  1  e  . 

MET  is  a  powerful  tool  for  eany  engineering  appi  l  c  at  i  ora .  It  is  idea, 

for  Modeling  co-site  antenna  environments  ir  which  the  interact  ior.  between 
antenna  and  environment  cannot  be  Ignored  .  In  aary  protlems.  however,  the 
estenslve  full  capability  of  NEC  is  not  really  required  becaure  the  antenna 
and  its  environment  are  not  very  complex  or  the  lnforstatior  sought  requires 

only  a  simplified  mcoel .  In  addition,  NEC  requires  the  support  of  and  access 
to  a  .  arge  main-frame  computer  system.  These  computer  systems  are  esper.sive 
and  not  always  readily  available  at  remote  field  activities.  Ever  wher.  the 

computer  facilities  are  available,  heavy  demand  usage  say  result  ir  s.ow  turn¬ 
around .  ever,  for  relatively  simple  or  small/  NEC  runs.  One  viat.e  solution 
is  a  "stripped  down"  version  of  NEC  that  would  retain  only  the  basic  solution 
and  the  most  frequently  used  options  and  which  could  be  implemented  or  a  mir.i- 
cr  microcomputer  with  an  advanced  FORTRAN  language  capability.  M1NINEC'  • 
offers  many  of  the  required  NEC  options,  but  makes  use  of  a  BASIC  language 
that  is  compatible  with  many  popular  microcomputers.  MININEC  j,  is  only  suit¬ 
able  for  small  problems  less  than  7C,  unknowns  and  wires,  depending  or  the 
computer  memory  end  BASIC  compiler/. 

COMPUTER  RECvb  IRENEMTS 

Occasionally  a  technology  develops  which  is  destined  to  produce  signifi¬ 
cant  changes  In  the  way  people  think  and  conduct  their  business.  For  many 
decades,  scientists  and  engineers  struggled  with  unmanageable  equations  and 
lata  using  trial  and  error  techniques,  employing  logarithmic  tables  and  in¬ 
adequate  slide  rule  calculations.  Then  came  the  digital  computer. 

In  the  ’ 9Fs0s  and  fcCs,  physically  .arge  and  eipensive  computing  machines 
that  were  relatively  slow,  with  limited  capability  compared  tc  today’s  stan¬ 
dards-  became  aval. able  to  a  few.  At  first,  stored  programs  were  accessible 
through  direct  connection  of  individual  terminals  a  short  distance  away.  The 
revolution  had  begun. 

In  the  ''Os.  technologists  rushed  to  convert  prover  algorithms  ir.tc  oinpu- 
ter  programs  or  tc  develop  new  algorithms  suitable  for  eff::iert  computer 

programming  for  use  as  analysis  and  synthesis  tools  by  the  scientific  com¬ 

munity.  These  tools,  for  the  most  part,  required  the  support  of  large  central 
machines.  Meanwhile,  slide  rules  were  being  replaced  by  hand-held  calculators 
with  tr  igonometr ic  functions,  some  of  which  could  be  programmed  for  simple 
repetitive  algorithms. 

Today,  large  central  processing  systems  are  being  supplemented  with  sma.  . 
powerful  mini-  and  microcomputers.  The  development  of  the  low  cost  micro¬ 
processor  chip  means  that  computers  with  capabilities  that  equal  or  exceed 
those  of  the  earlier  mair,  frame  machines  of  the  cr-i  are  now  available  ;r 

compact  size.  f-lzes  range  from  suitcase,  or  desktop,  machines  the  rr, .  r,  - 

compu'er  hr  file  ^at  met  machines  the  minicomputer  that  ar  te  expanded  r 
configured  to  meet  specialized  needs.  The  microcomputer  is  becoming  more  and 


■on  a  f  for  debit  •■  a  personal  coapu  inf  tool.  Tha  microcomputer ,  or  "hoae 
coaputer* ,  is  emerging  aa  today's  aoat  important  engineering  and  scientific 
tool,  al lotting  widespread  networking.  Anyone  with  a  alcrocoaputer  or  teralnal 
with  an  acoustic  coupler  and  telephone  has  access  to  a  wide  variety  of 
computing  facilities  around  the  country,  as  well  aa  an  alaoat  llaltless  source 
of  information. 

NININEC  has  been  written  with  the  alcrocoaputer  in  aind.  But,  It  can 
also  be  laplaaented  on  alnl-  or  larger  coaputera  that  have  the  BASIC  language 
capability.  However,  soae  changes  In  the  prograa  aay  be  required.  Program¬ 
ming  haa  been  kept  slaple.  with  few  machine -dependent  prograa  atateaenta,  ao 
that  it  will  be  compatible  with  aoat  BASIC  languagea. 

NEC  la  aultable  for  both  aaal  1  and  large  numerical  model  a.  The  upper 
Halt  la  determined  by  the  coat  factora  and  memory  size  of  the  aalnfraae  on 
irfiich  It  resides.  A  aodel  containing  up  to  2000  unknowns  (segamnts)  aeeaa  to 
be  the  practical  upper  Halt.  On  the  other  hand,  NININEC  la  aultable  only  for 
saMll  problems.  The  upper  limit  la  determined  by  the  memory  size  and  speed  of 
the  alcrocoaputer  eaployed.  Practical  Halts  aeaa  to  be  30  to  AO  unknowns 
(current  pulses)  when  using  interpreter  BASIC,  due  to  the  tlae  required  to 
obtain  a  solution.  However,  if  one  la  willing  to  wait  an  hour  or  aw>re  for  the 
solution,  s  aodel  with  65  to  lb  unknowns  is  possible.  Serious  antenna 
modeling  requires  the  use  of  a  BASIC  compiler.  In  addition,  a  aath 
co-processor  board  is  recommended.  Present  alcrocoaputer  memory  size  Halts 
MINIMEC  to  aodels  with  less  than  100  unknowns.  For  problems  of  100  or  more 
unknowns,  a  main  Cream  is  recaaaended ,  and  in  that  case,  the  use  of  NEC  Is  the 
natural  choice. 
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2.0  THE  THEORY  OF  MININEC 


The  MININEC  program  is  based  on  the  numerical  solution  of  an  integral 
equation  representation  of  the  electric  fields.  Discussion  of  similar  formu¬ 
lations  can  be  found  elsewhere,  for  example,  see  Harrington  (reference  5). 
The  real  advantage  is  that  the  solution  technique  as  implemented  in  MININEC 
results  in  a  relatively  compact  (i.e.,  short)  computer  code.  The  discussion 
that  follows  in  this  section  is  condensed  from  reference  2. 

2.1  THE  ELECTRIC  FIELD  INTEGRAL  EQUATION  AND  ITS  SOLUTION 

It  has  become  customary  in  solving  wire  antenna  problems  to  make  several 
assumptions  which  are  valid  for  thin  wires.  They  are  that  the  wire  radius,  a, 
is  very  small  with  respect  to  the  wavelength  and  the  wire  length.  Because  it 
is  necessary  to  subdivide  wires  into  short  segments,  the  radius  is  assumed 
small  with  respect  to  the  segment  lengths  as  well,  so  that  the  currents  can  be 
assumed  to  be  axially  directed;  i.e.,  there  are  no  azimuthal  components  of 
current . 

Figure  1  gives  the  geometry  of  a  typical,  arbitrarily  oriented  wire. 
Assume  that  the  wire  is  straight,  even  though  the  theory  applies  equally  to 
bent  configurations.  The  same  wire  is  also  shown  broken  into  segments  or 
subsect i ons . 


In  equations  (1),  (2),  and  (3)  below,  the  vector  and  scalar  potentials 
are  given  by 
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The  integral  equation  relating  the  incident  field,  Ejnc«  and  the  vector 
and  scalar  potentials  is 

“  X  “  A  A 

-Einc  •  S  a  - JuA  •  S  -  S  •  V*  .  ( 4 ) 

Equation  (4),  above,  is  solved  in  MININEC  by  using  the  following  procedure. 

The  wires  are  divided  into  equal  segments,  and,  as  shown  in  Figure  1,  the 
vectors  r  ,  n=0,  1,.  .  .N+1  are  defined,  with  respect  to  the  global  coordinate 
origin,  (/)  The  unit  vectors  parallel  to  the  wire  axis  for  each  segment  shown 
are  defined  as 
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Pulse  testing  and  pulse  expansion  functions  used  in  MININEC  are  defined  as 
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where  the  points  sn+1/2  designat®  segment  midpoints, 


s_  ,  s_ 

rv+ 1  n 


nel/?  2 

;r  terms  of  the  global  coordinates, 


(7) 


r  ,  ♦  r 
n+ 1  n 


(8) 


It  is  assumed  that  the  components  of  the  vectors  Einc  and  A  in  equation 

-  are  sufficiently  smooth  over  each  segment  that  their  respective  values  on 

-ji  *  segment  may  be  replaced  by  those  taken  at  the  point  s  .  The  pulse  func- 

*  rs  jf  6  are  then  used  as  testing  functions  on  (4),  resulting  in 
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The  currents  are  expanded  in  pulses  centered  at  the  Junctions  of  adjacent 
segments  as  illustrated  in  Figure  2(a).  Note  that  pulses  are  omitted  from  the 
wire  ends.  This  is  equivalent  to  placing  a  half  pulse  of  zero  amplitude  at 
each  end,  thus  imposing  the  boundary  condition  for  zero  current  at  unattached 
wire  ends.  The  current  expansion  can  be  written  as 

N 

I(s)  =  y  I  P  (s)  .  (11) 
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A  difference  approximation  is  applied  to  equation  (3)  to  compute  the 
charge.  Thus,  as  shown  in  Figure  2(b),  the  charge  can  be  represented  as 
pulses  displaced  from  the  current  pulses  by  a  half  pulse  width. 
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Figure  2  Wire  segmentation  scheme  cllujt'at  nq  equally  weighted  Dulses  'o>  .  u'lf’t  and  i  hjiqe 

Substituting  (11)  into  (10)  produces  a  system  of  equations  that  can  be 
expressed  in  matrix  form.  Each  matrix  element,  Z  ,  associated  with  the  n-th 
current  and  the  s  observation  point  involves  fiP.  lar  and  vector  potential 
terms  with  integral^  of  the  form 
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ru  *  ((»B  -  s’)2  ♦  4a2  sin2  |  ) 


Equation  (12)  do«s  not  land  itself  to  straightforward  integration  becauae  of 
the  singularity  at  r>0.  The  1/r  can  be  subtracted  froai  the  integrand  and  then 
added  as  a  separate  tens  to  yield 


-i  j  *  .  ji  j 


i  *  -Jkr 
1  e  ■  -  1 


The  first  term  of  (15)  can  be  rewritten  as  an  elliptic  integral  of  the  first 
kind  (reference  6). 


Ti'il-  •)  'S  /  ? 


where 


r  2  2i 1/2 

[(s^  -  s'r  ♦  4a* ' 

F(-^  , s )  has  an  approilaation  (reference  6) 
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This  equation  is  substituted  into  equation  (12)  and  written  as 
s 

V 

/  k( s-s' )  ds'  x  ♦  I2  ♦  I3  . 


(20) 


I1t  1^.  and  I3  are  defined  as 


n  / 


(21) 
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Similarly, 


This  integral  has  a  well  behaved  integrand  and  can  be  integrated  nuuerically . 
The  integration  is  broken  up  into  two  integrals  over  the  ranges  (a  s)  and 
(s,s  )  for  best  accuracy.  Gaussian  quadrature  is  used  for  the  numerical  inte¬ 
gration  (reference  7).  The  number  of  points  used  in  the  integration  routine 
is  automatically  selected  by  consideration  of  the  pulse  accuracy  required  for 
the  source  to  observation  distance.  The  final  integral  is 


(23) 


The  integrand  is  nonsingular  and  can  be  integrated  numerically.  To  obviate 
the  need  for  double  integration,  it  is  convenient  to  approximate  the  integral 
by  replacing  r  by  a  reduced  kernel  approximation  of  equation  ( 1  *4 ) . 

Thus 


The  integral  can  be  integrated  numerically  by  the  same  procedure  as  for  I^. 

Thus,  equation  (12)  with  its  singularity  problem  is  evaluated  by  adding 
1^  of  equation  (21),  of  equation  (22),  and  1^  of  equation  (24). 

This  approach  to  evaluate  (12)  is  accurate  for  a  wide  range  of  wire  radii 
but  breaks  down  when  the  radius  becomes  very  small.  For  very  small  radii, 
equation  (12)  may  be  expressed  as  a  single  integral  and  evaluated  using  two 
terms  of  a  Maclaurin  series,  after  Harrington  (reference  5).  This 
approximation  for  the  ♦  terms  is: 


(25) 
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Figure  3  demonstrates  the  range  of  validity  with  and  without  the  small  radius 
correction.  WlthoMt  the  correction,  MININEC  gives  acceptable  answers  for  wire 
radii  between  10~  and  10~  wave  lengths.  Note  that  MININEC  is  within  10%  or 
better  of  th^data  published  by  King  (references  8  and  9),  for  radii  between 
10"  and  10"  wave  len|ths.  The  small  radius  correction  provides  correct 
results  for  radii  of  10”  wave  lengths  or  smaller.  In  MININEC,  the  switch  to 
the  small  radius  approximations  occurs  automatically  for  radii  of  10~  and 
smaller. 


By  substitution,  the  matrix  equation  to  be  solved  is 


[Z  ]  [I  ]  =  IV] 
mn  n  m 


(26) 


where 
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and 
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=  E,  (s  ) 
inc  m 


( r  —  r  ) 

v  mi- 1/2  m-1/2' 


(28) 


U  ]  ia  a  square  matrix  and  [In]  and  t are  column  matrices  with  n=1,2,...N 
and  m*1,2,...N  for  N  total  unknowns  (N  is  the  total  number  of  current  pulses). 
The  extension  of  these  equations  to  two  or  more  coupled  wires  follows  the  same 
line  of  development  and  will  not  be  covered  here. 


The  column  vector  [V  ]  represents  an  applied  voltage  that  superimposes  a 
constant  tangential  elecfric  field  along  the  wire  for  a  distance  of  one 
segment  length  centered  coincident  with  the  location  of  the  current  pulses. 
Hence,  for  a  transmitting  antenna,  all  elements  of  [V  ]  are  set  to  zero  except 
for  the  element(s)  corresponding  to  the  segment(s)  located  at  the  desired  feed 
point(s).  For  an  incident  plane  wave,  all  elements  of  [V  ]  must  be  assigned  a 
value  depending  on  the  strength,  polarization  (or  orientation) ,  and  angle  of 
incidence  of  the  plane  wave.  The  applied  voltage  source  (transmit  case), 
however,  is  the  only  ready-made,  or  programmed,  option  in  MININEC. 

As  stated  above,  the  [Z .  ]  matrix  in  equation  (26)  is  filled  by  the 
evaluation  of  an  elliptic  Integral  and  use  of  Gaussian  quadrature  for 
numerical  integration.  The  solution  of  (26)  can  be  accomplished  by  using  any 
one  of  a  number  of  standard  matrix  solution  techniques.  MININEC(3)  uses  a 
triangular  decomposition  (LU  decomposition)  with  the  Gauss  elimination 
procedure  with  partial  pivoting  (reference  7). 

2.2  WIRE  JUNCTIONS 


The  theory  developed  thus  far  for  straight  wires  is  equally  applicable  to 
bent  wires.  However,  for  coding  simplicity  in  MININEC,  bent  wires  are  treated 
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in  the  same  way  as  the  junctions  of  multiple  numbers  of  wires.  That  is,  a 
bend  in  an  otherwise  straight  wire  is  treated  as  the  junction  between  two 
straight  wires. 

It  has  been  generally  accepted  that  the  currents  at  junctions  of  thin 
wires  conform  to  Kirchoff's  current  law  (reference  10).  Rather  than  expli¬ 
citly  enforcing  this  condition  in  MININEC,  an  overlapping  segment  scheme 
(reference  11)  is  employed  at  junctions  of  two  or  more  wires.  A  detailed 
discussion  of  this  approach,  including  arguments  for  validity,  appears  in  both 
references  8  and  9.  Only  those  aspects  essential  to  the  use  of  MININEC  are 
discussed  here. 

Consider  a  wire  having  no  connections  at  either  end.  The  wire  is  subd¬ 
ivided  into  segments  and  the  current  is  expanded  in  pulses  centered  at  adja¬ 
cent  segment  junctions  as  described  above  and  illustrated  in  Figure  ‘•(a).  The 
end  points  have  no  pulses,  or  alternatively  the  end  points  have  half  pulses 
with  zero  amplitude.  A  second  wire  is  to  be  attached  to  one  end  of  the  first. 
The  second  wire  is  subdivided  into  segments  with  pulses  for  currents  located 
as  in  the  first  case.  However,  a  full  pulse  is  located  at  the  attachment  end, 
with  half  the  pulse  extending  onto  wire  two,  and  half  onto  wire  one,  as  illus¬ 
trated  in  Figure  4(b).  The  half  on  wire  one  assumes  the  dimensions  (length 
and  radius)  of  the  half  segment  on  wire  one,  while  the  half  on  wire  two 
assumes  the  dimensions  appropriate  to  wire  two.  Wire  two  overlaps  onto  wire 
one  with  a  full  pulse  centered  at  the  junction  end.  Note  that  the  free  end  of 
the  wire  has  a  zero  half  pulse.  A  third  wire  may  be  assumed  to  also  overlap 
onto  wire  one,  as  illustrated  in  Figure  4(c).  It  can  be  shown  (see  references 
8  and  9)  that  for  a  junction  of  N  wires,  only  N-1  overlapping  pulses  are 
required  to  satisfy  Kirchoff’s  current  law.  Alternatively,  wire  three  could 
have  overlapped  onto  wire  two  (not  illustrated  here). 

The  convention  in  MININEC(3)  is  that  the  overlap  occurs  onto  the  earliest 
wire  specified  at  a  given  junction.  It  is  assumed  that  a  wire  can  overlap 
onto  another  wire,  provided  that  another  wire  was  previously  specified.  It 
cannot  overlap  onto  a  wire  not  yet  specified.  Either  end  of  a  wire  may 
overlap  onto  either  end  of  another  wire.  All  that  is  required  to  impose  the 
continuity  conditions  at  the  junction  is  that  there  be  N-1  overlaps  for  a 
Junction  of  N  wires. 

Current  reference  directions  are  assumed  to  be  based  on  the  order  in 
which  the  coordinates  of  a  wire  are  specified.  A  positive  wire  current  is 
from  the  end  first  specified,  end  one,  towards  the  other  end,  end  two.  By  use 
of  Kirchoff's  current  law  and  the  current  reference  direction,  the  currents  at 
the  Junction  can  be  found.  For  example,  suppose  the  wires  in  Figure  4(c)  are 
all  specified  from  left  to  right.  Let  the  pulse  amplitudes  for  the  first 
pulse  on  wires  two  and  three  be  I_  and  I  ,  respectively.  Then  the  currents 
out  of  the  Junction  into  wires  two  and  three  are  the  complex  amplitudes  of  the 
first  pulses,  the  overlapping  pulses,  on  wires  two  and  three,  respectively. 
Hence,  the  current  on  wire  one  into  the  junction  is  the  sum  of  these  currents; 
i  .e. ,  1 1  =  l2  *  1^. 


(a)  Wire  one  with  no  end  connection*. 


(b)  Wire  two  overlaps  onto  wire  one. 


Figure  4.  Illustration  of  the  overlap  scheme  used  at  multiple  wire  junctions. 


MININEC(J)  automatically  determines,  during  geometry  input,  whether  there 
is  a  connection  on  either  end  of  a  wire,  and  if  so,  to  which  wire  and  wire  end 
it  is  connected.  After  solving  for  the  current  pulse  amplitudes,  MININEC(3) 
then  computes  the  Junction  currents,  if  any,  for  each  wire  end.  The  final 
display  indicates  free  ends  by  the  letter  E  (for  free  end)  and  junction  ends 
by  the  letter  J.  The  geometry  and  currents  are  displayed  wire  by  wire. 

2.3  THE  GROUND  PLANE 


The  method  of  images  is  used  in  MININEC  to  solve  for  currents  in  wires 
located  over  a  perfectly  conducting  ground  plane. 


Consider  a  wire  structure  represented  by  N  segments.  In  the  presence  of 
a  perfectly  conducting  ground  plane,  by  image  theory,  the  structure  and  ground 
plane  may  be  replaced  by  the  original  structure  and  its  image.  Hence,  there 
are  now  2N  segments  and  2N  unknowns  to  be  determined.  Equation  (26)  can  be 
written  as 


*  ?  1 , 2N 


•?N,2N 


2N ,  2N 


‘  I.  * 
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(29) 


The  image  current,  I  .  ...I  ,  are  equal 

structure,  I  ...I  ,  so*  that  I  =  I 

IN  n  2N-n*1 


to  the  currents  on  the  original 
Half  the  equations  represented  in 
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(29)  contain  redundant  information  and  may  be  dlacarded.  It  may  be  reduced  to 
a  square  matrix  again  by  adding  appropriate  colunns;  i.e.,  by  using  the 
current  identity.  Hence,  (29)  becomes 


V 


*  * 


(30) 


where  Z|j 


Zij  *  Zi  2N-J+1  * 


For  a  wire  attached  to  ground,  a  current  pulse  is  automatically  added  to 
the  wire  end  point  connected  to  ground  so  that  current  continuity  with  its 
image  is  observed;  i.e.,  a  non-zero  half  pulse  is  placed  on  both  the  wire  end 
and  its  image.  The  voltage  in  equation  (30)  is  divided  by  two  in  this  case. 
Either  end  of  a  wire  may  be  attached  to  ground. 


2.4  LUMPED  PARAMETER  LOADING 


The  wire  structures  discussed  so  far  consist  of  perfectly  conducting 
wires.  If  an  impedance  due  to  a  fixed  load,  Z  =  R  ♦  JX,  is  added  to  the 

structure  so  that  its  location  coincides  with  chat  of  one  or  more  of  the 

non-zero-current  pulse  functions  (i.e.,  a  lumped  load  is  placed  on  the  wire  at 

the  junction  of  two  segments),  then  the  load  introduces  an  additional  voltage 

(a  voltage  drop)  equal  to  the  product  of  the  current  pulse  magnitude  and  Z^. 
Hence,  equation  (26)  becomes 

[Z*  ][I  ]  =  tV.l  (3D 

mn  n  m 


where  Z'  =  Z  for  m  i  n  and  Z’  =  2  ♦  Z,  for  m  r  n.  Hence,  a  specified 

impedance  represented  as  the  sun  of  a  resistance  and  a  reactance,  and  located 
on  a  wire  coincident  with  a  current  pulse  is  simply  added  to  the  diagonal 
impedance  element  or  self-term  corresponding  to  that  pulse.  A  distributed 
impedance  such  as  wire  conductivity  can  be  treated  in  the  same  way  by  use  of 
an  equivalent,  lunped-circuit ,  element- impedance  relationship. 

2.5  NEAR  FIELDS 


The  electric  near  fields  and  the  magnetic  near  fields  can  be  determined 
from  the  current  distribution  obtained  in  the  solution  of  equation  (26). 

The  near  electric  fields  are  computed  by  the  method  described  by  A.  T. 
Adams,  et  al . ,  (reference  12).  Using  MININEC ,  the  current  on  the  wire  struc¬ 
ture  is  approximated  using  the  computed  current  pulses.  To  determine  the 
electric  field  at  a  given  point  in  the  near  field,  a  small,  virtual  thin-wire 
dipole  is  placed  at  the  point  with  its  axis  parallel  to  the  appropriate  vector 
component.  The  open-circuit  voltage  at  the  near  field  point  can  be  calculated 
from  a  knowledge  of  the  current  distribution  over  the  wire  structure  and  the 
mutual  impedance  between  the  wire  structure  and  the  virtual  dipole.  In  other 
words , 
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The  virtual  dipole  la  open-circuited.  V  la  the  open-circuit  voltage.  1  are 
the  MININEC  computed  current  pulaea  of  tne  wire  atructure.  I  are  the  mutual 
impedance a  between  the  wire  atructure  and  the  virtual  dipole.  The  mutual 
impedances  are  calculated  using  the  MIN1NEC  method  of  equation  (21).  The 
electric  field  strength  along  the  direction  of  the  virtual  dipole  la  given  by 


£  .  _  - “ - 

d  length  of  dipole  . 


This  equation  is  evaluated  once  for  each  electric  field  vector  component  in 
the  x,  y  and  z  directions  at  the  near  field  point  of  interest.  In  MININEC,  a 
virtual  dipole  of  length  .001  wave  length  is  used. 

MININEC  calculates  the  three  vector  components,  E  ,  E  and  E  ,  as  real 
and  imaginary  terms,  from  which  the  magnitude  and  phase*  are^  determfned.  The 
average  value  is  determined  by 
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which  is  a  conservative  estimate  of  the  maximum  value.  The  maximum  or  peak 
electric  field  is  determined  by  the  method  described  by  Adams  and  Mendelovicz 
(reference  13).  The  peak  electric  field  is 
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and  where  C^,  C  ,  and  are  the  phase  angles  for  the  corresponding  field 
component.  *  ^ 


The  near  magnetic  fields  are  computed  by  a  comparable  method.  As  is  the 
case  for  electric  near  fields,  the  currents  or.  the  wires  are  approximated  ty 
the  current  pulses  of  the  MININEC  solution.  A  virtual,  thin-wire  dipole  is 
placed  at  the  near  field  point  with  its  axis  parallel  to  the  appropriate 
vector  component.  The  near  magnetic  field  is  then  calculated  using  the 
MININEC  current  distribution  and  the  difference  between  the  appropriate 
components  of  the  vector  potential. 

The  vector  potential  is  generally  defined  such  that 
1 
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expressed  In  rectangular  coordinates,  becomes 


where  1,  J ,  k  are  the  unit  vectors  parallel  to  the  x,  y.  z  coordinate  #x:s, 
respectively.  And  where  A  .  A  .  A  are  the  cor  respond  l  ng  components  ...  f  tre 
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where  r  is  the  position  vector  at  the  observation  point,  k  r  r  /  |r  i.  and  K  = 

o  o 

tck  :  —  n  .  The  integral  is  evaluated  in  closed  form  over  each  straight  wire 

A 

segnent  for  each  current  pulse  and  is  reduced  to  a  susmation  over  the  wire 
segiaents.  The  fields  are  then  evaluated  as  real  and  i«aginary  parts  of  the  e 
and  «  e  opponents  at  a  specified  radiai  distance.  If  the  radial  distance  is 
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defined  by  circular  boundaries  concentric  about  the  origin  or  linear 
boundaries  parallel  to  the  y-axis  and  spaced  along  the  positive  x-axis.  Thus, 
in  the  latter  case,  the  ground  surface  is  divided  into  "strips"  at  user 
defined  x-axis  intercepts.  In  the  former  case,  the  ground  surface  is  divided 
into  concentric  rings  at  user  specified  radii.  In  this  case,  the  first  ring, 
or  zone,  may  include  a  radial  wire  ground  screen.  For  both  circular  and 
linear  zone  grounds,  each  zone  may  have  a  different  surface  impedance  and  each 
zone  may  have  a  different  height  ( Z-coordinate)  relative  to  the  first  zone. 
In  MININEC,  the  nunber  of  zones  is  limited  by  an  array  dimension  and  is 
currently  set  to  5. 

In  the  Fresnel  reflection  coefficient  method,  the  far  field  is  obtained 
by  summing  the  contributions  of  a  direct  ray  and  a  reflected  ray  from  each 
current  pulse.  The  field,  due  to  the  reflected  ray,  is  modified  by  the 
Fresnel  plane  wave  reflection  coefficient,  which  depends  on  the  ground  surface 
impedance  at  the  bounce  point,  or  specular  point,  and  the  angle  of  incidence. 

The  Fresnel  reflection  coefficients  have  not  been  applied  to  the  MININEC 
current  calculation.  When  a  real  ground  is  specified,  the  currents  are  cal¬ 
culated  by  using  the  perfectly  conducting  image  theory  (described  in  section 
2.3).  The  real  ground  corrections  are  applied  to  the  far  field  calculations 
only.  This  compromise  is  designed  to  keep  MININEC  relatively  compact  and 
provide  accurate  results  whenever  the  ground  directly  beneath  the  antenna  is  a 
good  conductor. 

Following  along  the  lines  of  the  development  given  by  Burke  and  Poggio 
(reference  4),  a  wave  incident  upon  a  finite  ground  (i.e.,  a  real  ground) 
yields  a  reflected  field,  Er,  given  by 

Er  =  FH[  ( Ej  •  .')  P]  ♦  Ry  [Ej  -  (Ej  •  P)  P]  (45) 


er  *  Vi  ♦  (rh  -  V  (EI  ’  P)  P 


where  P  is  a  unit  vector  perpendicular  to  the  plane  of  incidence,  E^  is  the 
incident  field  and  R  and  R„  are  the  vertical  and  horizontal  reflection 


jfficients,  respectively.  ^ 


The  two  terms  in  square  brackets  in  equation  (44)  correspond  to  hori¬ 
zontally  and  vertically  polarized  waves.  The  reflected  field  is  obtained  by 
decomposing  the  incident  field  into  horizontally  and  vertically  polarized 
waves,  computing  a  reflected  wave  for  each,  and  recombining  the  two. 

The  vertical  and  horizontal  coefficients  are 
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-  .s  trw  angle  of  Incidence  and  Z  is  the  relative  impedance  of  the 
r  sur  'ace  relative  to  the  free  space  impedance). 


To  obtain  the  far  fields,  tha  integral  in  equation  ( *0 )  iaplles  the  suaMtion 
over  all  the  current  pulaea.  For  the  direct  field,  the  currents  and  vectors 
pointing  to  the  current  pulse  centers  are  calculated  and  stored  in  arrays  by 
MINIVEC  during  the  aatri^  solution  process.  The  incident  field  on  the  ground 
surface  is  also  computed  froa  a  s venation  over  the  currents,  but  this  requires 
the  coordinates  of  the  specular  point  (the  bounce  point).  For  the  geoaetry 


illustrated  in 

Figure  5.  the 

specular  point  is  given  by 
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The  value  of  i  or  r  Is  used  appropriately  for  the  caae  of  linear 

or  circular  HH&u,'§&undar  ft<sin?S  determine  In  uhlch  aedla  the  bounce  occurs. 
The  height  of  the  ground  at  this  point  is  used  to  locate  the  usage  of  the 
source . 

The  ground  surface  inpedance  in  any  zone  is  given  by 


The  surface  iapedance,  when  a  ground  screen  is  present  and  the  specular  point 
lies  on  the  ground  screen,  is  given  by  Wait  in  reference  ’  *  also  see 
reference  * ) .  The  iapedance  of  the  ground  screen  by  itself  is 

jf  r 
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where  N  is  the  nutter  of  wires  in  the  ground  screen  and  C  is  the  ■'adius  f 
each  wire.  The  effective  ground  iapedsnee  is  foraed  by  coasting  the  peralle. 
iapedance  of  the  ground  without  the  ground  screen  and  the  iapedance  of  the 
ground  screen  without  the  ground,  or 
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( where  Z  = 


Z^  if  no  ground  screen  is  present 


The  total  field  at  e  point  s.e.r 
reflected  fields  as  described.  ^an  the 
gain  option  is  selected,  the  eJ 
field  is  used  in  equation  to 


is  the  vector  sue  of  the  d 
range,  *•  .  .s  set  *  c  zero  or 
r  tern  is  set  tc  unity.  The  total 
calculate  the  power  gain. 
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The  tv  ar  «nt«nra  prot.em  generated  Ty  «  method  >..f  »o««nt  5 

omputer  program  .s.  at  t*st  ,  ar  appr  01  imat  1 1  r .  Hc%  the  solution  is  tv 

'eaiity  depends  ir.  part  r  the  m»*rical  methods  employed  :n  the  code  and 

ho*  eel.  these  methods  are  implemented  .  the  inherent  accuracy  l.e.,  the 
njeter  -f  significant  digits  of  the  coaiputer  .  >  hoe  well  the  antenna  being 

modeled  i_rf\r»s  •  the  .imitations  l.e.,  simplifying  assumptions  of  the 

e .  ec  t  r  magnet  ;,-j  f  rrru.at  i»-r  ..see  tv.  reate  the  computer  program,  and  n  the 
.ser  ’  s  eiperier.ee  .  Nonethe.ess .  highly  accurate  answers  car  be  obtained  by 
a  r  e  f  u  .  m  c  d  e  .  .  r  g  ,.  f  the  antenna  ccnfiguraticn,  taking  intc  account  the 
.  nhe^ert  .  iir  1 1  at  ions  .  t  the  computer  program. 

6e.iah.e,  accurate  arswets  are  obtained  when  the  user  has  accumulated 
Sv.f!.  .eft  eiper.erce  !  r  <m  frequent  arc:  systematic  eiercise  of  the  program  tc 
-eti.  grije  pn-t.em  areas  He  must  he  fully  aware  of  potential  difficulties 

•  t  r  ugh  „  ..t  *  te  »vOf .  :  r  g  ;r  .  ess,  from  initially  setting  up  a  problem,  tv.  inter¬ 
net. ng  ‘he  "es-lts.  It  ;s  r  ec  ommended  that  the  user  run  h'lfclUEi  for  a  masher 
f  e.emertar,  ;r  I  .ns.  omparir.g  the  results  tc  independent  solutions  or  real 
uv  r  .  :  measurement s ,  ^r.t  ;  .  he  h  as  the  confidence  tc  apply  the  code  tc  a  problem 
‘  r  mt  .  t  •  ‘  e  answer  .s  ..nwrowr.  "The  examples  ir  this  section  provide  a  good 
;  .  e  •  .  ft  ar* 


eve!  peer  t  f  ,-orf  i dance  ir  t  h,e  computer  soluticr.  is  a  process  of  Sla¬ 
very  f  the  .  imitati  .rs  f  the  computer  program.  This  entails  the  modeling 
:  *  r  umh  er  f  simple  antenna  structures  found  in  standard  teats  on  antenna 

theory.  ’he  results  may  he  used  as  guidelines  tc  model  more  complex  antenna 

geometries  Natural  first  choices  are  dlpcles  in  free  space  and  over  ground 

r  sen  pc  1 e s  ,  followed  hy  TEE  and  inverted  1  shaped  antennas,  etc.  For  each 
ir’ffv  •  y  p  e  ,  a  number  .  f  prchlems  are  selected  involving  different  wire 
,erg»r.  Ar  .  -ad;..  “TN1NF  ;s  "or  fer  each  antenna  prctleir.  a  number  of  times. 

•  a*  .  :  rg  •  he  segment  at  lor  :.e..  a  crvergerce  test  and  other  parameters  tc 

*he  i r  gram  limitations.  Insight  for  effective  application  of  the 
P r  g'  ar  ga.-e:  f  -  at  om.par  i  scr  ?  with  measured  data  or  analytic  solutions 

w* er  aval  1  ah  .e  Ir  »► : s  manner ,  modeling  guidelines  are  derived  and  updated 
‘  ’  w  s  i*p  .er  jrtfrr  a  pr lems  for  wh:  h  there  is  reliahle  measured  data  cr 
.  )  a.  e  p  t  ed  •  h  e<  r  e  t  l  a  .  !»t  a  . 

•  ;r  :  f  AN'FNKA 

*  '  e  •  •  e;  f  e*  ;  a;  t  e  *  a  v  i  r  .  f  ♦  h  e  dipole  ant  er.na  has  h  e  e  r  studied 

.rter.se.  y.  ard  The  .iterature  is  rich  with  examples  that  include  measured  and 

•  *  ei.ret  i  a.  lata,  -efererves  r  .  -.  ‘ 1  .  empar  iscr  of  *IN1KE<  resul  *  s  to 

ha* a  f  r  dip.  ,es  f  va-i  us  dimensions  establishes  validation  and 

.  '  •  •  v  1  as;.1  ‘  r  a.  <;el  ;  -  g  guidelines.  F .  r  example.  onvergence  tests 

•  r  .  «  r  .  .e-g**  dipc.es  -  e  v  e  a  Is  'he  accuracy  that  ar;  he  expected  and 

.  -  .  :e  -  .i  •»*;.  r».  -  ,  *  e-  .  r  f  r  se  let  t  iv.r  t  segment  at  i  -.,r  der.s  1 1  y  t  he 

-  rter  f  segments  per  wire  h  ased  or  wire  .er.gth  ir  wave  .engths, 

h.gi.re'  ’  •  ►  r  ugh  shvw  'he  resul  t  s  of  rcr  ver  gerce  tests  fer  an  elec- 

'  •  a .  .  .  ah  r  *  - ;  p..<  ;  e  much  shorter  *  har  the  first  res.  nance  length  ,  a  dipcle 

•  ea»  *  e  n  r„r  e  .  -*rd  a  ).;<.►  rear  art  ires,  nance  ,  r  espe.  t  i  vel  y  .  Each  ,  dipole 

•-  <•  •  •  .  a  ,  '  *  ;  •  ard  er.  t  e  r  fr  i  ver  .  Far  ’  a  1  ;  gur  e  s  f  ,  b  ,  and  T  give 

•  •  e  .  ar  .  at  .  r  !  adm  ; » tar.,  e  .  .  er  sus  the  numt  er  f  segment  j?  .  Fart  t  Figures 
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cal.lv  shout  dipole  convergence  test 


Figure  7.  Convergence  test  for  an  electrically  short  dipole  showing  the  percent  difference  in  admittance  between  MININEC  and  R.W.P.  King 
(references  8  and  91  (Part  b). 
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■igure  8  Convergence  test  for  a  half  wave  dipole  when  admittance  is  given  (Part  a) 


DIPOLE  CONVERGENCE  TEST 


i.  Convergence  test  for  a  half  wave  dipole  showing  the  percent  difference  in  admittance  between  MININEC  and  R.W.P.  King  (references  8 
(Part  b) 


(ml  1 1 1  MHOS)  ANTIRESONANT  DIPOLE  CONVERGENCE  TEST  B  CmllllMHOS) 


(V 


O 

+  -  -  - 


(V 


- ► 


I 


(0 


QD 

I 


in 


u 


in 

h 

z 

U1 

z 

u 

Ui 

in 

u. 

a 

PC 

UI 

o 

z 

3 

z 


> 


m 

♦ 

o 

I 

>- 

v 

Ui 

u 

z 

< 


z 

Q 

< 


« 

\ 


% 

\ 


\ 


% 

\ 


\ 

% 

\ 

®  V 


01 

Q 

I 

Z 

E 

m 

(M 

Q 

• 

♦ 


\ 


\ 


/ 


y 


/ 


/ 


\ 


\ 


\ 


/ 


o  / 
/ 


\ 


\ 


X 


X 


z 

a 

< 

o 

z 


0. 

cr 


\ 


X 


x 


— - 

10 


-  * - —  - 

(V  CD 

*  m 


4~ 

* 

co 


co 


27 


(V 


CD 

C0 


* 

CO 


a 

co 


441 

to 

_  CM 

w 

t- 

z 

N 

• 

Ui 

X 

E 

E 

z 

X 

13 

in 

1 

Ui 

® 

* 

a 

m 

CM 

in 

• 

• 

® 

>- 

-  Cd 

a 

(NJ 

® 

a 

u. 

a 

m 

a 

m 

a 

CM 

1 

• 

• 

Ui 

L> 

a: 

1 

X 

i 

i 

Z 

CD 

Ui 

< 

® 

<*- 

X 

X 

0 

h- 

• 

z 

3 

Z 


CO 


(\i 


Figure  10.  Convergence  test  for  an  antiresonant  dipole  when  admittance  is  given  (Part  a). 
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Figure  1 1 .  Convergence  test  for  an  antiresonant  dipole  showing  the  percent  difference  in  admittance  between  MININEC  and  R.W.P.  King 
(references  8  and  9).  (Partb). 


(7,  9,  and  11)  give  the  percent  difference  between  MININEC  and  the  values 
published  by  R.VI.P.  King,  references  8  and  9, 

The  electrically  short  dipole.  Figures  6  and  7,  and  the  half  wave  dipole, 
Figures  8  and  9,  show  definite  signs  of  convergence  and  stability.  No  sign  of 
convergence  is  seen  for  the  antiresonant  dipole.  Figures  10  and  11.  The 
authors  also  have  seen  similar  convergence  problems  near  antiresonance  for 
other  method  of  moments  codes  (notably  NEC).  Figures  6  through  11  can  be  used 
as  guidance  for  selection  of  the  segmentation  density.  An  antenna  is  modeled 
as  a  collection  of  wires.  Each  wire  is  divided  into  a  number  of  short 
segments  selected  by  the  user.  The  number  of  segments  to  achieve  a  desirable 
confidence  level  can  be  based  on  the  results  of  Figures  6  through  11, 
depending  on  the  wire  length  in  wave  lengths.  Using  this  data  does  not 
guarantee  convergence  or  the  percent  accuracy  for  a  more  complicated  antenna, 
but  it  does  provide  a  starting  point.  Convergence  testing  is  always  advisable. 

Given  the  convergence  properties  of  MININEC,  how  well  does  it  predict 
dipole  properties?  Figure  12  is  a  comparison  of  a  single,  30-segment  MININEC 
model  to  the  theory  given  by  King  (reference  9).  Shown  is  the  admittance- 
versus-dipole  length  for  both  MININEC  and  King.  The  difference  is  less  than 
.5  millimho  for  most  of  the  range,  with  the  greatest  difference  of  about  1.5 
millimhos  in  the  susceptance  at  a  dipole  length  of  .6^  wave  length.  For 
longer  or  shorter  antennas,  the  user  is  advised  to  perform  suitable 
convergence  tests. 

The  accuracy  of  the  method  of  moments  solution  depends  also  on  meeting 
the  thin  wire  criterion.  To  illustrate,  Figure  13  shows  the  variation  of 
admittance  versus  the  wire  radius^  The  data,  given  by  King  (reference  9)  are 
also  shown  for  radii  between  10”^  and  10~^  wave  lengths.  The  segment  to 
radius  ratio,  4/a,  is  25  at  10  and  2.5  at  10  .  For  thicker  wires  than 
10  ,  the  thin  wire  criteria  is  not  achieved  and  the  results  are  expected  to 
be  not  as  go^d .  The  data  show  valid  behavior  for  thin  wires  with  A/a  >  2.5  or 
radii  of  10  wave  lengths  and  smaller. 

Numerical  problems  may  occur  in  the  solution  when  quantities  become  toe 
small  for  the  inherent  accuracy  of  the  computer.  An  example  is  the  erroneous 
results  that  can  occur  for  very  short  segments.  Figure  1*4  shows  the  results 
of  a  test  designed  to  identify  the  short  segment  limit.  Shown  is  the 
admittance-versus-dipole  length  in  wave  lengths  for  a  10-segment  dipole  in 
free  space.  The  conductance  and.  susceptance  displays  the  proper  behavior  for 
a  dipole  lengt£  greater  than  10  wave  lengths.  This^  corresponds  to  a  segment 
length  of  10  wave  lengths  and  longer.  Below  10  ,  the  conductance  oscil¬ 
lates  about  the  expected  values  as  the  segment  length  is  reduced,  and  at  times 
displays  negative,  non-physical  values.  A  change  from  single  precision  to 
double  precision  extends  the  validity  range  to  even  shorter  segments,  but 
significantly  increases  the  solution  time  beyond  acceptable  run  times  for  a 
16-bit  microcomputer.  For  MININEC,  on  a  16-bit  machine  in  single  precision, 
the  segment  length  should  always  be  greater  than  10  wsve  lengths. 

Antennas  are  often  constructed  of  wires  and  towers  or  other  conductors 
with  vastly  different  radii.  Even  simple  dipoles  may  have  tapered  elements. 

A  typical  MININEC  model  may  therefore  involve  the  connection  of  wires  with 


ANTENNA  BEHAVIOR  OF  MINI NEC 


HALF  LENGTH  (WAVE  LENGTHS) 


D I  POLE  RADIUS  (WAVE  LENGTHS) 


admittance  v«.  electrical  length 


tllPClE  LENGTH  (WAVE  LENGTH) 


large  step  changes  in  radii.  The  stepped-radius  wire  junction  has  been 
extensively  studied  by  Glisson  and  Wilton  (reference  17).  Figure  15  is  the 
stepped  wire  geometry  used  in  their  study.  They  adapted  a  body  of  revolution 
computer  code,  PEC,  to  solve  very  accurately  for  currents  and  charges  along 
the  stepped  wire  antenna.  The  results  were  compared  to  NEC  in  this  study. 


Figure  15  Geometry  for  the  stepped  radius  antenna  problem. 


rigures  16,  17,  18,  and  19  show  the  current  distribution  predicted  by  PEC,  NEC 

and  MININEC  for  radii  steps  of  1:1.25,  1:5,  1:10  and  1:100,  respectively.  The 

NEC  data  are  from  their  report.  The  MININEC  results  follow  the  PEC  data 

surprisingly  well  for  all  step  ratios.  (We  believe  the  difference  between  NEC 
and  PEC  may  be  an  error  in  the  data  in  the  Glisson  report.  We  have  not 
observed  this  difference  in  NEC  data.)  Further  investigation  of  MININEC  for 
different  stepped  radius  problems  should  be  conducted.  Suggestions  include 
moving  the  feed  closer  to  the  step  and  switching  the  radii  a1  and  a^. 

Multiple  wire  antenna  structures  may  often  require  very  close  spacing. 
When  the  spacing  is  very  small,  the  currents  may  not  be  adequately  represented 
by  a  thin  filament  on  the  wire  axis  as  it  is  represented  in  MININEC.  Figures 

2C,  21,  22,  and  23  show  MININEC  data  for  a  parallel  wire  test  used  to  investi¬ 
gate  the  close  spacing  limit.  The  test  consists  of  evaluating  the  self  and 

mutual  admittance  between  two  parallel  half  wave  dipoles.  One  antenna  is 
driven  (i.e.,  the  source  is  in  the  center  pulse)  while  the  second  is  not.  The 
self  admittance  is  the  feed  point  current  on  the  first  wire  if  the  applied 

voltage  is  1  ♦  jO  volts  and  the  mutual  admittance  is  the  current  for  the 

center  pulse  on  the  second  wire. 

Figure  20  shows  the  self  admittance  compared  to  the  theory  by  R.W.P.  King 
reference  9)  for  dipole  center  to  center  spacings  between  .1  and  .5  wave 
length.  Figure  21  shows  a  similar  comparison  for  the  mutual  admittance  over 
the  same  range.  The  differences  between  MININEC  and  R.W.P.  King  are  mostly 
less  than  .2  millimho  and  are  no  greater  than  .4  millimho  in  the  worst  case 
over  the  range  shown  for  both  self  and  mutual  admittance. 


CURRENT  DISTRIBUTION  OF  MlNINEC  COMPARED  TO  NEC  ANO  PEC 
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Figure  18  Currents  for  a  stepped  radius  (unction  of  aj/a 


Figure  19  Currents  for  a  stepped  radius  junction  of  aj^  '  '00 


Figure  22.  Magnitude  of  the  mutual  admittance  between  closely  spaced  parallel  dipoles. 


Given  the  good  agreement  with  theory  down  to  a  spacing  of  .1  wave 
length,  how  does  MININEC  fare  for  closer  spacing?  Figures  22  and  23  show  the 
magnitude  and  phase  of  the  mutual  admittance  for  spacings  down  to  the  point  of 
contact  of  the  two  parallel  dipoles.  Keep  in  mind  the  good  agreement  between 
MININEC  and  theory  for  spacings  of  .1  wave  length  and  greater.  (Reference  9 
does  not  provide  data  for  spacings  less  than  .1,  no  comparison  is  shown.)  It 
can  be  seen  that  the  magnitude  and  phase  continue  smoothly  as  the  spacing  is 
reduced.  Although  these  data  are  not  conclusive,  it  can  be  implied  that 
MININEC  can  model  antenna  configurations  with  wire  spacings  less  than  .1  wave 
length.  Whenever  a  model  has  close  spacing,  however,  it  is  advisable  to 
examine  the  results  very  closely  to  ensure  proper  behavior. 

3.2  LOOP  ANTENNAS 

A  circular  wire  loop  antenna  may  be  modeled  by  connecting  a  nunber  of 
wires  to  form  a  polygon  approximation  to  the  circular  loop.  A  simple  model 
has  one  segment  per  wire,  with  each  wire  forming  one  side  of  the  polygon 
model,  so  that  the  number  of  sides  and  the  number  of  segments  are  equal.  For 
a  given  circumference,  the  number  of  wires,  and  hence  the  number  of  segments, 
can  be  increased  until  the  solution  stabilizes,  indicating  the  number  of  sides 
required  to  model  the  circular  loop.  Figure  24  shows  the  results  of  this 
procedure  for  a  loop,  one  wave  length  in  circumference.  The  polygon  model  is 
circumscribed  by  a  circle  whose  circumference  is  one  wave  length.  The  wire 
radius  (a  =  .00674  meter)  is  chosen  to  correspond  to  the  published  data  given 
by  R.W.P.  King  (reference  9).  At  best,  the  real  part  of  the  MININEC 
admittance  comes  to  within  31  of  King's  data  and  the  imaginary  part  approaches 
to  within  61.  For  22  segments  (and  22  sides)  the  percent  difference  in  real 
and  imaginary  admittance  is  about  equal,  and  less  than  61  for  each. 

Figure  25  compares  MININEC  and  R.W.P.  King  admittance  data  for  a  range  of 
loop  diameters  from  .1  to  2.0  wave  lengths.  The  MININEC  model  is  the  22  seg¬ 
ment  or  22  sided  polygon  loop.  The  agreement  is  excellent.  The  difference 
between  King  and  MININEC  is  no  greater  than  .4  millimho  over  the  entire  range. 
From  .1  to  .8  wave  length,  the  MININEC  data  and  King  data  are  virtually 
identical . 

Figures  26  and  27  show  MININEC  data  for  small  loops  with  a  circunference 
from  10  to  just  above  .4  wave  length.  Keep  in  mind  the  excellent  agreement 
with  King's  data  for  loops  of  .1  and  greater  (Figure  25).  The  real  and 
imaginary  parts  of  the  admittance  Figures  26  and  27,  respectively,  are  well 
behaved  for  loops  greater  than  10~  wave  lengths.  Below  10  ,  the  real  part 
of  the  admittance  becomes  unstable  due  to  numerical  problems  encountered  at 
the  limits  of  single  precision.  Note  that  at  22  segments,  the  segment  size  at 
a  circumference  of  10~  ,  is  very  nearly  the  same  short  segment  length  limit 
displayed  by  the  dipole  test  in  Figure  14.  The  data  in  Figures  25,  26,  and  27 
suggest  a  smellioop  limit  for  MININEC  (on  a  16-bit,  single  precision  micro¬ 
computer)  of  10  wave  lengths  in  circumference.  This  corresponds  to  a  loop 
18  inches  in  diameter  at  2  MHz  (about  the  size  of  a  basketball  goal). 

3.3  M0N0P0LES  AND  ANTENNAS  ABOVE  GROUND 

Simply  stated,  an  antenna  above  a  perfectly  conducting  ground  plane  is 
equivalent  to  the  original  antenna  and  its  mirror  image  in  free  space.  Hence, 
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MIN  I NEC  LOOP  ANTENNA  BEHAVIOR 


CIRCUMFERENCE  (BETA-b)  IN  WAVE  LENGTHS 
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LOOP  CIRCUMFERENCE  (BET A-b)  IN  WAVE  LENGTHS 


•MlllMHOS  MINI NE C  LOW  FREQUENCY  BEHAVIOR 


\m*i!  'ooiw  prints  tMfi  by  MININf-f  (P^»* 


all  the  modeling  results  and  guidelines  presented  so  far  are  directly  applic¬ 
able  to  monopoles.  Specifically,  the  convergence  properties  illustrated  in 
Figures  6  through  11  can  be  used  for  the  initial  selection  of  the  segmentation 
required  for  a  monopole.  However,  this  should  not  preclude  convergence 
testing  whenever  possible. 

Figure  28  illustrates  the  geometry  of  a  TEE-antenna.  The  antenna  is 
driven  or  fed  at  its  base  from  a  coaxial  termination  at  the  ground  plane.  The 
dimensions  for  two  TEE-antenna  designs  (K  h  =  .2  and  K  h  r  .5)  are  also  given. 
A  convergence  test  was  performed  for  each  antenna  using  the  segmentation 
scheme  in  the  table.  The  results  of  these  tests  are  given  in  Figures  29  and 
30  for  K  h  =  .2  and  K  h  =  .5,  respectively.  A  comparison  of  the  "best" 
results  1 8  the  measurements  of  Prasad  and  King  (reference  18)  for  MININEC  and 
several  other  codes  is  given  in  Figure  31. 

Figure  31  compares  five  computer  programs  including  MININEC  for  the  two 
TEE-antennas .  In  each  case,  the  programs  were  tested  for  convergence  and  the 
best  answer  with  respect  to  Prasad's  measurements  is  given.  NEC  is  the  code 
previously  described  (reference  M).  TCP  (Tr iangular-Galerkin  Procedure)  is 
the  code  written  by  Chao  and  Straight  (reference  11)  using  triangular 
expansion  and  testing  functions  in  a  Galerkin  procedure  (i.e.,  triangles  for 
both  testing  and  expansion  functions).  PSRT  (Piece-wise  Sinusoidal  Reaction 
Technique)  is  a  sinusoidal  Galerkin  code  written  by  Richmond  (reference  19). 
TWTD  (Thin  Wire  Time  Domain)  is  a  time  domain  method  of  moments  code  written 
by  Van  Elaricun  and  Miller  (reference  20).  TWTD  uses  subsection  collocation 
(quadradic  interpolation  with  point  matching)  to  solve  for  the  time-dependent 
induced  currents  and  time-dependent  radiated  fields.  Admittance  data  are 
obtained  from  a  discrete  Fourier  transform  of  the  source  current.  All  codes 
except  MININEC  are  in  FORTRAN  and  require  mainframe  (large)  computers.  The 
data  show  that  MININEC  can  provide  equally  accurate  answers. 

3.**  NEAR  FIELDS 

MININEC  can  calculate  the  near  fields  for  antennas  in  free  space  and  over 
perfectly  conducting  ground.  Only  antennas  over  perfectly  conducting  ground 
are  considered  in  this  section. 

Figure  32  shows  a  comparison  of  MININEC  to  NEC  and  to  measurements.  The 
data  are  the  near  electric  fields  of  a  10.67-meter  monopole  over  a  good  con¬ 
ducting  ground  screen  at  2  MHz.  The  fields  are  for  1  KW  radiated  power.  The 
measurements  were  made  using  an  E-field  sensor  (EFS-1)  manufactured  by  Instru¬ 
ments  for  Industry  (reference  21).  The  NEC  data  are  from  a  single  precision 
version  (NEC-1)  running  on  a  VAX  11780  computer.  The  agreement  between  both 
codes  and  the  measurements  is  acceptable  over  the  range  shown.  The  accuracy 
of  the  measurements  is  5  to  10S,  with  the  greatest  error  occurring  for 
distances  of  10  meters  and  greater,  due  to  the  effects  of  nearby  structures. 
The  differences  between  NEC  and  MININEC  are  due  in  part  because  the  NEC  data 
is  from  a  single  precision  version. 

Figures  33  and  3*1  are  a  comparison  of  MININEC  and  NEC  near  field  data  for 
a  quarter  wave  monpole  over  perfect  ground.  The  NEC  data  are  from  a  double 
precision  version  (NEC-3)  running  on  a  VAX  11780  computer.  Figure  33  gives 
the  vertical  and  radial  components  of  the  electric  field  and  Figure  3^  gives 
the  phi-component  of  the  magnetic  field.  The  MININEC  data  has  been  scaled  to 
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TEE  ANTENNA  CONVERGENCE  TEST 


l.  Convergence  test  for  a  TE E-antenna  with  K( 
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Figure  31.  Comparison  of  TEE-antenna  impedance 
computations  with  the  measured  values  of  Parsad. 


the  power  level  of  the  NEC  data.  The  NEC  and  MININEC  data  are  essentially 
identical.  Figures  35  and  36  show  the  percent  difference  between  the  NEC  and 
MININEC  fields  of  Figures  33  and  34.  The  greatest  difference  occurs  very 
close  to  the  monopole  within  a  segment  length. 

If  the  near  fields  are  calculated  along  the  surface  of  the  monopole,  the 
differences  between  MININEC  and  NEC  are  much  more  pronounced.  Figures  37  and 

38  show  the  electric  fields  along  the  wire  surface  of  the  monopole  and  Figure 

39  shows  the  magnetic  fields.  The  MININEC  data  were  scaled  to  the  same  power 
radiated  as  the  NEC  data.  The  differences  are  due  to  the  approximation  used 
by  MININEC  to  determine  the  fields  since  the  current  distribution  (not  illus¬ 
trated)  is  nearly  the  same  for  both  codes.  The  impedance  calculated  by  each 
code  is  a  measure  of  this  agreement.  MININEC  predicts  an  impedance  of  42.170 
♦  J  21.478  ohms  and  NEC  predicts  42.387  J  24.873  ohms. 

The  accuracy  of  the  MININEC  near  fields  has  been  illustrated.  MININEC 
near  fields  are  sufficiently  accurate  for  well  converged  solutions  at  dis¬ 
tances  greater  than  a  segment  length. 

3.5  FAR  FIELDS 

The  correct  pattern  shape  can  often  times  be  calculated  using  a  coarse 
approximation  to  the  antenna  current  distribution.  However,  since  the  antenna 
input  impedance  is  used  to  determine  gain,  it  is  necessary  to  use  a  well- 
converged  solution  to  obtain  accurate  gain  data.  Figure  40  illustrates  both 
these  points.  Shown  is  a  comparison  between  MININEC  and  the  classical  solu¬ 
tions  from  Schelkunoff  (reference  22)  and  Jasik  (reference  23).  The  coarse 
solutions  are  represented  by  the  Schelkunoff  and  Jasik  data.  Their  data  are 
obtained  by  assuming  sinusoidal  currents  as  noted.  The  MININEC  data  are 
obtained  by  reference  to  the  convergence  data  of  the  previous  sections  and  by 
adjusting  the  frequency  to  obtain  the  exact  resonance  condition  of  near  zero 
reactance.  The  agreement  in  gain  and  impedance  data  (when  available)  is 
fairly  good.  Data  given  by  Schelkunoff  for  antennas  longer  than  one  wave 
length  cannot  be  trusted  because  of  the  assumptions  he  employs  for  the  current 
distribution. 
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Figure  33.  Near  electric  fields  of  a  quarterwave  monopole  computed  by  MININEC  and  NEC. 


Figure  34.  Nea.  magnetic  ♦••Ms  of  a  quartervvave  monopole  computed  by  MININEC  and  NEC 
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VERTICAL  ELECTRIC  FIELD  vs.  VERTICAL  DISTANCE 
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Figure  37.  Vertical  component  of  the  electric  field  at  one  radius  distance  for  a  quarterwave  monopole. 
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Figure  39.  Phi-component  of  the  megnetic  field  at  one  radius  distance  for  a  quarterwave  monopole 
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3.6  MEMORY,  DISK  STORAGE  AND  RUN  TIME 


Computer  memory,  disk  storage  capacity,  and  solution  time  are  key 
limiting  factors  in  the  use  of  all  method  of  moments  thin  wire  antenna  codes 
because  of  the  need  to  store  and  manipulate  (solve  or  invert)  a  matrix  of 
complex  numbers.  These  limits  are  particularly  acute  when  using  a 
microcomputer.  Mega-byte  hard  disks  and  compilers  can  partially  alleviate 
these  limits. 

MININEC  has  been  written  specifically  for  use  in  the  personal  computer 
environment.  Hence  the  choice  of  the  BASIC  language,  and  the  choice  of  the 
simple  pulse  expansion  and  testing  functions  (to  keep  overhead  down).  Every 
effort  has  been  made  to  produce  a  fast  compact  computer  code.  Earlier  ver¬ 
sions  of  MININEC  were  written  to  minimize  program  size  (length).  The  present 
version,  however,  has  sacrificed  size  for  improved  internal  documentation, 
modularity  and  increased  capability. 

Figure  is  a  comparison  of  the  run  times  between  MININEC  version  2  and 
MININEC  version  3.  Both  codes  were  compiled  for  comparison  using  a  Microsoft 
BASIC  compiler.  Use  of  a  math  co-processor  will  significantly  reduce  the  run 
times.  The  co-processor  was  not  used  to  obtain  the  data  in  Figure  41.  For 
comparison,  some  other  attributes  of  the  codes  are  as  follows: 


MININEC(2)  MININEC(3) 


INTERPRETER 

COMPILED 

INTERPRETER 

COMPILED 

No.  of  lines 

543 

— 

1607 

— 

Max.  no.  of  wires 

10 

75 

10 

50 

Max.  no.  of 

segmt  nts 

50 

75 

30 

50 

Disk  storage 

(k  bytes) 

13 

57 

44 

108 

The  solution  time  and 

size  of  the  executable  is  a 

function  of  the 

compiler 

the  compiler/linker  options  used. 

The  significant  increase  in  speed  of  MININEC(3)  is  attributable  to  ar. 
improved  solution  routine.  The  matrix  fill  time,  the  time  to  compute  all 
terms  of  the  matrix,  is  virtually  the  same  for  both  versions  of  MININEC.  For 
large  problems,  the  fill  time  is  usually  longer  than  the  factor  time,  the  time 

to  solve  the  matrix.  Figure  42  co.., pares  the  matrix  fill  time,  factor  time  ami 

total  solution  time  for  MININEC(3)  to  solve  a  dipole  in  free  space.  Matrix 
fill  time  dominates  for  problems  above  20  segments. 

The  solution  time  not  only  increases  with  the  number  of  segments,  but 
also  increases  with  the  number  of  wire  junctions.  This  effect  can  be  seen  it; 
the  data  in  Figure  43.  Shown  is  the  total  solution  time  (fill  time  plus 

factor  time)  for  a  half  wave  dipole  and  a  loop  antenna  in  free  space.  The 
loop  is  an  extreme  case  in  which  the  number  of  wire  junctions  equals  the 
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Figure  41 .  A  comparison  of  the  run  times  of  MININEC  version  2  and  3. 
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43  Solution  t  me  tc  VlNINECi3l  to  *olve  a  dipole  loop  antenna,  monopole  end  TEE  antenna 


number  of  segments.  Also  shown  in  Figure  U3  is  the  total  solution  time  for  a 
quarter  wave  monopole  and  a  TEE-antenna  (previously  described)  over  perfect 
ground.  The  wire  junction  in  the  TEE-antenna  is  responsible  for  the  longer 
solution  time  compared  to  the  monopole.  The  solution  time  is  also  longer  for 
antennas  over  perfect  ground  compared  to  free  space.  For  comparison,  the 
monopole  is  equivalent  to  a  dipole  in  free  space  with  twice  the  nunber  of 
segments.  Even  so,  the  dipole  solution  time  is  a  little  shorter  than  the 
monopole. 

The  brute  force  approach  to  circunvent  memory  and  solution  time  limits  is 
to  seek  out  ever  larger,  faster  computers.  The  logical  extension  is  to 
rewrite  MININEC  in  a  more  powerful  language  such  as  FORTRAN  and  use  a  main¬ 
frame  computer.  But  if  a  mainframe  is  available,  one  really  should  be  using 
the  more  powerful  NEC  program.  NEC  is  written  in  FORTRAN  and  is  designed  for 
efficient  mainframe  use.  The  concept  of  MININEC  is  to  use  a  personal  computer 
(PC).  In  time,  as  PCs  become  more  powerful,  then  MININEC,  too,  will  expand  in 
capability.  In  the  meantime,  MININEC  is  appropriate  for  application  to  small 
problems  (less  than  100  segments).  Larger  problems  should  be  solved  with  NEC 
on  an  appropriate  size  computer. 


4.0  EXAMPLES  AND  USER  GUIDANCE 


This  section  is  intended  to  be  used  both  as  a  reference  and  as  a  means 
for  first-time  users  to  become  familiar  with  the  input  and  output  (I/O) 
options  of  MININEC.  The  first-time  user  should  spend  a  few  minutes  at  the 
computer  terminal  while  following  along  with  the  simple  two-wire  example  in 
section  4.1.  A  few  more  minutes  at  the  terminal  (perhaps  with  a  simpler  four 
segment  dipole)  exploring  all  of  the  MININEC  options  should  be  enough  to 
master  this  skill.  However,  the  art  of  actually  modeling  a  wire  antenna, 
i.e.,  composing  the  model  and  properly  interpreting  MININEC  output,  is 
acquired  through  considerable  study  of  antenna  theory  and  properties,  study  of 
the  data  of  section  3*0,  and  equally  important,  accumulating  experience  by 
using  MININEC. 

4.1  GETTING  STARTED 

First,  gather  up  all  the  known  information  on  the  antenna  to  be  modeled, 
including  measurements  or  reliable  analytical  data,  if  available.  It  is 
helpful  to  make  a  sketch  of  the  antenna  using  any  convenient  cartesian 
coordinate  system.  For  this  example,  please  consider  the  inverted  L-antenna 
in  Figure  41.  You  will  need  to  know  the  X,  Y,  Z  location  for  each  wire  end  of 
each  wire  relative  to  the  origin  of  your  choice.  And  you  will  need  to  know 
the  radius  of  each  wire.  All  dimensions  must  be  in  meters.  This  information 
for  the  example  can  be  found  in  Figure  44. 

You  will  need  to  decide  how  many  segments  to  use  on  each  wire.  This  may 
be  done  initially  by  calculating  the  length  of  each  wire  in  wave  lengths  at 
the  desired  frequency.  Then  refer  to  Figures  6  through  11  as  appropriate  for 
the  initial  choice.  This  will  not  guarantee  that  the  MININEC  solution  will  be 
well  converged,  but  it  provides  a  good  place  to  start.  A  convergence  test  is 
always  a  good  idea.  A  possible  segmentation  scheme  for  the  inverted  L-antenna 
is  suggested  in  Figure  44.  Alternatively,  using  your  prior  experience  on  a 
similar  antenna,  you  may  be  able  to  come  closer  to  the  converged  solution  the 
first  time.  For  example,  the  inverted  L-antenna  is  one  wire  short  of  being  a 
TEE-antenna  similar  to  the  ones  in  section  3.3  (see  Figures  2 8  through  31). 
For  the  purpose  of  this  discussion,  however,  we  will  use  a  minimal  number  of 
segments  in  order  to  keep  the  solution  time  reasonably  short. 

MININEC  is  designed  for  demand  mode  execution.  Therefore,  all  you  need 
to  do  is  answer  the  questions  when  prompted  or  provide  the  required  data  appro¬ 
priately.  In  general,  you  must  define  the  antenna  geometry  and  define  the 
environment.  MININEC  can  then  solve  for  the  currents.  Once  MININEC  has  the 
solution,  you  may  then  request  near  and  far  fields  (patterns). 

Now  run  MININEC.  The  first  question  you  must  answer  is: 

OUTPUT  TO  CONSOLE,  PRINTER,  OR  DISK  (C/P/D)? 

If  you  enter  C  for  console,  all  data  will  be  displayed  on  the  monitor.  If  you 
answer  P  for  printer,  the  questions  or  prompts  will  still  be  displayed  on  the 
monitor,  hut  all  other  data  will  be  printed.  If  you  answer  D  for  disk,  you 
will  then  he  prompted  to  supply  the  name  of  a  disk  file  for  storage  of  all 
MININEC  output.: 
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WIRE  NO 

WIRE  RADIUS 
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Figure  44  Inverted  L  antenna  data  for  the  example  m  Section  4  1 
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OUTPUT  TO  CONSOLE,  PRINTER,  OR  DISK  (C/P/D)?  D 
FILENAMES  (NAME. OUT)? 

As  with  the  printer,  all  questions  and  prompts  appear  on  the  monitor,  but  the 
results  go  to  the  file  you  specify. 

Next  you  must  supply  the  frequency  in  MHz.  Enter  the  desired  frequency 
or  simply  enter  return  to  select  the  default  value  of  299.8  MHz  (or  a  wave 
length  of  one  meter). 

FREQUENCY  (MHz)? 

WAVE  LENGTH  =  1  METER 

For  the  next  question,  you  must  select  free  space  or  a  ground  plane  by 
entering  1  or  -1,  respectively. 

ENVIRONMENT  (-*-1  FOR  FREE  SPACE,  -1  FOR  GROUND  PLANE)? 

If  you  select  >1  for  ground,  you  are  then  prompted  for  the  number  of  media, 
which  must  be  an  integer  from  zero  to  5. 

ENVIRONMENT  (+1  FOR  FREE  SPACE,  -1  FOR  GROUND  PLANE)?  -1 
NUMBER  OF  MEDIA  (0  FOR  PERFECTLY  CONDUCTING  GROUND)? 

A  zero  selects  a  perfectly  conducting  ground  plane.  See  Section  4.3  for 
further  details.  Suffice  to  say  at  this  point  that  selection  of  1  to  5  media 
does  not  effect  the  current  distribution  or  the  antenna  impedance,  but  only 
effects  far  field  calculations.  Please  select  a  zero  to  continue  along  with 
this  narrative. 

« 

Next,  you  must  specify  the  number  of  wires  .  If  you  specify  more  than 
the  number  allowed,  a  warning  message  appears  and  the  question  is  repeated. 

NO.  OF  WIRES?  100 

NUMBER  OF  WIRES  EXCEEDS  DIMENSION.  .  . 

NO.  OF  WIRES? 

A  zero  will  place  you  at  the  main  menu. 


•  Note  to  NEC  users: 

The  MININEC  post  processor  program  is  used  to  convert  a  NEC  input  data  set 
Into  a  MININEC  antenna  geometry  description.  These  data  are  automatically 
stored  in  the  MININEC. INP  file.  At  this  point,  if  the  MININEC. INP  file  is  not 
empty,  MININEC  will  use  that  data  for  the  geometry  description  and  you  will 
skip  the  normal  geometry  input.  See  Appendix  A  for  further  information. 


6? 


NO.  OF  HIRES?  0 


iHHHiiiimiHHii  MININEC  MENU  ••••«•••••»••••••••••••■ 

G  -  CHANGE  GEOMETRY  C  -  COMPUTE/DISPLAY  CURRENTS 

E  -  CHANGE  ENVIRONMENT  P  -  COMPUTE  FAR-FIELD  PATTERNS 

X  -  CHANGE  EXCITATION  N  -  COMPUTE  NEAR-FIELDS 

L  -  CHANGE  LOADS/NETS  Q  -  QUIT 

F  -  CHANGE  FREQUENCY 


COMMAND? 

To  recover  from  this  point,  select  G  for  change  geometry  and  answer  the 
questions  on  environment  again. 

NO.  OF  HIRES?  0 


G  -  CHANGE  GEOMETRY 
E  -  CHANGE  ENVIRONMENT 
X  -  CHANGE  EXCITATION 
L  -  CHANGE  LOADS 
F  -  CHANGE  FREQUENCY 


COMMAND?  G 

ENVIRONMENT  (*1  FOR  FREE  SPACE,  -1  FOR  GROUND  PLANE)?  -1 
NUMBER  OF  MEDIA  (0  FOR  PERFECTLY  CONDUCTING  GROUND)?  0 

NO.  OF  HIRES? 

Let's  assume  two  wires.  The  next  prompt  is  for  the  number  of  segments  on 
wire  t. 

HIRE  NO.  1 

NO.  OF  SEGMENTS? 

You  will  be  prompted  for  each  wire  in  turn.  If  you  answer  zero,  you  will 
return  to  the  question  for  the  number  of  wires.  This  is  a  convenient  escape 
mechanism,  sometimes  useful  when  you  change  your  mind. 

HIRE  NO.  1 

NO.  OF  SEGMENTS?  0 

NO.  OF  HIRES?  2 

HIRE  NO.  1 

NO.  OF  SEGMENTS? 

If  at  any  time  you  specify  too  many  segments  on  any  one  wire,  or  the  total 
nuaber  of  aegments  on  all  wires  specified  becomes  larger  than  the  maximum 
allowed  by  the  program  array  dimensions,  an  error  message  will  be  displayed, 
and  you  will  return  to  the  question  for  the  number  of  wires.  In  this  case,  to 
keep  this  session  short,  let's  choose  4  segments  for  wire  1. 


MININEC  MENU  eeeeeeseeeeesseeeeeee 

C  -  COMPUTE/DISPLAY  CURRENTS 
P  -  COMPUTE  FAR-FIELD  PATTERNS 
N  -  COMPUTE  NEAR-FIELDS 
Q  -  QUITE 


Next,  enter  the  X,  Y,  Z  coordinates,  in  meters,  for  end  one  of  the  first 
wire.  Then  enter  the  X,  Y,  Z  coordinates  for  end  two;  then  enter  the  radius, 
as  prompted. 

NO.  OF  SEGMENTS?  4 

END  ONE  COORDINATES  (X,Y,Z)?  0,0,0 

END  TWO  COORDINATES  (X,Y,Z)?  0,0,. 191 

RADIUS?  .004 


COORDINATES 

END 

NO.  OF 

X  Y 

Z 

RADIUS 

CONNECTION 

SEGMENTS 

0  0 

0 

-1 

0  0 

.091 

.004 

0 

4 

CHANGE  WIRE  NO. 

1  (Y/N)? 

After  entry  of  the  radius,  MININEC  responds  immediately  with  a  table  as  shown 
above.  The  table  gives  the  coordinates  of  the  wire  ends,  the  radius  and  the 
end  connection  information,  so  you  may  verify  that  you  have  entered  the  data 
correctly.  If  not,  you  have  the  opportunity  to  start  again  with  this  wire. 
Otherwise,  you  may  continue  to  the  next  wire.  The  end  connection  information 
is  useful  to  verify  that  a  connection  has  been  made.  In  this  case,  the  minus 
one  indicates  that  end  one  is  connected  to  ground.  A  connection  to  ground  is 
indicated  by  a  negative  integer,  whose  absolute  value  is  the  same  as  the  wire 
number.  The  zero  indicates  that  end  two  of  this  wire  is  not  yet  connected  to 
any  other  wire. 

Now,  continue  on  to  the  next  wire  and  give  the  appropriate  data: 

CHANGE  WIRE  NO.  1  (Y/N)?  N 

WIRE  NO.  2 

NO.  OF  SEGMENTS?  6 

END  ONE  COORDINATES  (X.Y.Z)?  0,0,. 191 

END  TWO  COORDINATES  (X.Y.Z)?  0,. 309,. 091 

RADIUS?  .004 


COORDINATES 

END 

NO.  OF 

X 

0 

Y  Z 

0  .191 

RADIUS 

CONNECTION 

1 

SEGMENTS 

0 

CHANGE 

.309  .191 

WIRE  NO.  2  (Y/N)? 

.004 

0 

6 

The  connection  data  indicates  that  end  one  of  wire  two  is  connected  to  the  top 
of  wire  one,  i.e.,  end  two  of  wire  one.  The  absolute  value  of  the  end  con¬ 
nection  Integer  is  the  wire  number  to  which  wire  two  Is  connected.  A  plus 
sign  Indicates  an  end  one  connected  to  an  end  two  (as  In  this  case),  or  an  end 
two  connected  to  an  end  one.  A  negative  sign  indicates  an  end  one  connected 
to  an  end  one,  or  an  end  two  connected  to  an  end  two.  And,  of  course,  a  zero 
means  no  connection. 

By  the  way,  if  you  happen  to  give  either  wire  end  a  negative  Z-coordinate 
when  a  ground  plane  is  specified,  you  will  get  an  error  message  and  will  have 
to  reconsider  your  entry.  Likewise,  MININEC  will  not  let  you  get  away  with  a 
zero  radius,  or  a  zero  wire  length. 


If  you  are  now  satisfied  with  the  wire  two  data  entry,  we  may  proceed. 
MININEC  will  produce  a  table  of  the  coordinates  for  the  location  of  each 
current  pulse  on  each  wire.  Note  that  the  radius  and  connection  data  are  also 
given  so  that  you  may  verify  your  data  entry. 

CHANGE  WIRE  NO.  2  (Y/N)?  N 

***•  ANTENNA  GEOMETRY 


WIRE  NO.  1  COORDINATES  CONNECTION  PULSE 


X 

Y 

Z 

RADIUS 

END 

1  END  2 

NO. 

0 

0 

0 

.004 

-1 

1 

1 

0 

0 

.04775 

.004 

1 

1 

2 

0 

0 

.0955 

.004 

1 

1 

3 

0 

0 

.14325 

.004 

1 

0 

4 

WIRE  NO 
X 

.  2  COORDINATES 

Y 

Z 

RADIUS 

CONNECTION 
END  1  END  2 

PULSE 

NO. 

0 

0 

.191 

.004 

1 

2 

5 

0 

.0515 

.191 

.004 

2 

2 

6 

0 

.103 

.191 

.004 

2 

2 

7 

0 

.1545 

.191 

.004 

2 

2 

8 

0 

.206 

.191 

.004 

2 

2 

9 

0 

.2575 

.191 

.004 

2 

0 

10 

CHANGE  GEOMETRY  (Y/N)? 

This  table  is  essential  to  proper 

location 

of  the 

feed 

point 

(or  source 

excitation  point) 

and  location  of  loads. 

You  now  have  one  last  chance  to  change  the  geometry  before  proceeding. 
CHANGE  GEOMETRY  (Y/N)?  N 


NO.  OF  SOURCES? 

You  must  now  decide  how  many  feed  points  to  use,  where  they  are  located,  and 
what  voltages  are  applied.  You  will  be  prompted  for  this  data  for  each 
source,  in  turn.  In  this  case,  let’s  keep  it  simple. 

NO.  OF  SOURCES?  1 

SOURCE  NO.  1: 

PULSE  NO.,  VOLTAGE  MAGNITUDE,  PHASE  (DEGREES)?  1,1,0 

Sources  are  always  co-located  with  current  pulse  functions.  Hence,  you  may 
have  to  refer  to  the  above  table  of  antenna  geometry  to  select  an  appropriate 
pulse  to  apply  the  source.  In  this  example,  the  source  is  at  the  ground 
plane,  i.e.,  pulse  number  one,  located  on  wire  one.  I  have  chosen  one  volt  at 
zero-degree  phase  angle  for  this  example. 


The  next  set  of  questions  provide  the  opportunity  to  add  impedance  load¬ 
ing  to  the  antenna.  To  keep  things  simple,  let's  avoid  loading  and  continue 
on.  Please  refer  to  section  4.5  for  more  detailed  information  on  loading 
options. 

NUMBER  OF  LOADS?  0 


G  -  CHANGE  GEOMETRY 
E  -  CHANGE  ENVIRONMENT 
X  -  CHANGE  EXCITATION 
L  -  CHANGE  LOADS/NETS 
F  -  CHANGE  FREQUENCY 


MININEC  MENU  eeaeeseasaaaaaeaaaaae 

C  -  COMPUTE/DISPLAY  CURRENTS 
P  -  COMPUTE  FAR-FIELD  PATTERNS 
N  -  COMPUTE  NEAR-FIELDS 
Q  -  QUIT 


COMMAND  ? 


You  are  now  faced  with  the  main  menu  again.  By  selection  of  an  appro¬ 
priate  command  letter,  you  may  change  the  geometry  (G) ,  change  the  environment 
(E),  change  the  source  or  excitation  (E),  change  the  loads  and  networks  (L), 
and  change  the  frequency.  When  satisfied  with  the  antenna  geometry  and  envi¬ 
ronment,  you  are  ready  to  determine  the  antenna  properties.  The  preferred 
choice  at  this  point  is  C,  compute  and  display  the  currents.  If  you  select  P 
for  patterns  or  N  for  near  fields,  MININEC  will  check  to  see  if  the  currents 
have  been  calculated;  if  not,  MININEC  will  compute  the  currents  before  you  can 
proceed.  So,  let's  select  C. 


COMMAND?  C 


BEGIN  MATRIX  FILL 

MATRIX  FILL  10*  COMPLETE  -  APPROX  TIME  REMAINING  1:48 

MININEC  responds  almost  immediately  with  an  estimate  of  the  time  in  minutes 
and  seconds  required  to  complete  filling  the  impedance  matrix.  At  intervals, 
the  estimate  will  be  updated  and  the  total  time  will  be  given  when  this  step 
is  completed.  Similarly,  MININEC  will  estimate  the  time  to  solve  the  matrix 
for  the  currents  and  in  turn  display  the  total  times. 

BEGIN  MATRIX  FILL 

MATRIX  FILL  100*  COMPLETE  -  APPROX  TIME  REMAINING  0:00 
FILL  MATRIX:  1:27 


FACTOR  MATRIX  100*  COMPLETE  -  APPROX  TIME  REMAINING  0:00 
FACTOR  MATRIX:  0.04 


When  the  solution  is  complete,  MININEC  computes  and  displays  the 
impedance  and  power  input  for  each  source  in  turn.  If  there  is  more  than  one 
source,  the  sum  total  power  input  will  also  be  displayed.  The  solution  is 
also  displayed  in  terms  of  the  current  distribution,  wire  by  wire. 


tiiimitiiHHUHi  SOURCE  DATA  ••••••••••••*••• 

PULSE  1  VOLTAGE  =  (1  .  C  J) 

CURRENT  =  (  9 . 85227 8E-04  .  i .  479977E-C.3  J' 
IMPEDANCE  =  (  311.6818  .  -468. '982  J) 

POWER  x  4 . 926  1 39E-04  WATTS 


CURRENT  DATA 


WIRE  NO. 

PULSE 

NC. 

1 

-» 

c 

3 

4 

J 


1  : 

REAL 

(AMPS) 

9.852278E-C4 
.00C962 
6.940962E-C4 
7 .86','T?9E-C4 
6.45'’  1SE-C4 


IMAGINARY 

(AMPS) 

1 .4?99't?E-C? 
.880544E-C4 
-2. 186 38 3E-C 3 
-3.248648E-C? 
-3.943692E-C3 


MAGNITUDE 
, AMPS' 

1  .•’'”,922E-C  - 
’  .24  A*’  3E-2 ' 
2. 362  ’4 jt— C  : 
3. ; 42^6  3E-0 3 
; .  9962Cc'E-2  • 


PHASE 
.  DEGREES 

kt .  ''46  ’ 

-  <9  .  32  "r 
-6  '  .  ’Ch" 


-c  - 


WIRE  NC.  2: 
PULSE  REAL 


NC. 

J 

6 

T 

8 

9 

*C 


, AMPS) 

6.45"  15E-C4 
=>.  160325E-04 
3Hc'68C5E-24 
2.4C80*^E-C4 
1 . 2522 i 8E-C4 
4. 12’ 46E-C* 


E 


IMAG INARY 
. AMPS' 

-  3 . 94  36926  --  • 
-4 . 466  36 4E-0  ' 
-4.c.198c'9E-C  ? 
-4.C93054E-C 

-  3 . 2C 1 494E-C  3 

-  ’  .991"'  14E-2  - 


«AGN  ITUDE 
AMPS' 

.  996  2  -'L-  i  - 
4 . 4  9t'29‘'  E  -2  - 
4 . 4  ;<4446E-2  ; 
- .  '22 '  32E-2  • 
3.2C3942E-2  ; 
' .89216  3F-2 ? 


PHASE 

DEGREES 

_  4  '  ' 

-  S  '  .  4  ^  4  4 

-6k .  4862 
-86.*'  !  ; 


-6  6  . 
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SAVE  CURRENTS  TC  A  -ICE  Y  S 
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current  law  has  heer  _sed  to  solve  fo*"  "e  •>  ••  ♦  •  e  s  •.  e  *'->•"  *  -  r 

the  ’unction.  And  a  J  is  substitute?  ’  c  •  *-e  .  -  •  .ri  c  ,  <’.•>  <  *  ' 
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surface  (or  media)  relative  to  the  first  surface.  A  negative  value  is  a  step 
down;  a  positive  value  is  a  step  up.  A  zero  value  means  a  flat  ground.  Note 
that  there  is  no  diffraction  coefficient  correction  in  MININEC  for  a  cliff 
edge.  There  is  also  no  correction  for  blockage  due  to  a  large  step  up.  i.e., 
the  antenna  cannot  be  in  a  well.  For  each  media  in  turn,  you  must  supply  the 
distance  in  meters  from  the  origin  to  the  media  boundary. 

When  you  choose  a  circular  boundary  with  two  or  more  media,  you  will  be 
prompted  for  the  number  of  radial  wires  in  the  ground  screen.  This  approxi¬ 
mation  is  really  accurate  for  dense  ground  screens  only.  One  hundred  or  more 
should  be  used  for  best  results.  Zero,  of  course,  is  also  a  valid  choice. 
You  must  also  supply  the  wire  radius  of  the  wires  in  the  ground  screen.  The 
length  of  the  wires  in  the  ground  screen  is  the  same  as  the  distance  to  the 
interface  between  media  one  and  two. 

CHANGE  EXCITATION 

The  X  option  on  the  MININEC  menu  provides  the  means  to  change  the  antenna 
feed  excitation  (i.e.,  number  of  feeds,  feed  location  and  magnitude  and  phase 
without  changing  the  geometry,  environment,  frequency  or  loading.  Changing 
only  the  excitation  does  not  require  re-filling  or  re-factoring  of  the  matrix. 
Hence,  after  solving  for  the  currents  for  the  initial  excitation,  you  car. 
rapidly  try  out  all  possible  source  locations  or  try  multiple  source  excita¬ 
tion.  Source  locations  must  coincide  with  the  locations  of  the  current  pulse 
functions.  You  cannot  put  source  excitations  on  non-existent  pulses.  At 
least  one  source  must  be  giver.  MININEC  will  not  allow  more  sources  than  the 
maximun  set  by  the  dimension  statement  (see  the  first  30  lines  of  MININEC  tc 
determine  this  limit). 

There  is  but  one  source  model  in  MININEC.  The  model  imposes  a  constant 
field  over  a  pulse  width,  with  amplitude  and  location,  coincident  with  a 
current  pulse,  chosen  by  the  user.  In  spite  of  the  simplicity,  this  source 
model  is  a  good  approximation  for  most  transmit  cases.  What  do  you  do  then  tc 
evaluate  the  performance  of  a  receive  only  antenna?  One  way  to  evaluate  the 
receive  properties  of  an  antenna  is  to  place  a  transmitting  dipole  at  a  great 
distance  (i.e.,  in  the  far  field  of  the  receive  antenna).  The  field  incident 
on  the  receive  antenna  is  a  good  approximation  to  a  plane  wave. 

CHANGE  LOADS 

The  L  option  on  the  MININEC  menu  provides  the  means  to  alter  or  aac 
lumped  parameter  loads  to  an  antenna.  Changing  the  load  requires  re-filling 
and  re-factoring  of  the  matrix.  Each  load  location  must  coincide  with  a 
current  pulse  expansion  function.  MININEC  will  not  allow  more  loads  than  the 
maximum  for  which  it  is  dimensioned  (see  the  first  30  lines  of  MININEC  u 
determine  this  limit).  There  are  two  kinds  of  loads;  impedance  loading  anc 
S-parameter  loading.  You  may  select  the  loading  type  when  prompted  afte- 
specifying  the  number  of  loads. 

The  simplest  load  type  is  impedance  loading.  You  must  supply  the  pulse 
location,  resistance  and  reactance  for  each  load,  in  turn.  The  reactarnr 
value  you  specify  will  not  change  when  you  change  the  frequency.  If  you 


*-■  merit**  •  -hange  appropriately  with  frequency,  you  must 
mem-m-  •  f  f«»r;  mm  you  change  frequency.  Alternatively,  you 
t  ■  ’»rmm  of  an  equivalent  S-paraswter  function. 

-•  •  *»**<«i  f  circuit  analysis  aakes  use  of  the  concept  of  a 
•  -parameter  Generally,  S  is  defined  as  S  =  G  ♦  ju 

•*  •  •  aquarv  y  arc  „  is  the  angular  frequency  (reference  24). 

;•  '-me  varying  signals,  G  =  0. 

.  *  «  •*  »i  i  *'<  u 1 1  can  always  be  expressed  as  a  function  of 

*  a.gabra.  the  iapedance  can  be  represented  as  a  ratio 
■«  *  -  per  mm  ter  iapedance  function  is  a  polynomial  in  S 
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4.6  CHANGE  FREQUENCY 


The  F  option  on  the  MININEC  Menu  provides  a  way  to  alter  the  Frequency. 
You  will  be  proapted  for  the  frequency  in  MHz.  MININEC  will  coapute  and  dis¬ 
play  the  wave  length  and  return  to  the  aenu.  Changing  the  frequency  will 
require  re-filling  and  re-factoring  of  the  lapedance  aatrli.  The  current 
geoaetry  will  be  used. 

4.7  COMPUTE/DISPLAY  CURRENTS 

The  C  option  on  the  MININEC  aenu  triggers  filling  and  factoring  of  the 
aatrix  for  the  antenna  configuration  most  recently  specified.  If  a  solution 
has  already  been  coaputed,  and  no  changes  have  been  aade  to  the  envlronaent. 
excitation,  loads  and  frequency,  then  the  C  option  will  slaply  display  the 
iapedance  and  current  distribution. 

4.8  COMPUTE  FAR  FIELD  PATTERNS 

The  P  option  on  the  MININEC  menu  is  used  to  specify  the  far  field  pattern 
calculation.  If  the  currents  have  not  already  been  coaputed,  they  will  be 
coaputed  before  you  can  proceed.  You  may  choose  to  coapute  the  patterns  in 
dBI,  i.e.,  in  dB  above  an  isotropic  radiator,  or  in  volts  per  meter,  i.e.,  the 
electric  field. 

When  you  choose  dBI,  you  are  prompted  for  the  zenith  angle  and  the 
azimuth  angle.  In  each  case  you  must  supply  three  numbers,  the  initial  angle, 
the  angle  increment  and  the  number  of  angles.  If  you  specify  zero  for  the 
number  of  angles,  one  is  assumed.  The  zenith  and  azimuth  angles  are  the  theta 
and  phi  angles,  respectively,  as  shown  in  Figure  45.  When  the  patterns  are 

calculated  in  dBI,  the  e-^k^/R  dependence  of  the  far  field  is  suppressed, 
i.e.,  the  pattern  is  at  infinite  range. 

z 


Figure  45.  Coordinate  system  for 
far  field  patterns. 


When  you  choose  to  compute  the  patterns  in  volts  per  meter  (i.e.,  the 
electric  field  strength),  the  power  (radiated)  level  is  displayed  and  you  are 
proapted  to  change  this  level.  The  fields  will  be  scaled  to  the  level  you 
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specify.  »ut,  you  auat  specify  the  rang*  tn  a*t*r*.  If  you  *p*cify  zero, 

-  lit  r 

tha  a  />  dependence  It  suppressed.  For  accuracy,  ba  * ra  tha  rang*  you 

tpaclfy  la  sufficient  for  tha  far  flald. 

For  both  dBI  and  volts  par  aatar  ,  you  will  ba  proaptad  to  save  tha 
pat  tarn  data.  Whan  you  answer  yaa,  you  arm  proaptad  for  a  fll*  naaa .  Tha 
pattarn  data  will  ba  saved  «a  ASCI  T  laages  tn  thla  flla.  You  aay  uaa  an 
adltor  or  tha  MMIiEC  poat  procaaaor  to  prapara  thla  data  for  plotting  ( *** 
Appandli  B) . 

*1.9  CONFUTE  NEAR  FIELDS 

Tha  N  option  on  tha  NIWIWEC  aanu  la  uaad  to  apaclfy  tha  location  of 
polnta  for  calculation  of  tha  near  alactrlc  and  aagnatlc  field*.  If  the 
currant*  have  not  a  ready  been  computed,  they  will  ba  caaiputad  before  you  can 
proceed.  Ycu  auat  chocs*  alactrlc  or  aagnatlc  field*.  Than  you  will  b a 
proaptad  for  the  flald  location  In  cartesian  coordinates.  You  are  proaptad 

for  tha  Initial  coordinate,  tha  lncraaant  and  tha  mmtoar  of  atapa  for  each  of 

tha  principle  dlr*<tlons.  All  dlaansions  auat  ba  In  aatar*.  Neat,  tha 
radiated  power  levi  la  displayed  and  you  ara  given  a  chance  to  change  thla 
level.  The  near  fields  will  be  scaled  to  tha  power  level  you  specify.  If  you 
specify  zero,  the  original  power  level  will  ba  uaad.  You  will  alao  ba 
proaptad  to  save  the  near  flald  data.  The  field  data  will  ba  stored  In  tha 

disk  file  you  specify.  Use  MW  POST  (sea  Appandli  B)  or  an  adltor  of  your 
choice  to  process  this  data  for  plotting. 

**.10  QUIT 

The  Q  option  on  the  HIWINEC  menu  provides  a  clean  and  efficient  ter¬ 
mination  of  the  MININEC  session. 
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APPENDIX  A 


A  PRE-PROCESSOR  FOR  MININEC 


MNPRE.BAS  converts  NEC  input  dsts  sets  to  HININEC  geometry  specifica¬ 
tions.  First  prepare  an  input  data  file  suitable  for  NEC  (see  reference  5  for 
instructions).  Then  run  HNPRE.BAS.  MNPRE.BAS  will  prompt  you  for  a  disk  file 
name  containing  the  NEC  data  set.  It  will  convert  the  geometry  portion  of  the 
NEC  data  set  into  the  geometry  specifications  required  by  MININEC  and  write 
the  results  into  a  binary  file  named  MININEC. INP.  When  you  run  MININEC,  it 
will  check  MININEC. INP  for  data.  If  empty,  data  entry  is  normal  keyboard 
entry.  If  not  empty,  MININEC  will  read  the  data  and  display  the  segmentation 
data.  The  MININEC  session  will  then  proceed  as  normal  starting  with  the 
quest  ion : 

CHANGE  GEOMETRY  (Y/N)? 

A  source  listing  of  MNPRE.BAS  follows: 


I1 

*  *  •! 
»*.  i1 

i*  v 


10  REM  $TITLE i  ’ IGNOMINI  /£' 

20  OEEINT  I-N 
30  MAXSEGS-75 
40  ON  ERROR  GOTO  320 

50  INPUT  "NEC  Dataaat  Naaa  { . NEC) " ? FS : FS-FS+" . NEC" 

60  OPEN  "MININEC . INP"  AS  12  LEN-30 

70  FIELD  *2,2  AS  S$ ,  4  AS  Xl$,4  AS  Y1S.4  AS  Zl$,4  AS  X2$,4  AS 
Y2$ , 4  AS  Z2$ ,  4  AS  RS 

80  PUT  #2  'Duaaiy  firat  racord  for  no.  of  wiraa 

90  OPEN  FS  FOR  INPUT  AS  «1 

100  ON  ERROR  GOTO  0 

110  IF  EOF ( 1 )  THEN  300 

120  LINE  INPUT  #  1  ,  L$ 

130  IF  LEFTS ( L$ , 2 ) <  > "GW"  THEN  110 
132  PRINT  L$ 

140  I-3-(MID$(LS, 3, 1 )•", " ) i J-INSTR( I ,  L$  ,  "  ,  "  ) 

150  GOSUB  230i IX-XiLSET  S$-MKI$(IX) 

160  NSEGS-NSEGS+IX: IF  NSEGS >MAXSEGS  THEN  280 
170  NWI AES -NW IRES +1 

180  GOSUB  2 30 i LSET  XI $-MKS$ ( X ) 1GO8UB  230>LSET  Y1 $ -MKS$ ( X )> GOSUB 
230iLSET  Z15-MKS$(X) 

190  GOSUB  2 30 1 LSET  X2S-MXSS (X) iOOSUB  2 30 » LSET  Y2S -MKS$ ( X )« GOSUB 
230iLSET  Z2$-MKS$(X) 

200  I-J>1 1 R-VAL(MID$(L$, I ) ) 1  IF  R>0  THEN  LSET  R$-MKS$(R)  ELSE  290 
210  PUT  • 2 
220  GOTO  110 

230  I-J  +  l « J-INSTR( I . LS. " , " ) 1  IF  J  THEN 
X-VAL(MID$(L$. I . J-I ) ) » RETURN 

240  PRINT"Not  anouqh  fialda  in  thia  GW  card." 

250  CLOSE  1 K I LL  "MININEC . INP" 

260  PRINT  1  PR I NT "Any  kay  for  Initial  Optiona" 

270  WHILE  INKEYS-"" 1 WENDi SYSTEM 

280  PRINT“Sagaant  limit  axcaadad  ( "MAXSEGS " ) " : GOTO  250 
290  PRINT"Wira  radiua  nuat  ba  poaitiva  numbar . " t GOTO  250 
300  LSET  S$ -MKI $ ( NWI RES ) 1  PUT  *2,1 
310  CLOSE: SYSTEM 

320  PRINT  F$ "  fila  not  found.  <Ctri-Braak>  to  axit.":PRINT 
330  CLOSE  1 KI LL  "MININEC . INP" 

340  RESUME  50 
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APPENDIX  B 


A  POST-PROCESSOR  MININEC 


HNPOST.BAS  procaaaaa  MININEC  output  data  for  plotting  ualng  tha  CRAPS 
prograai  (aa#  NOSC  TD  820,  "CRAPS:  Graphical  Plotting  Syataa"  by  R.T.  Laird, 
July  1985).  You  aay  atora  MININEC  curranta,  naar  flalda  and  pattarn  data  In  a 
file  of  your  cholca.  Elactlon  to  atora  output  data  In  dlak  flits  la  accoa- 
pllshad  during  a  MININEC  aaaalon.  MNPOST.BAS  will  raad  thaaa  fllaa  and  proaipt 
you  for  tha  data  to  ba  plottad.  MNPOST.BAS  will  racognlsa  tha  typa  of  data 
and  dlaplay  It  for  your  convanianca.  Aftar  you  hava  choaan  tha  data  for  plot¬ 
ting,  tha  alnlaia  and  aailaua  valuaa  ara  coaputad  and  dlaplayad.  You  will  ba 
proaptad  to  adjuat  tha  acala  liaita  or  uaa  thaaa  valuaa.  Than  tha  data  will 
ba  writtan  to  a  flla  you  daalgnata  In  tha  foraiat  raqulrad  by  CRAPS.  A  prograai 
Hating  of  MNPOST.BAS  followa: 


1  REM  ****  MIN I NEC  POST  PROCESSOR  ****  NOSC  CODE  822  JCL  7-86 

2  CLS 

3  COLOR  2,0 

4  DIM 

EX( 300, 4 ) , EY( 300, 4 ) , EZ ( 300, 4 ) , EP ( 300 ) , N( 100 ) , P( 400 ) , X( 400, 3 ) , Z ( 16 
00,2) 

5  PRINT  H  +♦++  MINI  NEC  POST- PROCESSOR  ++*+" 

6  PRINT  : INPUT  "MININEC  OUTPUT  DATA  FILE  (n*»«.OUT)  " ?F$ 

7  T$ -RIGHTS (F$, 4) 

8  IF  LEFTS ( T$ , 1 )  -  M.M  THEN  10 

9  FS-FS+"  .OUT*' 

10  OPEN  FS  FOR  INPUT  AS  *1 

11  INPUT  • 1 , I , PO , G$ 

12  IF  G$-"C"  THEN  GOSUB  19 

13  IF  ( G$ ■ "DM  OR  G$ ■ "V" )  THEN  GOSUB  190 

14  IF  (G$-"E"  OR  G $ ■ “ H " )  THEN  GOSUB  121 

15  PRINT  : INPUT  “CONTINUE  (Y/N)  ";G$ 

16  IF  LEFTS ( G$, 1 )-"N"  THEN  349  ELSE  PRINT  "  “ 

17  GOTO  5 

18  REM  -  PROCESS  CURRENT  DATA  - 

19  J-0 

20  K-0 

21  FOR  L«1  TO  I 

22  J-J  +  l 

23  EX( L. 1 )-J 

24  K-K+l 

25  INPUT  I 1 , EY ( J , 1 ) , EY { J , 2 ) , EY ( J , 3 ) , EY ( J , 4 ) 

26  IF  EY(J,  1)+EY(J,2)+EYU,3)+EY(J,4)»4  THEN  29 

27  J-J+l 

28  GOTO  24 

29  N(L)-K-1 

30  J-J-l 

31  K-0 

32  EX( L, 2 ) -J 

33  NEXT  L 

34  CLOSE  #1 

35  REM  -  DISPLAY  DATA  - 

36  P-0 

37  K*  1 

38  L«1 

39  Nl-1 

40  GOSUB  110 

41  PRINT  : INPUT  "DISPLAY  ANOTHER  WIRE  (Y/N)  "?C$ 

42  IF  LEFTS (C$, 1 )«"N"  THEN  46 

43  L-L+l 

44  IF  L> I  THEN  38 

45  GOTO  40 

46  PRINT 

47  PRINT  "  ++++  SELECT  PLOT  DATA  ++++" 

48  PRINT:  PRINT  "ORDINATE  DATA:" 

49  P-P+l 

50  W-l 

51  IF  1-1  THEN  55 

52  PRINT  : PRINT  "NUMBER  OF  WIRES  -  ";I 


53  PRINT  i INPUT  "USE  DATA  PROM  WHICH  WIRE  M;W 

54  IP  (W<1  OR  W>I)  THEN  52 

55  P(P)-N(W) 

56  PRINT  i  PRINT  “DATA  TYPE «  I  -  REAL" 

57  PRINT  "  2  -  IMAGINARY" 

58  PRINT  "  3  -  MAGNITUDE" 

59  PRINT  "  4  -  PHASE" 

60  INPUT  "ENTER  CHOICE  ";C 

61  IP  (C<0  OR  04)  THEN  56 

62  PRINT  :  PRINT  "ABSCISSA  DATA: 

63  PRINT  :  PRINT  "  1  -  PULSE  POSITION  ON  WIRE" 

64  PRINT  "  2  -  DISTANCE  (■)  ALONG  WIRE" 

65  INPUT  "ENTER  CHOICE  ";J 

66  IP  (J<1  OR  J>2 )  THEN  62 

67  ON  J  GOSUB  80,91 

68  PRINT  iPRINT  "ADD  MORE  DATA  PROM  FILE  "rF$?"  (Y/N)“ 

69  INPUT  C$ 

70  IF  LEFTS (C$. 1 )«"N"  THEN  75 

71  Z(K,1)-1.234 

72  Z(K,2)  — 1.234 

73  K-K+l 

74  GOTO  38 

75  Z(K,  1)— 1.234 

76  Z(K,2)  — 1.234 

77  GOSUB  276 

78  RETURN 

79  REM  - ABSCISSA  TYPE  (J-l) - 

80  N1-EX(W,1) 

81  N2-EX(W.2) 

82  N3-0 

83  FOR  L-Nl  TO  N2 

84  Z(K,1)«N3 

85  Z(K, 2)»EY(L,C) 

86  K-K+l 

87  N3-N3+1 

88  NEXT  L 

89  RETURN 

90  REM  -  ABSCISSA  TYPE  ( J-2)  - 

91  Nl-EX(W.l) 

92  N2-EX(W,2) 

93  PRINT  :  PRINT  "COORDINATES  OF  WIRE  "  ;W 

94  INPUT  "  END1  (X,Y,Z)  ";X1,Y1,Z1 

95  INPUT  "  END2  (X,Y,Z)  " ?X2,Y2,Z2 

96  S- (X2-X1 ) * ( X2-X1 ) + ( Y2-Y1 ) * ( Y2-Y1 ) + ( Z2-Z1 ) * ( Z2-Z1 ) 

97  IF  S>0  THEN  100 

98  PRINT  "WIRE  LENGTH  »  0" 

99  GOTO  93 

100  SaSQR( S ) / (N(W)-l) 

101  X-0 

102  FOR  L«N1  TO  N2 

103  Z(K,1)«X 

104  Z{K,2)»EY<L,C) 

105  X«X+S 

106  K«K+1 
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107  NEXT  L 
106  RETURN 

109  REM  -  DISPLAY  CURRENTS  - 

110  PRINT  s  PRINT  "TOTAL  POWER  RADIATED  -  “ r PO : "  WATTS" 

111  PRINT  i PRINT  "WIRE  NUMBER  “rL:"s" 

112  PRINT  "REAL" , "IMAGINARY" . "MAGNITUDE" . “PHASE" 

113  PRINT  "  NO.  ", "(AMPS)". “(AMPS)". "(DEGREES)" 

114  N 1 “EX ( L , 1 ) 

115  N2-EX(L,2) 

116  FOR  J-Nl  TO  N2 

117  PRINT  USING  "#.###• . r EY ( J , 1 ) . EY ( J . 2 ) . EY ( J . 3 ) r s PRI NT 

USING  "#♦##. #*";EY(J, 4) 

118  NEXT  J 

119  RETURN 

120  REM  -  PROCESS  NEAR  FIELD  DATA  - 

121  FOR  L-l  TO  I 

122  INPUT  #1,  EX(L, 1 ) . EX(L, 2) , EX(L, 3) , EX(L,4) 

123  INPUT  #1,  EY (  L ,  1  )  , EY (  L ,  2  )  , E Y ( L , 3 ) , EY (  L ,  4  ) 

124  INPUT  #1.  EZ ( L , 1) , EZ ( L, 2 ) , EZ ( L, 3) , EZ ( L , 4 ) 

125  INPUT  #1,  EP (  L )  ,  P  (  L ) 

126  INPUT  #1,  X(L,l),X(L,2)fX(L,3) 

127  NEXT  L 

128  CLOSE  #1 

129  REM  -  DISPLAY  DATA  - 

130  L-l 

131  GOSUB  177 

132  PRINT  : INPUT  "DISPLAY  ANOTHER  FIELD  POINT  (Y/N)  "?C$ 

133  IF  LEFTS ( CS. 1)«"N"  THEN  GOTO  137 

134  L-L+l 

135  IF  L> I  THEN  L-l 

136  GOTO  131 

137  LI -1  : L2- I 

138  PRINT  SPRINT  "  ++++  SELECT  PLOT  DATA 

+  +  +  +" 

139  PRINT  : INPUT  "ABSCISSA  DATA  (X/Y/Z)  "?A$ 

140  A$*LEFT$ (AS  # 1 ) 

141  M-0 

142  IF  A$ * "X"  THEN  M-l 

143  IF  A$«"Y"  THEN  M-2 

144  IF  A$="Z"  THEN  M-3 

145  IF  M-0  THEN  GOTO  138 

146  N-0 

147  PRINT  SPRINT  "ORDINATE  DATAs  X  =  X-COMPONENT" 

148  PRINT  "  Y  -  Y-COMPONENT" 

149  PRINT  "  Z  -  Z-COMPONENT" 

150  PRINT  ”  P  =  MAXIMUM  OR  PEAK  VALUE" 

151  INPUT  "ENTER  CHOICE  OF  FIELD  COMPONENT  (X/Y/Z/P)  "  ;  D$ 

152  D$=LEFTS (D$ , 1 ) 

153  IF  D$="P"  THEN  GOTO  164 

154  PRINT  "DATA  TYPEs" 

155  PRINT  "  1  -  REAL" 

156  PRINT  "  2  -  IMAGINARY" 

157  PRINT  "  3  -  MAGNITUDE" 

158  PRINT  "  4  -  PHASE" 


159 

INPUT  CHOICE: " 

•  J 

160 

IF 

D$-“X"  THEN 

GOTO 

16S 

161 

IF 

D$-"Y"  THEN 

GOTO 

166 

162 

IF 

DS-"Z  THEN 

GOTO 

16" 

N-N*  1 

: Zl L,  1) 

«  X  (  N  ,  M  ' 

:  Z  i 

1. .  2  ' 

=  EP 

i  N  ' 

:  NEXT  1 

N-N*  1 

:  Z  i  L  ,  1  > 

-  X  (  N  ,  M  1 

.  z  , 

L  .  2  1 

i  -EX 

N  .  0 

:  NEXT 

N-N*  1 

:  Z  (  L  .  11 

-Xt N . M 1 

:  Z  ! 

1. ,  2 

'  'EY 

iN,.’ 

1  : NEXT 

N  -N>  1 

: Z( L.  1' 

-X ( N . M ) 

:  Z  ( 

L .  2 

1  -  EZ 

(  N  .  0 

I  : NEXT 

163  GOTO  147 

164  FOR  L-Ll 
: GOTO  168 

1 6 5  FOR  L-Ll 
: GOTO  168 

166  FOR  L-Ll 
: GOTO  168 


168  L2 -L2 ♦ 1 

169  PRINT  : PRINT  "ADD  MORE  DATA  FROM  FILE  • F$ •  (Y  N '  • 

170  INPUT  CS 

171  if  LEFTS  (CS  .  1  >  -"N"  THEN  GOTO  173 

172  Z( L2. 1 ) -l . 234  : Z ( L2 , 2 ) - - l . 2 34  :L1-L2-M  :L2*L2*I  : GOTO  138 

173  Z(L2,  1  )--l  .234  : Z ( L2 . 2 ) -- 1 . 2 34 

174  GOSUB  276 

175  RETURN 

176  REM  -  DISPLAY  NEAR  FIELD  DATA  - 

177  PRINT  : PRINT  "  ***  NEAR  FIELD  DATA  FROM 

FILE  rFS:'"  ***" 

170  PRINT  SPRINT  ’■  FIELD  POINT:  X  -  ":X(L,1):",  Y 

";X'L.2) : " ,  Z  «  ":X(L,3) 

179  PRINT  "VECTOR" , "REAL" . "IMAGINARY" , "MAGNITUDE" , "PHASE- 

180  IF  G$-"E"  THEN  PRINT 

"COMPONENT" , " (V/M) " , " (V/M) “ . " (V/M) " . " ( DEG ) " 

181  IF  G$«"H"  THEN  PRINT 

"COMPONENT" , " (AMPS/M) " , " (AMPS/M) " , " (AMPS/M) " , ” ( DEG ) " 


182  PRINT  "  X  sPRINT  USING  "#.#### 

" ;EX(L, 1 ) ,EX(L, 2) , EX(L, 3) ; sPRINT  USING  "####.#*" • EX ( L , 4 ) 

183  PRINT  "  Y  sPRINT  USING  "#.#### 

" ;EY(L, 1 ) ,EY(L. 2) ,EY(L, 3)  ; sPRINT  USING  "####.##" r  EY ( L , 4 ) 

184  PRINT  "  Z  " ; sPRINT  USING  "#.#### 

" ;EZ(L, 1 ) ,EZ(L, 2) ,EZ(L, 3) ; sPRINT  USING  "####.##"; EZ ( L , 4 ) 

185  IF  G$*"E"  THEN  PRINT  "MAXIMUM  OR  PEAK  FIELD  =  " r  EP ( L ) ; "  V/M" 

186  IF  G$  =  "H"  THEN  PRINT  "MAXIMUM  OR  PEAK  FIELD  =  " r  EP ( L ) ; " 

AMPS/M" 

187  PRINT  "RADIATED  POWER  =  ";P(L)r"  WATTS" 

188  RETURN 

189  REM  -  PROCESS  PATTERN  DATA  - 

190  IF  G$="V"  THEN  GOTO  197 

191  REM  INPUT  DATA  IN  DB 

192  FOR  L=1  TO  I 

193  INPUT  #1,  X(L,1),X(L,2), EX { L , 1 ) , EX ( L , 2 ) , EX ( L , 3 ) 

194  NEXT  L 

195  REM  INPUT  DATA  IN  V/M 

196  GOTO  201 

197  INPUT  #1,  RO 

198  FOR  L=1  TO  I 

199  INPUT  #1,  X(L,1),X{L,2), EX ( L , 1 ) , EX ( L , 2 ) , EY ( L , 1 ) , EY ( L , 2 ) 

200  NEXT  L 

201  CLOSE  #1 


4T  4*r‘  4 


202  REM  - DISPLAY  DATA - 

20  3  J- I  :  K  »  1  0  -IF  K  > I  THEN  K-I 

204  IF  G$  * " D"  THEN  GOSUB  232 

205  IF  G$  * " V "  THEN  GOSUB  238 

206  PRINT  : INPUT  "DISPLAY  MORE  PATTERN  DATA  (Y  N)  ";C$ 

207  IF  LEFT  S  (  C  $  ,  1  )  *  "  N ''  THEN  GOTO  212 

208  IF  K-I  THEN  GOTO  203 

209  J-K  : K -K ♦  1  0  :IF  J>I  THEN  J-K-9 
2  10  IF  K  >  I  THEN  K-I 

211  GOTO  204 

212  LI  -  1  : L2 - I 

213  PRINT  : PRINT  "  +♦♦  +  SELECT  PLOT  DATA 


♦  *  *  ♦ 

214  PRINT  : INPUT  "ABSCISSA  DATA:  (Theta/Phi)  ":A$ 

215  A$ -LEFTS ( A$ , 1 ) 

216  M-0 

217  IF  A$-"T"  THEN  M- 1 

218  IF  AS* " P"  THEN  M-2 

219  N-0 

220  IF  M-0  THEN  GOTO  213 

221  IF  G$="D"  THEN  GOSUB  252 

222  IF  G$  = " V"  THEN  GOSUB  262 

223  L2-L2+1 

224  PRINT  : PRINT  "ADD  MORE  DATA  FROM  FILE  ' " ; FS ; (Y/N)  " ; 

225  INPUT  C$ 

226  IF  LEFTS (C$,1)="N“  THEN  GOTO  228 

227  Z ( L2 , 1 ) -1 -234  : Z ( L2 , 2 ) — 1 . 234  :L1*L2+1  :L2=L2+I  :GOTO  213 

228  Z(L2, 1  )  — 1 .234  :  Z  (  L2 , 2  )  — 1 . 2  34 

229  GOSUB  276 

230  RETURN 

231  REM  -  DISPLAY  PATTERN  DATA  (DBI)  - 

232  PRINT  : PRINT  "  RADIATION  PATTERN  DATA" 

233  PRINT  .-PRINT 

"ZENITH" , "AZIMUTH" , "VERTICAL" , "HORIZONTAL" , "TOTAL" 

234  PRINT  " (THETA) (PHI )", "PATTERN  ( dB) PATTERN 
(dB) ", "PATTERN  (dB)" 

235  FOR  L=J  TO  K  :PRINT  X ( L , 1 ) , X ( L , 2 ) , EX ( L , 1 ) , EX ( L , 2 ) , EX ( L , 3 ) 

: NEXT  L 

236  RETURN 

237  REM  -  DISPLAY  PATTERN  DATA  (V/M)  - 

238  PRINT  : PRINT  "  RADIATION  PATTERN  DATA" 

239  PRINT  :PRINT  "RADIAL  DISTANCE  =  " r RO  ;  "  METERS" 

240  PRINT  "POWER  LEVEL  =  " ; PO ; "  WATTS" 

241  PRINT  : PRINT"  ZENITH  AZIMUTH","  E( THETA) 

E ( PHI ) " 

242  PRINT  "(THETA)  { PHI ) " , "MAG ( V/M)  PHASE( DEG) " , "MAG( V/M) 

PHASE (DEG) " 

243  FOR  L=J  TO  K 


244  PRINT  USING 

245  PRINT  USING 

246  PRINT  USING 

247  PRINT  USING 

248  PRINT  USING 

249  NEXT  L 


###.##  ”;X(L,1),X(L,2); 

##.### - "  r  EX  ( L ,  1 )  ? 

####.##  " ; EX( L, 2 ) ; 

##.*## - -  ;  EY  ( L ,  1 )  r 

####.##"; EY ( L , 2 ) 
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250  RETURN 

251  REM  -  SAVE  PATTERN  DATA  ( BDI )  - 

252  PRINT  : INPUT  "ORDINATE  DATA  ( Vertical/Hor i zontal /Tota 1 )  " 

253  D$=LEFT$ ( D$ , 1 ) 

254  J=0 

255  IF  D$="V"  THEN  J=1 

256  IF  D$="H"  THEN  J=2 

257  IF  D$="T"  THEN  J=3 

258  IF  J=0  THEN  GOTO  260 

259  FOR  L=L1  TO  L2  :N=N+1  : Z ( L, 1 ) =X ( N, M )  : Z ( L , 2 ) =EX ( N , J )  : NEX 

260  RETURN 

261  REM  -  SAVE  PATTERN  DATA  (V/M)  - 

262  PRINT  : PRINT  "ORDINATE  DATA:  T  =  E-THETA  COMPONENT' 

263  PRINT  "  P  =  E-PHI  COMPONENT" 

264  INPUT  "ENTER  CHOICE  OF  COMPONENT  (T/P)  " ; D$ 


265  D$=LEFT$ ( D$ , 1 ) 

266  PRINT  : PRINT  "DATA  TYPE: 


1  -  MAGINTUDE" 


267  PRINT  " 

268  INPUT  "CHOICE:" ;J 

269  IF  D$="T"  GOTO  272 

270  IF  D$="P"  GOTO  273 

271  GOTO  274 

272  FOR  L=L1  TO  L2  :N=N+1 
: GOTO  274 

273  FOR  L=L1  TO  L2  :N=N+1 


2  -  PHASE" 

:Z(L,1)=X(N,M) 
:  Z  (  L ,  1 )  =X  ( N ,  M ) 


Z  (  L  ,  2  )  =EX  (  N  ,  0 
Z  (  L ,  2  )  =  E  Y  ( N ,  C 


n  f  v 


:  NF' 


274  RETURN 

275  REM  - FIND  MAX  &  MIN  - 

276  PRINT  : PRINT  "PLOT  FORMAT:  1  -  ONE  ORDINATE  SC  ALE 

277  PRINT  "  2  -  TWO  ORDINATE  SCALES " 

278  PRINT  "  3  -  POLAR  OR  SMITH  CHART 

279  INPUT  "ENTER  CHOICE  ";C  :PRINT 

280  IF  (C<1  OR  C> 3 )  THEN  GOTO  276 

281  IF  C=2  THEN  PRINT  " -  TWO  CURVES  IS  ASSUMED  - 

282  Ll=l  : L2=I 

283  IF  G$="C"  THEN  300 

284  XL=9 . 999999E+3 5  :XH=-9 . 999999E+35 

285  YL=9 . 999999E+35  :YH=-9 . 999999E+35 

286  FOR  L=L1  TO  L2 

287  IF  Z ( L, 1 ) <XL  THEN  XL=Z(L,1) 

288  IF  Z ( L , 1 ) >XH  THEN  XH-Z(L.l) 

289  IF  Z ( L, 2 ) <YL  THEN  YL=Z(L,2) 

290  IF  Z ( L , 2 ) > YH  THEN  YH - Z ( L , 2 ) 

291  NEXT  L 

292  L2-L2+1 

293  IF  Z(L2,1)— 1.2  34  THEN  GOTO  32  1 

294  IF  C><2  THEN  GOTO  298 

295  Y1 *YL  : Y2-YH 

296  L1-L2+1  : L2-L2+I 

297  GOTO  285 

298  LI  *L2  + 1  :L?*L2-*-I  :GOTO  Jf 

299  REM  ---  CURRENTS  MAX  k  min 

300  IF  C<  >  2  THEN  304 

301  IF  P-2  THEN  304 

302  PRINT  :  PRINT  BUT.  THF. RF  k R F 


308  FOR  L-Ll  TO  L 2 

309  IT  Z( L, 1 ) <XL  THBM  XL-I(L.l) 

310  IP  X(  L,  1 ) >XH  THBM  XH-I(L.l) 

311  IP  X(L,2)<YL  THBM  YL-KL.2) 

312  IP  Z( L,  2 ) >YH  THBM  YH-KL.2) 

313  MBXT  L 

314  IP  J-2  THBM  318 

315  IP  C< >2  THBM  318 

316  Yl-YLtY2-YH 

317  YL*9 . 999999E+35  «YH— 9. 999999E+35 

318  LI  «L2+2 

319  NEXT  J 

320  L2-K 

321  PRINT  >PRINT  M  ", "MINIMUM", "MAXIMUM" 

322  PRINT  “ABSCISSA" , XL. XH 

323  IP  C-2  THEN  PRINT  "ORDINATE" , Y1 , Y^ , "  ( LBPT  SIDE)" 

324  PRINT  "ORDINATE", YL.YH 

325  IP  C-3  THEN  GOTO  337 

326  PRINT  tINPUT  "CHANGE  ABSCISSA  RANGE  (Y/N)  “ ;C$ 

327  IP  LEFTS ( CS , 1 ) “"N"  THEN  GOTO  330 

328  INPUT  "NEW  LOWER  LIMIT  -  ";XL 

329  INPUT  "NEW  UPPER  LIMIT  -  ";XH 

330  PRINT  {INPUT  "CHANGE  ORDINATE  RANGE  (Y/S)  "  *C$ 

331  IP  LEFTS (C$. 1)-"N"  THEN  GOTO  337 

332  IF  C>  <2  THEN  GOTO  335 

333  INPUT  "NEW  LOWER  LIMIT  (LEFT  SIDE)  -  ";Y1 

334  INPUT  “NEW  UPPER  LIMIT  (LEFT  SIDE)  ■  "?Y2 

335  INPUT  "NEW  LOWER  LIMIT  -  " ? YL 

336  INPUT  "NEW  UPPER  LIMIT  *  " ; YH 

337  PRINT  {INPUT  "PLOT  DATA  FILE  ( nai« . DAT )  ":FS 

338  T$»RIGHT$(F$,4) 

339  IF  LEFTS (T$,l)  •  THEN  341 

340  F$ssF$+"  . DAT" 

341  OPEN  F$  FOR  OUTPUT  AS  #1 

342  IF  C»1  THEN  PRINT  *1,  XL r " , " ; XH ; " , " ; YL ? " , " ; YH 

343  IF  C»2  THEN  PRINT  #1,  XL; " , " ;XH; " , " rYl ; " , " ? Y2 r " , " ; YLr " , " r YH 

344  FOR  L*1  TO  L2 

345  PRINT  #1,  Z(L,1);",";Z(L,2) 

346  NEXT  L 

347  CLOSE  #1 

348  RETURN 

349  SYSTEM 
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APPINDIX  c 


MIN  IN  EC  PM6MM  LIST  INC 


NININCC  Compilation 

Futtit  run  tlaes  for  MININIC3  have  boon  achieved  using  87BASIC/ 
INLINE (TN) ,  MlcroNay's  BASIC  compiler  post  procossor  uhieh  gonorstos  ln-llno 
0067  code  for  sll  floating  point  otprossions.  Tho  following  toblo  glvos  som 
idea  of  tho  dlfforoneo  In  run  tlaes  for  tho  matrix  fill  In  tho  sample  problaa 
in  NOSC  TO  516,  Appondli  B. 


BASICA  In tor pro tor 
IB*  BASIC  Compiler 
87BASIC  Coapllor 
87BASIC/ INLINE 


A  1/2  hours 
22  alnutos 
8  alnutos 
A  alnutos 


87BASIC/ INLINE  is  available  froa  MicroWay,  P.0.  Box  79.  Kingston,  Moss. 
0236A  Phono  (617)  7A6-73A1.  Tho  eurront  prlco  of  tho  pockago  Is  two  hundrod 
dollars. 

MINIMEC3.BAS  asy  bo  run  with  tho  BASICA  Intorprotor,  but  tho  asxlaua 
nusber  of  pulsos  aust  bo  roducod  to  A2  with  10  wlros.  A  asxlaua  of  50  wires 
and  50  pulses  asy  bo  used  with  tho  IBM  BASIC  Coapllor.  Tho  other  compilers 
allow  70  pulsos. 

A  program  listing  dimensioned  for  tho  IBM  BASIC  Coapllor  follows: 
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MIHXBfBC ( 3 ) 


NOSC  CODE  822  (JCL  CHANGE  6) 


9-25-86 

2  DEFINT  I,J,K,N 

3  DIH  K1 (6, 2) ,Q(14) 

4  REM  -  MAXIMUM  NUMBER  OF  SEGMENTS  (PULSES  ♦  2  *  WIRES)  -  150 

5  M8-150 

6  DIM  X(150),Y(150),Z(150) 

7  REM  -  MAXIMUM  NUMBER  OF  WIRES  -  50 

8  MW- 50 

9  DIM  A( 50) , CA( 50) , CB( 50) ,CG( 50) , J1 ( 50) , J2( 50, 2 ) ,N( 50, 2 ) , S( 50) 

10  REM  -  MAXIMUM  NUMBER  OF  LOADS  -  11 

11  ML- 11 

12  REM  -  MAXIMUM  ORDER  OF  S-PARAMETER  LOADS  -  8 

13  MA-8 

14  DIM  LA ( 2 , 1 1 , 8 ) , LP ( 1 1 ) , LS ( 1 1 ) 

15  REM  -  MAXIMUM  NUMBER  OF  MEDIA  «  6 

16  MM-6 

17  REM  -  H  MUST  BE  DIMENSIONED  AT  LEAST  6 

18  DIM  H(6) , T ( 6 ) ,U(6) , V( 6 ) , Z1 ( 6 ) , Z2 ( 6 ) 

19  REM  -  MAXIMUM  NUMBER  OF  PULSES  -  50 

20  MP-50 

21  DIM  C%(50, 2) , Cl ( 50 ) , CR( 50 ) ,P( 50) ,W%(50) 

22  DIM  ZR( 50, 50) , ZI ( 50, 50) 

23  REM  -  ARRAYS  E.L  &  M  DIMENSIONED  TO  MW+MP-100 

24  DIM  E( 100) , L( 100) ,M( 100) 

25  COLOR  2,0 

26  GOTO  1499 

27  REM  **********  KERNEL  EVALUATION  OF  INTEGRALS  12  &  13 

28  IF  K<0  THEN  33 

29  X3-X2+T* ( V1-X2 ) 

30  Y3-Y2+T* ( V2-Y2 ) 

31  Z3=Z2+T* ( V3-Z2 ) 

32  GOTO  36 

33  X3=V1+T* (X2-V1 ) 

34  Y3-V2+T* ( Y2-V2 ) 

35  Z3=V3+T* ( Z2-V3) 

36  D3=X3*X3+Y3*Y3+Z3*Z3 

37  REM  -  MOD  FOR  SMALL  RADIUS  TO  WAVELENGTH  RATIO 

38  IF  A( P4 ) <  =SRM  THEN  D=SQR( D3 ) : GOTO  49 

39  D-D3+A2 

40  IF  D>0  THEN  D=SQR( D) 

41  REM  -  CRITERIA  FOR  USING  REDUCED  KERNEL 

42  IF  161=0  THEN  49 

43  REM  -  EXACT  KERNEL  CALCULATION  WITH  ELLIPTIC  INTEGRAL 

44  B=D3/ ( D3+4*A2 ) 

45  W0=C0+B* (Cl+B* (C2+B* (C3+B*C4 ) ) ) 

46  W1 =C5+B* (C6+B* ( C7+B* ( C8+B*C9 ) ) ) 

47  VO=(WO-W1*LOG( B) ) *SQR( 1-B) 

48  T3=T3+ ( VO+LOG ( D3/ ( 64*A2 ) )/2)/P/A(P4)-l/D 

49  B1=D*W 

50  REM - EXP(-J*K*R)/R 

51  T3=T3+C0S(B1)/D 
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-4V*  * 


52  T4*T4-SIH(B1)/D 

53  RETURN 


54  RB4  *****  PSI ( PI , P2 , P3 )  -  T1  +  J  *  T2  ********** 

55  REM  -  ENTRIES  REQUIRED  FOR  NEAR  FIELD  CALCULATION 


56  Xl-XO+Pl*T5/2 

57  Yl-Y0+Pl*T6/2 

58  Zl-Z0+Pl*T7/2 

59  X2*X1-X(P2) 

60  Y2-Y1-Y(P2) 

61  Z2-Z1-K*Z(P2) 

62  V1-X1-X(P3) 

63  V2«Y1-Y(P3) 

64  V3»Z1-K*Z(P3) 

65  GOTO  135 

66  I4-INT(P2) 

67  15*14+1 

68  X2*X0-(X(I4)+X(I5) )/2 

69  Y2*Y0- ( Y( 14 ) +Y ( 1 5 ) ) / 2 

70  Z2*Z0-K* ( Z ( 14 ) +Z ( 1 5 ) )/2 

71  V1*X0-X(P3) 

72  V2*Y0-Y(P3) 

73  V3*Z0-K*Z(P3) 

74  GOTO  135 

75  X2*X0-X(P2 ) 

76  Y2*Y0-Y(P2) 

77  Z2-Z0-K*Z(P2) 

78  I4-INT(P3) 

79  15*14+1 

80  V1*X0-(X( I4)+X( 15) )/2 

81  V2*Y0-(Y( I4)+Y( 15) )/2 

82  V3*Z0-K* (Z( I4)+Z( 15) )/2 

83  GOTO  135 


84  REM  -  ENTRIES  REQUIRED  FOR  IMPEDANCE  MATRIX  CALCULATION 

85  REM  -  S (M)  GOES  IN  (X1,Y1,Z1)  FOR  SCALAR  POTENTIAL 

86  REM  -  MOD  FOR  SMALL  RADIUS  TO  WAVE  LENGTH  RATIO 


87  FVS*1 

88  IF  K<1  THEN  94 

89  IF  A( P4 ) >SRM  THEN  94 

90  IF  ( P3*P2+1  AND  PI* ( P2+P3 ) /2 )  THEN  91  ELSE  94 

91  Tl=2*LOG(S(P4)/A(P4) ) 

92  T2*-W*S(P4) 

93  RETURN 

94  I4*INT(P1 ) 

95  15*14+1 

96  XI* (X( I4)+X( 15) )/2 

97  Y1*(Y( I4)+Y( 15) )/2 

98  Z1*(Z( I4)+Z( 15) )/2 

99  GOTO  113 


100  REM  -  S(M)  GOES  IN  (X1,Y1,Z1)  FOR  VECTOR  POTENTIAL 

101  REM  -  MOD  FOR  SMALL  RADIUS  TO  WAVE  LENGTH  RATIO 


102  FVS-0 

103  IF  K< 1  THEN  109 

104  IF  A( P4 ) > -SRM  THEN  109 
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105  IF  (I-J  AND  P3-P2+.5)  THEN  106  ELSE  109 

106  Tl-LOG(S(P4)/A(P4)) 

107  T2— W*S(P4)/2 

108  RETURN 

109  X1-X<P1) 

110  Yl-Y(Pl) 

111  Zl-Z(Pl) 

112  REM  -  S(U)-S(M)  GOES  IN  (X2.Y2.Z2) 

113  I4-INT(P2) 

114  IF  I4-P2  THEN  120 

115  I5-I4+1 

116  X2«(X(I4)+X(I5) )/2-Xl 

117  Y2=(Y(I4)+Y(l5))/2-Yl 

118  Z2-K*(Z(I4)+Z(I5) )/2-Zl 

119  GOTO  124 

120  X2=X(P2)-X1 

121  Y2«Y(P2)-Y1 

122  Z2=K*Z(P2)-Z1 

123  REM  -  S( V) -S(M)  GOES  IN  (V1.V2.V3) 

124  I4«INT(P3) 

125  IF  I4-P3  THEN  131 

126  1 5*14+1 

127  V1«(X(I4)+X(I5) )/2-Xl 

128  V2* ( Y( 14) +Y( 15 ) ) / 2-Y1 

129  V3-K*(Z(I4)+Z(I5) )/2-Zl 

130  GOTO  135 

131  V1»X(P3)-X1 

132  V2*Y(P3)-Y1 

133  V3=K*Z(P3)-Z1 

134  REM  -  MAGNITUDE  OF  S(U)  -  S(M) 

135  DO=X2*X2+Y2*Y2+Z2*Z2 

136  REM  -  MAGNITUDE  OF  S(V)  -  S(M) 

137  IF  D0>0  THEN  DO-SQR(DO) 

138  D3»V1*V1+V2*V2+V3*V3 

139  IF  D3>0  THEN  D3*SQR(D3) 

140  REM  -  SQUARE  OF  WIRE  RADIUS 

141  A2=A(P4)*A(P4) 

142  REM  -  MAGNITUDE  OF  S(V)  -  S(U) 

143  S4= ( P3-P2 ) *S( P4) 

144  REM  -  ORDER  OF  INTEGRATION 

145  REM  -  LTH  ORDER  GAUSSIAN  QUADRATURE 

146  T1=0 

147  T2*0 

148  161=0 

149  F2=l 

150  L=7 

151  T= ( D0+D3 ) /S ( P4 ) 

152  REM  -  CRITERIA  FOR  EXACT  KERNEL 

153  IF  T> 1 . 1  THEN  165 

154  IF  C$  =  "N"  THEN  165 

155  IF  J 2 ( W% ( I ) , 1 ) = J 2 ( W% ( J ) , 1 )  THEN  160 

156  IF  J2(W%( I ) , 1 )=J2(W%(J) , 2)  THEN  160 

157  IF  J2(W%( I) , 2)=J2(W%( J) , 1)  THEN  160 


158  IF  J2(W% ( I ) , 2 ) -J2 ( W% ( J) #  2 )  THEN  160 

159  GOTO  165 

160  IF  A( P4 ) >SRM  THEN  162 

161  IF  FVS-1  THEN  91  ELSE  106 

162  F2-2* (P3-P2 ) 

163  161  -  ( l-LOG( S4/F2/8/M  P4 ) ) ) /P/A( P4) 

164  GOTO  167 

165  IF  T>6  THEN  L-3 

166  IF  T>10  THEN  L«1 

167  I5-L+L 

168  T3-0 

169  T4-0 

170  T«(Q(L)+.5)/F2 

171  GOSUB  28 

172  T«( .5-Q(L) )/F2 

173  GOSUB  28 

174  L-L+l 

175  T1«T1+Q(L)*T3 

176  T2«T2+Q(L)*T4 

177  L-L+l 

178  IF  L<15  THEN  168 

179  T1-S4* (T1+I6 l ) 

180  T2-S4*T2 

181  RETURN 

182  REM  **********  COMPLEX  SQUARE  ROOT  ********** 

183  REM  -  W6+I*W7-SQR(Z6+I*Z7) 

184  T6-SQR( (ABS(Z6)+SQR(Z6*Z6+Z7*Z7) )/2) 

185  T7«ABS(Z7)/2/T6 

186  IF  Z6<0  THEN  191 

187  W6-T6 

188  W7-T7 

189  IF  Z7<0  THEN  W7— T7 

190  RETURN 

191  W6-T7 

192  W7-T6 

193  IF  Z7<0  THEN  W7—T6 

194  RETURN 

195  REM  **********  IMPEDANCE  MATRIX  CALCULATION  ********** 

196  IF  FLG-1  THEN  428 

197  IF  FLG=2  THEN  477 

198  REM  -  BEGIN  MATRIX  FILL  TIME  CALCULATION 

199  OT$-TIME$ 

200  Q$-"MATRIX  FILL  " 

201  PRINT 

202  PRINT  "BEGIN  ";Q$ 

203  REM  -  ZERO  IMPEDANCE  MATRIX 

204  FOR  1=1  TO  N 

205  FOR  J=1  TO  N 

206  ZR( I , J) -0 

207  ZI(I,J)=0 

208  NEXT  J 

209  NEXT  I 

210  REM  -  COMPUTE  ROW  I  OF  MATRIX  (OBSERVATION  LOOP) 
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211  FOR  1-1  TO  N 

212  I 1 -ABS ( C% (1,1)) 

213  1 2 -ABS ( C% ( 1 , 2 )  ) 

214  F4-SGN (Cl (1,1)) *8(11) 

215  F5-SGN(C%(I,2))*S(I2) 

216  REM - R(M  ♦  1/2)  -  R(M  -  1/2)  HAS  COMPONENTS  (T5,T6,T7) 

217  T5«F4*CA(I1)+F5*CA(I2) 

218  T6-F4*CB(  II )  +F5*CB(  12 ) 

219  T7-F4*CG( II )+F5*CG( 12) 

220  IF  C%(I,1)—  C%(I,2)  THEN  T7 -S ( I 1 ) * ( CG ( I 1 ) -fCG (12)) 

221  REM  -  COMPUTE  COLUMN  J  OF  ROW  I  (SOURCE  LOOP) 

222  FOR  J-l  TO  N 

223  J1-ABS(CI(J,1)) 

224  J2-ABS ( C% ( J , 2 ) ) 

225  F4-SGN ( C% ( J ,  1 )  ) 

226  F5*SGN(C%( J, 2) ) 

227  F6-1 

228  F7-1 

229  REM - IMAGE  LOOP 

230  FOR  K-l  TO  G  STEP  -2 

231  IF  C% ( J , 1 ) < > -C% ( J ,  2 )  THEN  235 

232  IF  K<0  THEN  332 

233  F6-F4 

234  F7-F5 

235  F8-0 

236  IF  K<0  THEN  248 

237  REM  -  SET  FLAG  TO  AVOID  REDUNANT  CALCULATIONS 

238  IF  II <>I2  THEN  246 

239  IF  ( CA (II) +CB (II)) -0  THEN  241 

240  IF  C«(X,1)<>C«(I,2)  THEN  246 

241  IF  J1<>J2  THEN  246 

242  IF  (CA( J1)+CB(J1) )-0  THEN  244 

243  IF  C%(J, 1) <>C%(J,2)  THEN  246 

244  IF  Il-Jl  THEN  F8-1 

245  IF  I-J  THEN  F8-2 

246  IF  ZR( I , J ) < > 0  THEN  317 

247  REM  -  COMPUTE  PSI (M,N,N+l/2 ) 

248  P1=2*W% ( I ) +1-1 

249  P2*2*W% ( J) +J-1 

250  P3-P2+.5 

251  P4-J2 

252  GOSUB  102 

253  U1-F5*T1 

254  U2-F5*T2 

255  REM  -  COMPUTE  PSI ( M , N-l /2 , N) 

256  P3-P2 

257  P2=P2- . 5 

258  P4=J1 

259  IF  F8<2  THEN  GOSUB  102 

260  VI =F4*T1 

261  V2=F4*T2 

262  REM  -  S ( N+l/ 2 ) *PS I ( M, N, N+l/2 )  +  S (N-l/2 ) *PSI (M,N-l/2 , N) 

263  X3=U1 *CA ( J2 ) +V1 *CA{ J 1 ) 
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264  Y3«U1*C»(J2)*V1*CB(J1) 

265  Z3-(F7*U1*CG(J2)+F6*V1*CG(J1) )*K 

266  MEM  -  MEAL  PART  OF  VECTOR  POTENTIAL  CONTRIBUTION 

267  D1-W2*(X3*T5+Y3*T6+Z3*T7) 

268  X3-U2*CA(J2)+V2*CA(J1) 

269  Y3-U2*CB(J2)+V2*CB<J1) 

270  Z3-(F7*U2*CG(J2)+F6*V2*CQ(J1 ) )*K 

271  REM  -  IMAGINARY  PART  OF  VECTOR  POTENTIAL  CONTRIBUTION 

272  D2«W2*(X3*T5+Y3*T6+Z3*T7) 

273  REM - COMPUTE  PS  I  (M+l  /  2 ,  N,  N-f  1 ) 

274  P1-P1+.5 

275  IF  F8-2  THEN  Pl-Pl-1 

276  P2-P3 

277  P3-P3+1 

278  P4-J2 

279  IF  F8< > 1  THEN  283 

280  U5«F5*U1+T1 

281  U6-F5*U2+T2 

282  GOTO  291 

283  GOSUB  87 

284  IF  F8<2  THEN  288 

285  U1»(2*T1-4*U1*F5)/S(J1) 

286  U2-(2*T2-4*U2*F5)/S(J1) 

287  GOTO  314 

288  U5-T1 

289  U6-T2 

290  REM  -  COMPUTE  PSI (M-l/2 , N, N+l ) 

291  Pl-Pl-1 

292  GOSUB  87 

293  U1«(T1-U5)/S(J2) 

294  U2-(T2-U6)/S(J2) 

295  REM  -  COMPUTE  PSI (M+l/2 , N-l , N) 

296  Pl-Pl+1 

297  P3-P2 

298  P2*P2-1 

299  P4-J1 

300  GOSUB  87 

301  U3-T1 

302  U4-T2 

303  REM  -  COMPUTE  PSI (M-l/2 , N-l , N) 

304  IF  F8<1  THEN  308 

305  T1-U5 

306  T2-U6 

307  GOTO  311 

308  Pl-Pl-1 

309  GOSUB  87 

310  REM  -  GRADIENT  OF  SCALAR  POTENTIAL  CONTRIBUTION 

311  U1»U1+(U3-T1)/S(J1) 

312  U2-U2+(U4-T2)/S(J1) 

313  REM  -  SUM  INTO  IMPEDANCE  MATRIX 

314  ZR( I,J)-ZR(I,J) +K* ( Dl+Ul ) 

315  ZI(I,J)-ZI(I,J)+K*(D2+U2) 

316  REM  -  AVOID  RE DUN ANT  CALCULATIONS 


317  IF  J<I  THEM  332 
314  IE  F8-0  THEM  332 

319  ZR(J, I)-ZR( I,J) 

320  ZI(J, I)-ZI(I.J) 

321  REM  -  SEGMENTS  ON  SAME  WIRE  8AME  DISTANCE  APART  HAVE  SAME 

Z 

322  Pl-J+1 

323  IP  PI >N  THEN  332 

324  IP  C»(P1, 1) <>C%(P1,2)  THEN  332 

325  IF  C%(P1,2)-C%(J,2)  THEN  328 

326  IP  C% (PI , 2 ) <  >-C% ( J,  2 )  THEN  332 

327  IP  (CA( J2)+CB( J2) ) <  >0  THEN  332 

328  P2-I+1 

329  IP  P2>N  THEN  332 

330  ZR( P2 ,  PI ) **ZR(  I ,  J) 

331  ZI(P2,P1)«ZI(I,J) 

332  NEXT  K 

333  NEXT  J 

334  PCT-I/N 

335  GOSUB  1601 

336  NEXT  I 

337  REM  -  END  MATRIX  PILL  TIME  CALCULATION 

338  T$-TIME$ 

339  GOSUB  1591 

340  PRINT  #3,"  " 

341  PRINT  #  3 , "PI LL  MATRIX  :  MrT$ 

342  REM  **********  ADDITION  OP  LOADS  ********** 

343  IP  NL-0  THEN  377 

344  F5-2*P*F 

345  FOR  1-1  TO  NL 

346  IF  L$-"N"  THEN  366 

347  REM  -  S-PARAMETER  LOADS 

348  Ul-0 

349  U2*0 

350  Dl-0 

351  D2=0 

352  S*1 

353  FOR  J«0  TO  LS(I)  STEP  2 

354  U1»U1+LA( 1, I,J)*S*F5“J 

355  D1*D1+LA(2, I, J)*S*F5AJ 

356  L-J+l 

357  U2»U2+LA{1, I,L)*S*F5“L 

358  D2*D2+LA( 2,I,L)*S*F5“L 

359  S— S 

360  NEXT  J 

361  J  *LP ( I ) 

362  D»D1*D1+D2*D2 

363  LI-(U2*D1-D2*U1)/D 

364  LR»(U1*D1+U2*D2)/D 

365  GOTO  369 

366  LR«LA( 1,1,1) 

367  LI«LA(2,I,1) 

368  J*LP( I ) 
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369  F 2-1/M 

370  IF  C% ( J , 1 ) <  > -C  (  J  ,  2  )  THEM  372 

371  IF  K <0  THEM  F 2-2/M 

372  ZR( J, J)-ZR( J. J)+F2*LI 

373  ZI( J, J)-ZI ( J, J)-F2*LR 

374  NEXT  I 


375  REM  **•*••***•  IMPEDANCE  MATRIX  FACTORIZATION 

376  REM  -  BEGIN  MATRIX  FACTOR  TIME  CALCULATION 


377  OT$-TIME$ 

378  Q$-“ FACTOR  MATRIX- 

379  PRINT 

380  PRINT  -BEGIN  ";Q$: 

381  X-N 

382  PCTN-X* ( X- 1 ) * ( X+X- 1 ) 


FOR  K-l  TO  N-l 

REM - SEARCH  POR  PIVOT 

T«ZR(K,K)*ZR(K,K)+ZI(K,K)*ZI (K,K) 

Il-K 

FOR  I-K+l  TO  N 

T1 -ZR( I , K ) *ZR( I ,K)+ZI ( I , K) *ZI ( I , K) 

IF  T1 <T  THEN  392 
1 1  - 1 
T-Tl 
NEXT  I 

REM  -  EXCHANGE  ROWS  K  AND  II 

IF  Il-K  THEN  403 
FOR  J»1  TO  N 
T1-ZR(K, J) 

T2-ZI (K, J) 

ZR( K, J ) «ZR(  1 1 ,  J  ) 

ZI ( K, J ) -ZI ( I 1 , J) 

ZR(I1,J)-T1 
ZI ( II, J)-T2 
NEXT  J 
P(K)«I1 

REM  -  SUBTRACT  ROW  K  FROM  ROWS  K+l  TO  N 

FOR  I-K+l  TO  N 

REM  -  COMPUTE  MULTIPLIER  L(I,K) 

T1-(ZR(I,K)*ZR(K,K)+ZI(I,K)*ZI(K,K))/T 
T2- ( ZI ( I,K)*ZR(K,K)-ZR(I,K)*ZI(K,K) )/T 
ZR( I , K) «T1 
ZI ( I , K) -T2 

REM  -  SUBTRACT  ROW  K  FROM  ROW  I 

FOR  J-K+l  TO  N 

ZR(I,J)=ZR(I,J)-(ZR(K,J)*T1-ZI(K,J)*T2) 

ZI ( I, J)-ZI( I, J)-(ZR(K,J)*T2+ZI(K, J)*T1 ) 

NEXT  J 
NEXT  I 
X-N-K 

PCT* 1 -X* ( X- 1 ) * { X+X- 1 ) / PCTN 
GOSUB  1601 
NEXT  K 

REM  -  END  MATRIX  FACTOR  TIME  CALCULATION 


422  T$-TIME$ 

423  GOSUB  1591 

424  PRINT 

425  PRINT  #3,  "FACTOR  MATRIX:  ";T$ 


426  REM  **********  SOLVE  ********** 

427  REM - COMPUTE  RIGHT  HAND  SIDE 


428  FOR  1-1  TO  N 

429  CR(I)»0 

430  CI(I)-0 

431  NEXT  I 

432  FOR  J«1  TO  NS 

433  F2-1/M 

434  IF  C%(E(J),1)  — C%(E(J),2)  THEN  F2-2/M 

435  CR( E ( J ) )«F2*M( J) 

436  CI(E(J)  )— F2*L(J) 

437  NEXT  J 

438  REM  -  PERMUTE  EXCITATION 

439  FOR  K»1  TO  N-l 

440  I 1-P( K) 

441  IF  Il-K  THEN  448 

442  Tl-CR(K) 

443  T2-CI(K) 

444  CR(K)-CR(I1) 

445  CI(K)-CI(I1) 

446  CR( I 1 ) *T1 

447  Cl ( II ) -T2 

448  NEXT  K 

449  REM  -  FORWARD  ELIMINATION 

450  FOR  1*2  TO  N 

451  Tl-0 

452  T2*0 

453  FOR  J*1  TO  1-1 

454  T1*T1+ZR( I, J)*CR(J)-ZI( I, J)*CI(J) 

455  T2»T2+ZR( I, J)*CI(J)+ZI( I, J)*CR(J) 

456  NEXT  J 

457  CR( I )*CR( I )-Tl 

458  Cl ( I ) *CI ( I ) -T2 

459  NEXT  I 

460  REM  -  BACK  SUBSTITUTION 

461  FOR  I-N  TO  1  STEP  -1 

462  Tl-0 

463  T2-0 

464  IF  I-N  THEN  469 

465  FOR  J-I  +  l  TO  N 

466  T1«T1+ZR( I, J) *CR( J)-ZI ( I , J) *CI ( J) 

467  T2«T2+ZR( I,J)*CI(J)+ZI( I,J)*CR(J) 

468  NEXT  J 

469  T«ZR( I , I) *ZR( I , I )+ZI  {  I ,  I  ) *ZI  (  I ,  I ) 

470  T1*CR( I ) -T1 

471  T2-CI ( I ) -T2 

472  CR(I)-(T1*ZR(I,I)+T2*ZI{I,I))/T 

473  Cl ( I )-(T2*ZR( I, I )-Tl*ZI ( I , I ) )/T 

474  NEXT  I 


475  FLG=2 

476  REM  **********  SOURCE  DATA  ********** 

477  PRINT  #3,"  " 

478  PRINT  #3,B$;"  SOURCE  DATA  " ;B$ 

479  PWR=0 

480  FOR  1=1  TO  NS 

481  CR=CR( E ( I ) ) 

482  CI=CI{E( I) ) 

483  T=CR*CR+CI*CI 

484  T1=(L(I)*CR+M(I)*CI)/T 

485  T2=(M(I)*CR-L(I)*CI)/T 

486  02= ( L( I ) *CR+M( I)*CI)/2 

487  PWR=PWR+02 

488  PRINT  *3, "PULSE  " ; E ( I ), "VOLTAGE  =  ( " yL( I ) y " , " yM( I ) y " J) " 

489  PRINT  #3,"  "."CURRENT  =  ( " ;CRy " f " yCI y " J) " 

490  PRINT  #3,"  "."IMPEDANCE  =  ( " y T1 y " , " ; T2 ; " J ) " 

491  PRINT  #3,"  "."POWER  =  ";02;"  WATTS" 

492  NEXT  I 

493  IF  NS> 1  THEN  PRINT  #3,"  " 

494  IF  NS> 1  THEN  PRINT  #3,"TOTAL  POWER  =  " y PWR; "WATTS" 

495  RETURN 

496  REM  **********  PRINT  CURRENTS  ********** 

497  GOSUB  196 

498  S$="N" 

499  PRINT  #3,  "  " 

500  PRINT  #3 , B§ ; "  CURRENT  DATA  "rB$ 

501  FOR  K=1  TO  NW 

502  IF  S$="Y"  THEN  507 

503  PRINT  #3,  "  " 

504  PRINT  #3,  "WIRE  NO. 

505  PRINT  #3,  "PULSE", "REAL", "IMAGINARY", "MAGNITUDE", "PHASE1 

506  PRINT  #3,  "  NO. " , " (AMPS) " , " (AMPS) "."(AMPS) "."(DEGREES)" 

507  N1=N( K, 1 ) 

508  N2=N ( K, 2 ) 

509  I=N1 

510  C=C%(I,1) 

511  IF  ( N1 =0  AND  N2=0 )  THEN  C=K 

512  IF  G=1  THEN  515 

513  IF  ( J 1 ( K) =— 1  AND  N1>N2)  THEN  N2=N1 

514  IF  J1(K)=-1  THEN  525 

515  E%=1 

516  GOSUB  572 

517  121=111 

518  J  2 1 = J 1 1 

519  GOSUB  607 

520  IF  S$="N"  THEN  PRINT  #3, 

1$, Ill ; TAB (29) ;J11 ; TAB (43) ; SI; TAB (57) ;S2 

521  IF  S$="Y"  THEN  PRINT  #1 , I 1 l y " , " ; J1 1 ? " , " ; SI ; " , " ; S2 

522  IF  N1=0  THEN  532 

523  IF  C=K  THEN  525 

524  IF  I$="J"  THEN  N1=N1+1 

525  FOR  I=N1  TO  N2-1 

526  1 2 1 =CR( I ) 
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527  J21=CI(I) 

528  GOSUB  607 

529  IF  S$="N"  THEN  PRINT  #3, 

I , CR( I ) ; TAB ( 29) ; Cl ( I ) ;TAB(43) ;S1?TAB(57) ;S2 

530  IF  S$="Y"  THEN  PRINT  #1 , CR( I ) ; " , " ;CI ( I ) ; " , " ? SI ; “ , " ; S2 

531  NEXT  I 

532  I =N2 

533  C=C%( 1,2) 

534  IF  (N1=0  AND  N2=0)  THEN  C=K 

535  IF  G=1  THEN  537 

536  IF  J1 ( K) =1  THEN  543 

537  E%=2 

538  GOSUB  572 

539  IF  (N1=0  AND  N2=0)  THEN  549 

540  IF  N1 >N2  THEN  549 

541  IF  C=K  THEN  543 

542  IF  I$="J"  THEN  549 

543  12 1 =CR(N2 ) 

544  J2 1 =CI (N2 ) 

545  GOSUB  607 

546  IF  S$=“N"  THEN  PRINT  #3, 

N2 , CR( N2 ) ;TAB( 29) ; Cl (N2 ) ; TAB ( 43 ) ; SI r TAB ( 57 ) ; S2 

547  IF  S$  =  ,'Y,,  THEN  PRINT  *1 ,  CR(  N2 );","?  Cl  ( N2  )?  "  ,  "  r  SI  ;  "  ,  "  ;  S2 

548  IF  J1(K)=1  THEN  554 

549  121=111 

550  J2 1 =J1 1 

551  GOSUB  607 

552  IF  S$="N"  THEN  PRINT 

1 3 , I $ , 1 1 1 ; TAB (29) ;J11 ;TAB(43) ;S1;TAB(57) ?S2 

553  IF  S$="Y"  THEN  PRINT  #1 , II l ; " , " ; J1 1 ; " f " ? SI ; " , " ; S2 

554  IF  S$="Y"  THEN  PRINT  #1,"  1  ,  1  ,  1  ,  1" 

555  NEXT  K 

556  IF  S$="Y"  THEN  569 

557  PRINT 

558  INPUT  "SAVE  CURRENTS  TO  A  FILE  (Y/N)  ";S? 

559  IF  S$="N"  THEN  570 

560  IF  S$  <  > " Y"  THEN  557 

561  PRINT  #3,"  " 

562  INPUT  "FILENAME  (NAME. OUT)  ";F? 

563  IF  LEFT? ( RIGHT$ ( F$ , 4 ) , 1 ) = " . "  THEN  564  ELSE  F$=F$+" .OUT" 

564  IF  0$>"C"  THEN  PRINT  #3, "FILENAME  (NAME. OUT):  " ;F$ 

565  OPEN  F$  FOR  OUTPUT  AS  *1 

566  PRINT  #3,"  " 

567  PRINT  # 1 , NW ; " , " ; PWR ;",C " 

568  GOTO  501 

569  CLOSE  #1 

570  RETURN 

571  REM  -  SORT  JUNCTION  CURRENTS 

572  I$="E" 

573  111=01 

574  Jll=01 

575  IF  ( C=K  OR  C=0)  THEN  580 

576  I$="J" 


577  Ill -CR( I ) 

578  Jll-CI(I) 

579  REM  -  CHECK  FOR  OTHER  OVERLAPPING  WIRES 

580  FOR  J-l  TO  NW 

581  IF  J-K  GOTO  604 

582  LI =N ( J ,  1 ) 

583  L2»N(J,2) 

584  IF  E%-2  THEN  590 

585  CO*C% ( LI , 1 ) 

586  CT»C%(L2,2) 

587  L3-L1 

588  L4=L2 

589  GOTO  594 

590  CO*C% ( L2 , 2 ) 

591  CT-C%(L1,1) 

592  L3-L2 

593  L4-L1 

594  IF  CO— K  THEN  596 

595  GOTO  599 

596  111*111 -CR(L3 ) 

597  J1 1 *J1 1 -Cl (L3 ) 

598  I$-"J" 

599  IF  CT-K  THEN  601 

600  GOTO  604 

601  111*111 +CR( L4 ) 

602  Jll-Jll+CI(L4) 

603  I$-"J" 

604  NEXT  J 

605  RETURN 

606  REM - CALCULATE  SI  AND  S2 

607  131-121*121 

608  J3 l -J2 1 *J2 1 

609  IF  ( 13 l >0  OR  J31 >0)  THEN  612 

610  Sl-01 

611  GOTO  613 

612  Sl*SQR(I3l+J3l ) 

613  IF  I21><0  THEN  616 

614  S2-0 1 

615  RETURN 

616  S2*ATN( J21 /I 2 l }/P0 

617  IF  12 1 >0  THEN  RETURN 

618  S2*S2+SGN( J2 1 ) *180 


619  RETURN 

620  REM  **********  FAR  FIELD  CALCULATION 

621  IF  FLG<2  THEN  GOSUB  196 

622  02-PWR 

623  REM  -  TABULATE  IMPEDANCE 


624  IF  NM-0  THEN  634 

625  FOR  1*1  TO  NM 

626  Z6*T( I ) 

627  Z7--V(I)/(2*P*F*8 . 85E-06) 

628  REM  -  FORM  IMPEDANCE-1 /SQR( DIELECTRIC  CONSTANT) 

629  GOSUB  184 
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630  D=W6 * W6+W7  * W7 

631  Z1(I)=W6/D 

632  Z2(I)=-W7/D 

633  NEXT  I 

634  PRINT  #3,"  “ 

635  PRINT  #3,B$;"  FAR  FIELD  " ;B$ 

636  PRINT  #3,"  " 

637  REM  -  INPUT  VARIABLES  FOR  FAR  FIELD  CALCULATION 

638  INPUT  "CALCULATE  PATTERN  IN  DBI  OR  VOLTS/METER  (D/V)";P$ 

639  IF  P$="D"  THEN  655 

640  IF  P$o"V"  THEN  638 

641  Fl*l 

642  PRINT 

643  PRINT  "PRESENT  POWER  LEVEL  =  ";PWR;"  WATTS" 

644  INPUT  "CHANGE  POWER  LEVEL  (Y/N)  "  ;A$ 

645  IF  A$="N"  THEN  650 

646  IF  A$ < > "Y"  THEN  644 

647  INPUT  "NEW  POWER  LEVEL  (WATTS)  ";02 

648  IF  0$>"C"  THEN  PRINT  #3,"NEW  POWER  LEVEL  =  ";02 

649  GOTO  644 

650  IF  (02<0  OR  02*0)  THEN  02*PWR 

651  F1*SQR(02/PWR) 

652  PRINT 

653  INPUT  "RADIAL  DISTANCE  (METERS)  "  ;  RD 

654  IF  RD<0  THEN  RD-0 

655  A$*" ZENITH  ANGLE  :  INITIAL, INCREMENT, NUMBER" 

656  PRINT  A$; 

657  INPUT  ZA,ZC,NZ 

658  IF  NZ-0  THEN  NZ-1 

659  IF  0$>"CM  THEN  PRINT  *3,A$;n:  " ; ZA; " , " ; ZC; " , " ;NZ 

660  A$  * " AZ I MUTH  ANGLE:  INITIAL, INCREMENT, NUMBER" 

661  PRINT  A$; 

662  INPUT  AA, AC , NA 

663  IF  NA*0  THEN  NA*1 

664  IF  0$>"C"  THEN  PRINT  #3,AS?":  " ; AA; - , " ; AC; " , " ;NA 

665  PRINT  #3,"  M 

666  REM  **********  FILE  FAR  FIELD  DATA  ********** 

667  INPUT  "FILE  PATTERN  (Y/N)";SS 

668  IF  S$*"N"  THEN  676 

669  IF  S$<>"Y"  THEN  667 

670  PRINT  *3,"  " 

671  INPUT  "FILENAME  (NAME . OUT) " ; F$ 

672  IF  LEFTS ( RIGHTS ( F$ , 4 ) , 1 ) ■ "  .  "  THEN  673  ELSE  F$*F$+" .OUT" 

673  IF  0$>"C“  THEN  PRINT  #3, "FILENAME  (NAME. OUT):  ";F$ 

674  OPEN  F$  FOR  OUTPUT  AS  #1 

675  PRINT  #1 , NA*NZ ;  "  ,  " ;02 ; " , “ ; P$ 

676  PRINT  #3,  "  " 

677  K9 l * • 0 16678/PWR 

678  REM  -  PATTERN  HEADER 

679  PRINT  # 3 , B$ ; "  PATTERN  DATA  " ; B$ 

680  IF  P$*"V"  GOTO  685 

681  PRINT  #3, "ZENITH" , "AZIMUTH" , "VERTICAL" , "HORIZONTAL" , "TOTAL" 

682  AS*" PATTERN  (DB)" 
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683  PRINT  #3,”  ANGLE ","  ANGLE “ , A$ , A$ , A$ 

684  GOTO  692 

685  IF  RD>0  THEN  PRINT  #3 , TAB( 15 )  ? "RADIAL  DISTANCE  =  " ; RD ; " 
METERS" 

686  PRINT  #3,TAB( 15) ; "POWER  LEVEL  *  " ; PWR*F1*F1 ; "  WATTS" 

687  PRINT  #3, "ZENITH  AZIMUTH","  E( THETA) 

E( PHI ) “ 

688  A$*M  MAG(V/M)  PHASE (DEG)" 

689  PRINT  #3,"  ANGLE  ANGLE" ,A$,A$ 

690  IF  S$*"Y"  THEN  PRINT  #1,RD 

691  REM  -  LOOP  OVER  AZIMUTH  ANGLE 

692  Q1*AA 

693  FOR  11*1  TO  NA 

694  U3*Q1*P0 

695  V1*-SIN(U3) 

696  V2*COS(U3 ) 

697  REM  -  LOOP  OVER  ZENITH  ANGLE 

698  Q2*ZA 

699  FOR  12*1  TO  NZ 

700  U4*Q2*P0 

701  R3*COS(U4 ) 

702  T3*-SIN(U4) 

703  T1*R3*V2 

704  T2»-R3*V1 

705  Rl— T3*V2 

706  R2*T3*V1 

707  X1»0 

708  Yl-0 

709  Z1*0 

710  X2*0 

711  Y2-0 

712  Z2*0 

713  REM - IMAGE  LOOP 

714  FOR  K*1  TO  G  STEP  -2 

715  FOR  1*1  TO  N 

716  IF  K>0  THEN  718 

717  IF  C% ( I , 1 ) *-C% (1,2)  THEN  812 

718  J*2*W% ( I ) -1+1 

719  REM  -  FOR  EACH  END  OF  PULSE  COMPUTE  A  CONTRIBUTION  TO 

E-FIELD 

720  FOR  F5*l  TO  2 

721  L*ABS(C%( I, F5) ) 

722  F3*SGN(C%(I,F5) )*W*S(L)/2 

723  IF  C% ( 1 , 1 ) < > -C% (1,2)  THEN  725 

724  IF  F3<0  THEN  811 

725  IF  K-l  THEN  728 

726  IF  NM<  >0  THEN  747 

727  REM - STANDARD  CASE 

728  S2-W*(X(J)*R1+Y(J)*R2+Z(J) *K*R3) 

729  Sl*COS(S2) 

730  S2-SIN(S2) 

731  B1-F3* (S1*CR( I)-S2*CI(I)) 

732  B2-F3*(S1*CI ' I )+S2*CR( I ) ) 
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733  IF  C%(I,1)=-C%(I,2)  THEN  742 

734  X1=X1+K*B1*CA(L) 

735  X2=X2+K*B2*CA(L) 

736  Y1=Y1+K*B1*CB(L) 

737  Y2=Y2+K*B2*CB(L) 

738  Z1=Z1+B1*CG(L) 

739  Z2=Z2+B2*CG(L) 

740  GOTO  811 

741  REM - GROUNDED  ENDS 

742  Z1=Z1+2*B1*CG(L) 

743  Z2=Z2+2*B2*CG(L) 

744  GOTO  811 

745  REM  -  REAL  GROUND  CASE 

746  REM  -  BEGIN  BY  FINDING  SPECULAR  DISTANCE 

747  T4=100000 1 

748  IF  R3=0  THEN  750 

749  T4=-Z(J)*T3/R3 

750  B9“T4*V2+X( J) 

751  IF  TB=1  THEN  754 

752  B9=SQR(B9*B9+(Y(J)-T4*V1)“2) 

753  REM  -  SEARCH  FOR  THE  CORRESPONDING  MEDIUM 

754  J2=NM 

755  FOR  J 1 =NM  TO  1  STEP  -1 

756  IF  B9>U(J1)  THEN  758 


757  J2=J1 

758  NEXT  J1 

759  REM  -  OBTAIN  IMPEDANCE  AT  SPECULAR  POINT 

760  Z4=Z1(J2) 

761  Z5=Z2(J2) 

762  REM  -  IF  PRESENT  INCLUDE  GROUND  SCREEN  IMPEDANCE  IN 

PARALLEL 


763  IF  NR=0  THEN  775 

764  IF  B9>U( 1 )  THEN  775 

765  R=B9+NR*RR 

766  Z8=W*R*LOG(R/(NR*RR) )/NR 

767  S8=-Z5*Z8 

768  S9=Z4*Z8 

769  T8=Z4 

770  T9=Z5+Z8 

771  D=T8*T8+T9*T9 

772  Z4* ( S8*T8+S9*T9 ) /D 

773  Z5«(S9*T8-S8*T9)/D 

774  REM  -  FORM  SQR( 1-Z ~2*SIN' 2 ) 

775  Z6*l- ( Z4*Z4-Z5*Z5 ) *T3*T3 

776  Z7=-(2*Z4*Z5)*T3*T3 

777  GOSUB  184 

778  REM  -  VERTICAL  REFLECTION  COEFFICIENT 

779  S8*R3- ( W6*Z4-W7*Z5 ) 

780  S9«-(W6*Z5+W7*Z4) 

781  T8*R3+(W6*Z4-W7*Z5) 

782  T9«W6*Z5+W7*Z4 

783  D»T8*T8+T9*T9 

784  V8= ( S8*T8+S9*T9 ) / D 
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785  V9«(S9*T8-S8*T9)/D 

786  REM  -  HORIZONTAL  REFLECTION  COEFFICIENT 

787  S8-W6-R3*Z4 

788  S9«W7-R3*Z5 

789  T8-W6+R3*Z4 

790  T9«W7+R3*Z5 

791  D-T8*T8+T9*T9 

792  H8«(S8*T8+S9*T9)/D-V8 

793  H9-(S9*T8-S8*T9)/D-V9 

794  REM  -  COMPUTE  CONTRIBUTION  TO  SUM 

795  S2«W*(X(J)*Rl+Y(J)*R2-(Z(J)-2*Ht J2) )*R3) 

796  SI-COS (S2) 

797  S2-SIN(S2) 

798  B1-F3* (S1*CR( I )-S2*CI ( I ) ) 

799  B2-F3* (Sl*CI ( I ) +S2*CR( I ) ) 

800  W6«B1*V8-B2*V9 

801  W7«B1*V9+B2*V8 

802  D-CA( L) *V1+CB( L ) *V2 

803  Z6«D*(B1*H8-B2*H9) 

804  Z7«D*(B1*H9+B2*H8) 

805  X1»X1-(CA(L)*W6+V1*Z6) 

806  X2»X2-(CA(L)*W7+V1*Z7) 

807  Yl-Yl- ( CB( L) *W6+V2*Z6 ) 

808  Y2-Y2- ( CB( L ) *W7+V2*Z7 ) 

809  Z1-Z1+CG(L)*W6 

810  Z2-Z2+CG(L)*W7 

811  NEXT  F5 

812  NEXT  I 

813  NEXT  K 

814  H2-(X1*T1+Y1*T2+Z1*T3)*G0 

815  H1-(X2*T1+Y2*T2+Z2*T3)*G0 

816  X4-(X1*V1+Y1*V2)*G0 

817  X3«(X2*V1+Y2*V2}*G0 

818  IF  P$- "D"  THEN  826 

819  IF  RD-0  THEN  841 

820  H1-H1/RD 

821  H2-H2/RD 

822  X3-X3/RD 

823  X4-X4/RD 

824  GOTO  841 

825  REM - PATTERN  IN  DB 

826  PI— 999 

827  P2-P1 

828  P3-P1 

829  T1-K91* (H1*H1+H2*H2) 

830  T2-K9 1  * ( X3*X3+X4*X4 ) 

831  T3-T1+T2 

832  REM  -  CALCULATE  VALUES  IN  DB 

833  IF  T1 > IE-30  THEN  Pl-4 . 343*LOG(Tl ) 

834  IF  T2> IE-30  THEN  P2-4 . 343*LOG(T2 ) 

835  IF  T3> IE-30  THEN  P3-4 . 343*LOG(T3) 

836  PRINT  #3 , Q2 ; TAB( 15) ;Q1 ; TAB( 29 ) ; PI ; TAB( 43 ) r  P2  r  TAB( 57 ) ; P3 

837  IF  S$-" Y"  THEN  PRINT  # 1 , Q2 ; " f " ; Q1 ; " , " ; PI ; " , " ; P2 ; " , " r  P3 
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838  GOTO  865 


839  REM  -  PATTERN  IN  VOLTS/METER 

840  REM  -  MAGNITUDE  AND  PHASE  OF  E( THETA) 


841  S1=0 

842  IF  ( HI =0  AND  H2=0)  THEN  844 

843  S1=SQR(H1*H1+H2*H2) 

844  IF  HlxO  THEN  847 


845  S2=0 

846  GOTO  850 

847  S2*ATN(H2/H1)/P0 

848  IF  HI <0  THEN  S2=S2+SGN(H2 ) *180 

849* REM  -  MAGNITUDE  AND  PHASE  OF  E(PHI) 

850  S3=0 

851  IF  (X3=0  AND  X4=0)  THEN  853 

852  S3=SQR(X3*X3+X4*X4) 


853  IF  X3>  <0  THEN  856 

854  S4=0 

855  GOTO  858 

856  S4=ATN(X4/X3)/P0 

857  IF  X3<0  THEN  S4*S4+SGN(X4)*180 


858  PRINT  #3, USING 

859  PRINT  #3, USING 

860  PRINT  #3, USING 

861  PRINT  #3, USING 

862  PRINT  #3, USING 


*##.#*  " ;Q2 , Q1 ; 

*#.### - "  ;S1*F1 ; 

*#*.**  " 7 S2 ; 

**.##* - "  ?S3*F1 ; 

###•##" ;S4 

863  IF  S$«“Y"  THEN  PRINT 

#1 , Q2  7 " « “ ;Q1 * " # " *  S1*F1 ? " # " ; S2 ; " , " 7  S3*F1 ? " , "S4 

864  REM  -  INCREMENT  ZENITH  ANGLE 

865  Q2*Q2+ZC 

866  NEXT  12 

867  REM  -  INCREMENT  AZIMUTH  ANGLE 

868  Q1*Q1+AC 

869  NEXT  II 

870  CLOSE  *1 

871  RETURN 

872  REM  **********  NEAR  FIELD  CALCULATION  ********** 

873  REM  -  ENSURE  CURRENTS  HAVE  BEEN  CALCULATED 

874  IF  FLG<2  THEN  GOSUB  196 

875  02 *P WR 

876  PRINT  #3,"  *' 

877  PRINT  #3 , B$  7 "  NEAR  FIELDS  " ; B$ 

878  PRINT  #3,"  " 

879  INPUT  “ELECTRIC  OR  MAGNETIC  NEAR  FIELDS  (E/H)  ";N$ 

880  IF(N$*"H"  OR  N$«"E")  GOTO  882 

881  GOTO  879 

882  PRINT 

883  REM  -  INPUT  VARIABLES  FOR  NEAR  FIELD  CALCULATION 

884  PRINT  "FIELD  LOCATION ( S ): " 

885  A$«" -COORD I NATE  (M):  INITIAL, INCREMENT, NUMBER  " 

886  PRINT  "  X" ; A$ ; 

887  INPUT  XX , XC , NX 

888  IF  NX-0  THEN  NX-1 

889  IF  0$>"CM  THEN  PRINT  #3 , "X" 7 A$ ? " :  "  7XX7  "  ,  "  ?XC?  “  ,  "  rNX 
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890  PRINT  "  Y" ;  A$ ; 

891  INPUT  YY,YC,NY 

892  IF  NY*0  THEN  NY=1 

893  IF  0$>"C"  THEN  PRINT  #3,  "Y"  ; A$  ;  "  s  "  ;  YY ;  "  ,  "  ; YC NY 

894  PRINT  "  Z";A$; 

895  INPUT  ZZ , ZC , NZ 

896  IF  NZ-0  THEN  NZ*1 

897  IF  0$>"C"  THEN  PRINT  #3, "Z" ; A$ ; “ :  " ; ZZ ; " , " ; ZC ; " , " ; NZ 

898  Fl-1 

899  PRINT 

900  PRINT  “PRESENT  POWER  LEVEL  IS  ";PWR;"  WATTS" 

901  INPUT  "CHANGE  POWER  LEVEL  (Y/N)  ";A$ 

902  IF  A$*"N”  THEN  907 

903  IF  A$ <  > " Y"  THEN  901 

904  INPUT  "NEW  POWER  LEVEL  (WATTS)  "  ;02 

905  IF  0$>"C"  THEN  PRINT  #3,"  "sPRINT  #3,"NEW  POWER  LEVEL  (WATTS) 
*  "  ;  02 

906  GOTO  901 

907  IF  (02<0  OR  02*0)  THEN  02-PWR 

908  REM  -  RATIO  OF  POWER  LEVELS 

909  F1*SQR(02/PWR) 

910  IF  N$»"H"  THEN  F1-F1/S0/4/P 

911  PRINT 

912  REM  -  DESIGNATION  OF  OUTPUT  FILE  FOR  NEAR  FIELD  DATA 

913  INPUT  "SAVE  TO  A  FILE  (Y/N)  ";S$ 

914  IF  S$*"N"  THEN  922 

915  IF  S$o"Y"  THEN  913 

916  INPUT  "FILENAME  (NAME. OUT)  ";F$ 

917  IF  LEFT$ ( RIGHTS ( F$ , 4 ) , 1 ) * " • "  THEN  918  ELSE  F$*F$+" .OUT" 

918  IF  0$>"C"  THEN  PRINT  #3,"  "sPRINT  #3, "FILENAME  (NAME. OUT) 


"  ;  F$ 


919 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 

930 

931 

" ;  ZZ 

932 

933 

934 

935 

936 

937 

938 

939 


OPEN  F$  FOR  OUTPUT  AS  #2 
PRINT  * 2 , NX*NY*NZ ; " , " ; 02 ; " , " ; N$ 

REM  -  LOOP  OVER  Z  DIMENSION 

FOR  IZ-1  TO  NZ 

REM  -  LOOP  OVER  Y  DIMENSION 

FOR  IY*1  TO  NY 

REM  -  LOOP  OVER  Z  DIMENSION 

FOR  IX*1  TO  NX 

REM - NEAR  FIELD  HEADER 

PRINT  #3,"  " 

IF  N$*"E"  THEN  PRINT  #3,B$;"NEAR  ELECTRIC  FIELDS"; 
IF  N$-"H"  THEN  PRINT  #3,B$;"NEAR  MAGNETIC  FIELDS"; 
PRINT  #3,TAB(10)  .-"FIELD  POINT:  "  ;  "X  *  ";XX;"  Y  *  " 


PRINT  #3, 
IF  N$*"E" 
IF  N$*"H" 
PRINT  *3, 
Al-0 
A3-0 
A4*0 

REM - 


"  VECTOR" , "REAL" , "IMAGINARY" , "MAGNITUDE" 
THEN  A$*"  V/M  " 

THEN  A$* "  AMPS/M  " 

"  COMPONENT  " , A$ , A$ , A$ , "  DEG " 


LOOP  OVER  THREE  VECTOR  COMPONENTS 


B$ 

B$ 

; YY; "  Z  = 
, "PHASE" 
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940  FOR  1=1  TO  3 

941  X0=XX 

942  Y0=YY 

943  Z0=ZZ 

944  IF  N$="H"  THEN  954 

945  T5=0 

946  T6=0 

947  T7=0 

948  IF  1=1  THEN  T5=2*S0 

949  IF  1=2  THEN  T6=2*S0 

950  IF  1=3  THEN  T7=2*S0 

951  U7*0 

952  U8=0 

953  GOTO  964 

954  FOR  J8*l  TO  6 

955  Kl ( J8, 1 ) *0 

956  Ki ( J8, 2 ) *0 

957  NEXT  J8 

958  J9=l 

959  J8=-l 

960  IF  1=1  THEN  X0»XX+J8*S0/2 

961  IF  1=2  THEN  Y0»YY+J8*S0/2 

962  IF  1=3  THEN  Z0-ZZ+J8*S0/2 

963  REM  -  LOOP  OVER  SOURCE  SEGMENTS 

964  FOR  J*1  TO  N 

965  J 1 =ABS ( C% ( J , 1 ) ) 

966  J2=ABS ( C% ( J , 2) ) 

967  J3-J2 

968  IF  J1>J2  THEN  J3-J1 

969  F4*SGN ( C% ( J , 1 ) ) 

970  F5=SGN(C%( J, 2) ) 

971  F6*l 

972  F7*l 

973  U5=0 

974  U6=0 

975  REM - IMAGE  LOOP 

976  FOR  K»1  TO  G  STEP  -2 

977  IF  C% ( J , 1 ) <  > -C% ( J , 2 )  THEN  983 

978  IF  K<0  THEN  1044 

979  REM  -  COMPUTE  VECTOR  POTENTIAL  A 

980  F6-F4 

981  F7*F5 

982  REM  -  COMPUTE  PSI(0,J,J+.5) 

983  P1«0 

984  P2=2*J3+J-1 

985  P3-P2+.5 

986  P4-J2 

987  GOSUB  75 

988  U1=T1*F5 

989  U2«T2*F5 

990  REM - COMPUTE  PS  I ( 0 , J- . 5 , J ) 

991  P3-P2 

992  P2-P2-.5 
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993  P4=J1 

994  GOSUB  66 

995  V1=F4*T1 

996  V2=F4*T2 

997  REM  -  REAL  PART  OF  VECTOR  POTENTIAL  CONTRIBUTION 

998  X3=U1*CA(J2)+V1*CA(J1) 

999  Y3=U1*CB(J2)+V1*CB(J1) 

1000  Z3=(F7*U1*CG(J2)+F6*V1*CG(J1) )*K 

1001  REM  -  IMAGINARY  PART  OF  VECTOR  POTENTIAL  CONTRIBUTION 

1002  X5=U2*CA(J2)+V2*CA( Jl) 

1003  Y5=U2*CB(J2)+V2*CB(J1) 

1004  Z5={F7*U2*CG(J2)+F6*V2*CG(J1) )*K 

1005  REM  -  MAGNETIC  FIELD  CALCULATION  COMPLETED 

1006  IF  N$="H"  THEN  1038 

1007  D1=(X3*T5+Y3*T6+Z3*T7)*W2 

1008  D2=(X5*T5+Y5*T6+Z5*T7)*W2 

1009  REM  -  COMPUTE  PSI ( . 5 , J , J+l ) 

1010  Pl= . 5 

1011  P2=P3 

1012  P3=P3+1 

1013  P4*J2 

1014  GOSUB  56 

1015  U1=T1 

1016  U2*T2 

1017  REM  -  COMPUTE  PSI ( - . 5 , J , J+l ) 

1018  PI— PI 

1019  GOSUB  56 

1020  U1*(T1-U1)/S( J2) 

1021  U2* (T2-U2 ) /S ( J2 ) 

1022  REM  -  COMPUTE  PSI ( . 5 , J-l , J ) 

1023  PI— PI 

1024  P3-P2 

1025  P2-P2-1 

1026  P4-J1 

1027  GOSUB  56 

1028  U3«T1 

1029  U4»T2 

1030  REM - COMPUTE  PS  I ( - . 5 , J-l , J ) 

1031  PI— PI 

1032  GOSUB  56 

1033  REM  -  GRADIENT  OF  SCALAR  POTENTIAL 

1034  U5* ( U1+ ( U3-T1 )/S(Jl)+Dl) *K+U5 

1035  U6« ( U2+ (U4-T2 ) /S( Jl ) +D2 ) *K+U6 

1036  GOTO  1044 

1037  REM  -  COMPONENTS  OF  VECTOR  POTENTIAL  A 

1038  K I ( 1 , J9 ) *K I ( 1 , J9)+(X3*CR( J)-X5*CI ( J) )*K 

1039  K J ( 2 , J9 ) »K 1 ( 2, J9)  +  (X5*CR( J)+X3*CI ( J) )*K 

1040  K1 (3. J9)-Kl ( 3, J9)+( Y3*CR( J)-Y5*CI ( J) )*K 

1041  Kl ( 4, J9 ) *K l (4,J9)+(Y5*CR(J)+Y3*CI (J) ) *K 

1042  Kl (5,J9)*Kl ( 5 , J9) +( Z3*CR( J )-Z5*CI(J))*K 

1043  Kl (6, J9)-KI (6, J9)+(Z5*CR( J)+Z3*CI (J) )*K 

1044  NEXT  K 

1045  IF  N$«MH"  THEN  1048 
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1046  U7-U5*CR(J)-U6*CI(J)+U7 

1047  U8-U6*CR(J)+U5*CI(J)+U8 

1048  NEXT  J 

1049  IF  N$*"E"  THEN  1071 

1050  REM  -  DIFFERENCES  OF  VECTOR  POTENTIAL  A 

1051  J8-1 

1052  J9-J9+1 

1053  IF  J9«2  THEN  960 

1054  ON  I  GOTO  1055,1060,1065 

1055  H(3)*K1 (5,1)-K1 (5,2) 

1056  H(4)«K1 (6,1)-KI (6,2) 

1057  H(5)*Kl (3,2)-Kl (3.1) 

1058  H(6)-K1 (4,2)-Kl (4,1) 

1059  GOTO  1093 

1060  H(l)*Kl (5,2)-Kl (5,1) 

1061  H(2)*Kl (6,2)-Kl (6,1) 

1062  H ( 5 ) *H ( 5 ) -K 1(1,2) +K 1(1,1) 

1063  H ( 6 ) *H ( 6 ) ~K 1(2,2) +K 1(2,1) 

1064  GOTO  1093 

1065  H(1)*H(1)-KI(3,2)+K1(3,1) 

1066  H ( 2 ) *H ( 2 ) -K 1(4,2) +K 1(4,1) 

1067  H(3)*H(3)+K»(1,2)-K1(1,1) 

1068  H ( 4 ) *H ( 4 ) +K 1(2,2) -K 1(2,1) 

1069  GOTO  1093 


1070  REM  -  IMAGINARY  PART  OF  ELECTRIC  FIELD 

1071  U7*M*U7/S0 

1072  REM  -  REAL  PART  OF  ELECTRIC  FIELD 

1073  U8— M*U8/S0 

1074  REM  -  MAGNITUDE  AND  PHASE  CALCULATION 


1075  S1*0 

1076  IF  (U7*0  AND  U8*0)  THEN  1078 

1077  S1*SQR(U7*U7+U8*U8) 

1078  S2-0 

1079  IF  U8<  >0  THEN  S2*ATN( U7/U8 ) /PO 

1080  IF  U8>0  THEN  1082 

1081  S2-S2+SGN(U7) *180 

1082  IF  I«1  THEN  PRINT  *3,"  X  ”, 

1083  IF  1*2  THEN  PRINT  #3,"  Y  ", 

1084  IF  1*3  THEN  PRINT  #3,"  Z  ", 

1085  PRINT 

#3 , TAB( 1 5 ) ; FI *U8 ; TAB( 29 ) ; FI *U7 ? TAB( 43 ) r FI *S1 ; TAB( 57 ) ;S2 

1086  IF  S$*"  Y"  THEN  PRINT  #2,  F1*U8;  "  ,  "  ;F1*U7  ?  "  ,  "  ;F.l*Sl  ;  "  ,  "  rS2 

1087  REM  -  CALCULATION  FOR  PEAK  ELECTRIC  FIELD 

1088  S1-S1*S1 

1089  S2-S2*P0 

1090  A1*A1+S1*C0S( 2*S2 ) 

1091  A3»A3+S1*SIN(2*S2) 

1092  A4*A4+S1 

1093  NEXT  I 

1094  IF  N$*"E"  THEN  1117 

1095  REM  -  MAGNETIC  FIELD  MAGNITUDE  AND  PHASE  CALCULATION 

1096  FOR  1*1  TO  5  STEP  2 

1097  Sl-0 


1098  IF  (H( I ) *0  AND  H(I+1)=0)  THEN  1100 

1099  S1=SQR( H( I ) *H( I ) +H( 1+1 ) *H( 1+1 )  ) 

1100  S2-0 

1101  IF  H( I ) < >0  THEN  S2=ATN( H( I+1)/H(I) )/P0 

1102  IF  H( I ) >0  THEN  1104 

1103  S2*S2+SGN( H( 1+1 ) ) *180 

1104  IF  1  =  1  THEN  PRINT  *3,"  X  ", 

1105  IF  1=3  THEN  PRINT  #3,"  Y  ", 

1106  IF  1=5  THEN  PRINT  #3,“  Z  ", 

1107  PRINT 

#  3 , TAB (15) ;F1*H( I ) ;TAB( 29) ; F1*H ( 1+1 ) ; TAB( 43 ) ;F1*S1 ;TAB( 57) ;S2 

1108  IF  S$=" Y"  THEN  PRINT 

#2,  F1*H( I ) ; “ , " ;F1*H( I+l ) ; " , " ; F1*S1 ; " , " r S2 

1109  REM  -  CALCULATION  FOR  PEAK  MAGNETIC  FIELD 

1110  S1=S1*S1 

1111  S2=S2*P0 

1112  A1=A1+S1*C0S(2*S2) 

1113  A3=A3+S1*SIN(2*S2) 


1114  A4=A4+S1 

1115  NEXT  I 

1116  REM  -  PEAK  FIELD  CALCULATION 

1117  PK-SQR ( A4/2+SQR { A1 * A1 +A3  * A3 ) / 2 ) 

1118  PRINT  #3,"  MAXIMUM  OR  PEAK  FIELD  =  ";F1*PK;A$ 

1119  IF  ( S$="Y"  AND  NS=“E")  THEN  PRINT  *2 , FI *PK ; " , " ; 02 

1120  IF  (S$*"Y"  AND  N$="H")  THEN  PRINT  #2 , FI *PK ; " , " ; 02 

1121  IF  S$*"Y"  THEN  PRINT  #2, XX; " , " ; YY; “ , " ; ZZ 

1122  REM  -  INCREMENT  X  DIMENSION 

1123  XX-XX+XC 

1124  NEXT  IX 

1125  REM  -  INCREMENT  Y  DIMENSION 

1126  YY-YY+YC 

1127  NEXT  IY 

1128  REM - INCREMENT  Z  DIMENSION 

1129  ZZ-ZZ+ZC 

1130  NEXT  IZ 

1131  CLOSE  #2 

1132  RETURN 

1133  REM  **********  FREQUENCY  INPUT  ********** 

1134  REM - SET  FLAG 

1135  PRINT 


1136  INPUT  "FREQUENCY  (MHZ)";F 

1137  IF  F*0  THEN  F*299.8 

1138  IF  0$>"C”  THEN  PRINT  *3,  "  ":PRINT  #3,  "FREQUENCY  (MHZ):":F 

1139  W-299.8/F 

1140  REM  - VIRTUAL  DIPOLE  LENGTH  FOR  NEAR  FIELD  CALCULATION 

1141  S0«.001*W 

1142  REM  -  1  /  (4  *  PI  *  OMEGA  *  EPSILON) 

1143  M=4 . 7 7783 3 5 2# *W 

1144  REM  -  SET  SMALL  RADIUS  MODIFICATION  CONDITION 

1145  SRM= • 0001 *W 

1146  PRINT  #3,  "  WAVE  LENGTH  =  ”;W;"  METERS" 

1147  REM  -  2  PI  /  WAVELENGTH 

1148  W=2*P/W 


C-23 


1149  W2-W*W/2 

1150  FLG-0 

1151  RETURN 


1152  REM  **********  GEOMETRY  INPUT  ********** 

1153  REM  -  WHEN  GEOMETRY  IS  CHANGED,  ENVIRONMENT  MUST  BE 

CHECKED 


1154  GOSUB  1371 

1155  PRINT 

1156  IF  INFILE  THEN  1162 

1157  INPUT  "NO.  OF  WIRES" ;NW 

1158  IF  NW*0  THEN  RETURN 

1159  IF  NW<  *MW  THEN  1162 

1160  PRINT  "NUMBER  OF  WIRES  EXCEEDS  DIMENSION..." 

1161  GOTO  1157 

1162  IF  0$>**C"  THEN  PRINT  #3,"  “:PRINT  #3,  "NO.  OF  WIRES:" ;NW 

1163  REM  -  INITIALIZE  NUMBER  OF  PULSES  TO  ZERO 

1164  N*0 

1165  FOR  1*1  TO  NW 

1166  IF  INFILE  THEN  GOSUB  1559  .-GOTO  1192 

1167  PRINT 

1168  PRINT  "WIRE  NO.";I 

1169  INPUT  "  NO.  OF  SEGMENTS": SI 

1170  IF  S1*0  THEN  1155 

1171  AS*"  END  ONE  COORDINATES  (X,Y,Z)" 

1172  PRINT  AS; 

1173  INPUT  X1.Y1.Z1 

1174  IF  G<0  AND  Z1<0  THEN  PRINT  "Z  CAN  DT  BE  NEG AT I VE " : GOTO  1172 

1175  AS-"  END  TWO  COORDINATES  (X,Y,Z)M 

1176  PRINT  AS; 

1177  INPUT  X2 , Y2 , Z2 

1178  IF  G<0  AND  Z2<0  THEN  PRINT  "Z  CANNOT  BE  NEGAT I VE ” : GOTO  1176 

1179  IF  X1-X2  AND  Y1*Y2  AND  Z1*Z2  THEN  PRINT"ZERO  LENGTH 
WIRE." :GOTO  1168 

1180  A$*"  RADIUS" 

1181  PRINT  "  "AS; 

1182  INPUT  A( I ) 

1183  IF  A  (  I ) <  *0  1  THEN  1181 

1184  REM  -  DETERMINE  CONNECTIONS 

1185  IF  0$>"C"  THEN  PRINT  #3,"  : P R I NT  #3, "WIRE  NO . " : I 

1186  GOSUB  1301 

1187  PRINT  "CHANGE  WIRE  NO.  "  ;  I  ;  "  (Y/N)  ": 

1188  INPUT  AS 

1189  IF  AS*" Y"  THEN  1167 

1190  IF  AS  <  > " N "  THEN  1187 

1191  REM  -  COMPUTE  DIRECTION  COSINES 

1192  X3-X2-X1 

1193  Y3-Y2-Y1 

1194  Z3-Z2-Z1 

1195  D-SQR(X3*X3+Y3*Y3+Z3*Z3 ) 

1196  CA( I ) -X3/D 

1197  CB( I ) *Y3/D 

1198  CG ( I ) -Z3/D 

1 199  S( I ) * D / S 1 
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1200  REM  -  COMPUTE  CONNECTIVITY  DATA  (PULSES  N1  TO  N) 

1201  N1=N+1 

1202  N( I , 1 ) =N1 

1203  IF  (Sl=l  AND  11=0)  THEN  N(I,1)=0 

1204  N=N1+S1 

1205  IF  11=0  THEN  N=N-1 

1206  IF  12=0  THEN  N=N-1 

1207  IF  N>MP  THEN  PRINT  "PULSE  NUMBER  EXCEEDS 
DIMENSION" : CLOSE: GOTO  1157 

1208  N( I , 2 ) =N 

1209  IF  (Sl=l  AND  12=0)  THEN  N(I,2)=0 

1210  IF  N<N1  THEN  1249 

1211  FOR  J=N1  TO  N 

1212  C%(J,1)=I 

1213  C% ( J , 2 ) =1 

1214  W%(J)=I 

1215  NEXT  J 

1216  C% ( N1 , 1 ) =1 1 

1217  C% (N, 2 ) =12 

1218  REM  -  COMPUTE  COORDINATES  OF  BREAK  POINTS 

1219  Il=Nl+2* ( 1-1 ) 

1220  13=11 

1221  X(I1)=X1 

1222  Y( II ) =Y1 

1223  Z(I1)=Z1 

1224  IF  C% ( N1 , 1 ) =0  THEN  1232 

1225  I 2=ABS ( C% ( N1 , 1 ) ) 

1226  F3=SGN(C%(N1,1) )*S(I2) 

1227  X(I1)=X(I1) -F3*CA( 12 ) 

1228  Y( I 1 ) =Y (II) -F3*CB( 12 ) 

229  IF  C% ( N1 , 1 ) =— I  THEN  F3=-F3 

1230  Z ( I 1 ) =Z ( II ) -F3*CG( 12 ) 

1231  13=13+1 

1232  I6=N+2* I 

1233  FOR  14=11+1  TO  16 

1234  J=I4-I3 

1235  X( I4)=X1+J*X3/S1 

1236  Y( I4)=Y1+J*Y3/S1 

1237  Z( I4)=Z1+J*Z3/S1 

1238  NEXT  14 

1239  IF  C% ( N , 2 ) =0  THEN  1247 

1240  I 2=ABS ( C% ( N, 2 ) ) 

1241  F3=SGN(C%(N, 2) )*S( 12) 

1242  13=16-1 

1243  X( I6)=X( I3)+F3*CA( 12) 

1244  Y ( 16 ) *Y ( 1 3 ) +F3*CB (12) 

1245  IF  I=-C% ( N, 2 )  THEN  F3=-F3 

1246  Z( I6)=Z( I3)+F3*CG( 12) 

1247  GOTO  1257 

1248  REM  -  SINGLE  SEGMEN  0  PULSE  CASE 

1249  I 1 =Nl+2* { I - 1 ) 

1250  X(I1)=X1 

1251  Y( I 1 ) = Y1 
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1252  Z(I1)=Z1 

1253  11=11+1 

1254  X( II ) =X2 

1255  Y( II ) =Y2 

1256  Z(I1)=Z2 

1257  NEXT  I 

1258  REM  **********  GEOMETRY  OUTPUT  ********** 

1259  PRINT  #3,  "  " 

1260  PRINT  #3,  "  ****  ANTENNA  GEOMETRY  ****•• 

1261  IF  N>0  THEN  1266 

1262  PRINT 

1263  PRINT  "NUMBER  OF  PULSES  IS  ZERO. ...  RE-ENTER  GEOMETRY" 

1264  PRINT 

1265  GOTO  1157 

1266  K=1 

1267  J=0 

1268  FOR  1=1  TO  N 

1269  I1=2*W% ( I ) -1+1 

1270  IF  K>NW  THEN  1281 

1271  IF  K=J  THEN  1281 

1272  J=K 

1273  PRINT  #3,"  " 

1274  PRINT  #3, "WIRE  NO.  " ;K;"  COORDINATES" ,,, "CONNECTION  PULSE" 

1275  PRINT  #3 , "X" , "Y" , "Z"  , "RADIUS" , "ENDl  END 2  NO." 

1276  IF  (N(K, 1 ) >  <0  OR  N(K,2)><0)  THEN  1281 

1277  PRINT  #3, -  -  0" 

1278  K=K+1 

1279  IF  K>NW  THEN  1288 

1280  GOTO  1272 

1281  PRINT 

#3.X(I1) ;TAB( 15) ;Y( II ) ;TAB( 29) ; Z ( il ) ;TAB( 43 ) ;A(W%(I) ) ;TAB(57) ; 

1282  PRINT  #3,  USING  "###  ##*  ** " ? C% ( I , 1 ) . C% ( I , 2 ) , I 

1283  IF  ( I=N ( K, 2 )  OR  N( K, 1 ) =N( K, 2 )  OR  C%(I,2)=0)  THEN  K=K+1 

1284  IF  C%(I,1)=0  THEN  C% ( I , 1 ) =W% ( I ) 

1285  IF  C% ( I , 2 ) =0  THEN  C% ( I , 2 ) =W% ( I ) 

1286  IF  ( K=NW  AND  N(K,1)=0  AND  N(K,2)=0)  THEN  1272 

1287  IF  (I=N  AND  K<NW)  THEN  1272 

1288  NEXT  I 

1289  PRINT 

1290  CLOSE  Is  IF  INFILE  THEN  INFILE=0:IF  0$>"C"  THEN  1295 

1291  INPUT  "  CHANGE  GEOMETRY  (Y/N)  ";A$ 

1292  IF  A$="Y"  THEN  1155 

1293  IF  A$ < > "N"  THEN  1291 

1294  REM  -  EXCITATION  INPUT 

1295  GOSUB  1432 

1296  REM  -  LOADS /NETWORKS  INPUT 

1297  GOSUB  1457 

1298  FLG=0 

1299  RETURN 

1300  REM  **********  CONNECTIONS  ********** 

1301  E ( I ) =X1 

1302  L( I ) =Y1 

1303  M( I ) =Z1 


1304  E(I+NW)«X2 

1305  L(I+NW)-Y2 

1306  M{I+NW)«Z2 

1307  G%»0 

1308  11-0 

1309  12-0 

1310  J1(I)»0 

1311  J2(I,1)— I 

1312  J2(I,2)— I 

1313  IP  G-l  THEN  1325 

1314  REM  -  CHECK  FOR  GROUND  CONNECTION 

1315  IF  Zl-0  THEN  1317 

1316  GOTO  1320 

1317  II  — I 

1318  J1(I)— 1 

1319  GOTO  1342 

1320  IF  Z2-0  THEN  1322 

1321  GOTO  1325 

1322  12  —  I 

1323  J1(I)*1 

1324  Gi=l 

1325  IF  1-1  THEN  1360 

1326  FOR  J-l  TO  1-1 

1327  REM  -  CHECK  FOR  END1  TO  END1 

1328  IF  (Xl-E(J)  AND  Yl-L(J)  AND  Zl-M(J))  THEN  1330 

1329  GOTO  1335 

1330  II— J 

1331  J2(I,1)«J 

1332  IF  J2( J, 1 ) — J  THEN  J2(J,1)*J 

1333  GOTO  1342 

1334  REM  -  CHECK  FOR  END1  TO  END2 

1335  IF  (X1*E( J+NW)  AND  Y1»L(J+NW)  AND  Z1»M(J+NW))  THEN  1337 

1336  GOTO  1341 

1337  I1=J 

1338  J2( I , 1 ) =J 

1339  IF  J2(J,2)— J  THEN  J2(J,2)-J 

1340  GOTO  1342 

1341  NEXT  J 

1342  IF  G%=1  THEN  1360 

1343  IF  1=1  THEN  1360 

1344  FOR  J=1  TO  1-1 

1345  REM  -  CHECK  END2  TO  END2 

1346  IF  (X2=E( J+NW)  AND  Y2=L(J+NW)  AND  Z2=M( J+NW) )  THEN  1348 

1347  GOTO  1353 

1348  12— J 

1349  J2( I , 2 ) =J 

1350  IF  J2( J,  2)— J  THEN  J2(J,2)=J 

1351  GOTO  1360 

1352  REM  -  CHECK  FOR  END2  TO  END1 

1353  IF  (X2«E(J)  AND  Y2«L(J)  AND  Z2=M(J))  THEN  1355 

1354  GOTO  1359 

1355  I2-J 

1356  J2( I , 2 ) =J 
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1357  IP  J2( J, 1 )*-J  THEM  J2(J,1)-J 

1358  GOTO  1360 

1359  MEXT  J 

1360  PRIMT  #3 ,  "  COORDINATES" , "  "."  "."END 

NO.  OP" 

1361  PRINT  #3,"  X","  Y","  Z". "RADIUS  CONNECTION 

SEGMENTS" 

1362  PRINT  #3 . XI ; TAB( 15) ;  Y1 ; TAB( 29) ; Z1 ; TAB( 57 ) ;  1 1 


1363  PRINT 

#3 . X2 ;TAB( 15) ;Y2;TAB(29) ;Z2;TAB(43) ; A( I ) ;TAB( 57 ) ; 12 ; TAB( 71 ) ; SI 

1364  RETURN 

1365  REM  **********  ENVIROMENT  INPUT  ********** 

1366  PRINT 

1367  PRINT  "  ****  WARNING  ***♦" 

1368  PRINT  "REDO  GEOMETRY  TO  ENSURE  PROPER  GROUND 
CONNECT ION/ D I S CONNECT I ON " 

1369  PRINT 

1370  REM  -  INITIALIZE  NUMBER  OF  RADIAL  WIRES  TO  ZERO 

1371  NR*0 

1372  REM  -  SET  ENVIRONMENT 


1373  PRINT  #3,"  " 

1374  A$* "ENVIRONMENT  (+1  POR  FREE  SPACE,  -1  FOR  GROUND  PLANE)" 

1375  PRINT  A$ ; 

1376  INPUT  G 

1377  IF  0$>"C"  THEN  PRINT  *3,A$;":  ";G 

1378  IF  G*1  THEN  1430 

1379  IF  Go-1  THEN  1375 

1380  REM - NUMBER  OF  MEDIA 

1381  A$»"  NUMBER  OF  MEDIA  (0  FOR  PERFECTLY  CONDUCTING  GROUND)" 

1382  PRINT  A$; 

1383  INPUT  NM 

1384  IF  NM<*MM  THEN  1387 

1385  PRINT  "NUMBER  OF  MEDIA  EXCEEDS  DIMENSION..." 

1386  GOTO  1382 

1387  IF  0$>"C"  THEN  PRINT  #3,A$;":  " ; NM 

1388  REM  -  INITIALIZE  BOUNDARY  TYPE 

1389  TB=1 

1390  IF  NM*0  THEN  1430 

1391  IF  NM*1  THEN  1398 

1392  REM  -  TYPE  OF  BOUNDARY 

1393  A$*"  TYPE  OF  BOUNDARY  ( 1-LINEAR,  2-CIRCULAR)" 

1394  PRINT  "  " ;  A$ ; 

1395  INPUT  TB 

1396  IF  0$>"C"  THEN  PRINT  *3,A$;" :  " ;TB 

1397  REM  -  BOUNDARY  CONDITIONS 

1398  FOR  1*1  TO  NM 

1399  PRINT  "MEDIA"; I 

1400  A$="  RELATIVE  DIELECTRIC  CONSTANT,  CONDUCTIVITY" 

1401  PRINT  "  ";A$; 

1402  INPUT  T(I),V(I) 

1403  IF  0$>"C"  THEN  PRINT  #3, A?;":  " ;T( I ) " , "V( I ) 

1404  IF  I>1  THEN  1416 

1405  IF  TB=1  THEN  1416 
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1406  A$«"  HUMBER  OP  RADIAL  WIRES  IH  GROUND  SCREEN" 

1407  PRINT  "  " ;A$; 

1408  INPUT  NR 

1409  IP  0$>“C"  THEN  PRINT  #3,A$r":  " ;NR 

1410  IP  NR-0  THEN  1416 

1411  A$-"  RADIUS  OP  RADIAL  WIRES" 

1412  PRINT  "  "  ?A$  r 

1413  INPUT  RR 

1414  IP  0$>"C"  THEN  PRINT  f3,A$;":  "  ?  RR 

1415  REM  -  INITIALIZE  COORDINATE  OF  MEDIA  INTERPACE 

1416  U(I)«1000000l 

1417  REM  -  INITIALIZE  HEIGHT  OP  MEDIA 

1418  H(I)«0 

1419  IP  I-NM  THEN  1424 

1420  A$«"  X  OR  R  COORDINATE  OF  NEXT  MEDIA  INTERFACE" 

1421  PRINT  "  ";A$; 

1422  INPUT  U(I) 

1423  IP  0$>"C“  THEN  PRINT  #3,A$r“ :  M?U(I) 

1424  IF  1-1  THEN  1429 

1425  A$*"  HEIGHT  OP  MEDIA" 

1426  PRINT  "  " ; A$ ; 

1427  INPUT  H( I ) 

1428  IP  0$>"C"  THEN  PRINT  #3,A$r":  ";H(I) 

1429  NEXT  I 

1430  RETURN 

1431  REM  **********  EXCITATION  INPUT  ********** 

1432  PRINT 

1433  A$-"NO.  OF  SOURCES  " 

1434  PRINT  A$ ; 

1435  INPUT  NS 

1436  IF  NS< 1  THEN  NS-1 

1437  IF  NS < *MP  THEN  1440 

1438  PRINT  "NO.  OF  SOURCES  EXCEEDS  DIMENSION  ..." 

1439  GOTO  1434 

1440  IF  0$>"C"  THEN  PRINT  *3,"  ":PRINT  #3,  A$ ; " :  ";NS 

1441  FOR  1-1  TO  NS 

1442  PRINT 

1443  PRINT  "SOURCE  NO.  " ; I ; " s " 

1444  A$«" PULSE  NO. ,  VOLTAGE  MAGNITUDE,  PHASE  (DEGREES)" 

1445  PRINT  A$; 

1446  INPUT  E ( I ) , VM , VP 

1447  IF  E(I)<»N  THEN  1450 

1448  PRINT  "PULSE  NUMBER  EXCEEDS  NUMBER  OF  PULSES..." 

1449  GOTO  1445 

1450  IF  0$>"C”  THEN  PRINT  #3,A$;":  " ; E( I ) " , "VM" , "VP 

1451  L(I)«VM*COS(VP*PO) 

1452  M(I)-VM*SIN(VP*PO) 

1453  NEXT  I 

1454  IF  FLG>2  THEN  FLG-1 

1455  RETURN 

1456  REM  **********  LOADS  INPUT  ********** 

1457  PRINT 

1458  INPUT  "NUMBER  OF  LOADS  ”;NL 


C-29 


1459  IF  NL<*ML  THEN  1462 

1460  PRINT  "NUMBER  OF  LOADS  EXCEEDS  DIMENSION..." 

1461  GOTO  1458 

1462  IF  0$>"C"  THEN  PRINT  *3, "NUMBER  OF  LOADS" ?NL 

1463  IF  NL< 1  THEN  1494 

1464  INPUT  "S-PARAMETER  (S-jw)  IMPEDANCE  LOAD  (Y/N)";L$ 

1465  IF  L$<  > "Y"  AND  L$<>"N"  THEN  1464 

1466  A$* "PULSE  NO. , RESISTANCE. REACTANCE" 

1467  IF  L$«"Y"  THEN  A$«  "PULSE  NO.,  ORDER  OF  S-PARAMETER 
FUNCTION" 

1468  FOR  I»1  TO  NL 

1469  PRINT 

1470  PRINT  "LOAD  NO.  "  ;  I ;  "  «  " 

1471  IF  L$«"Y"  THEN  1478 

1472  PRINT  A$; 

1473  INPUT  LP {  I )  , LA (1,1,1), LA (2,1,1) 

1474  IF  LP { I ) >N  THEN  PRINT  "PULSE  NUMBER  EXCEEDS  NUMBER  OF 
PULSES...":  GOTO  1472 

1475  IF  0$>"C"  THEN  PRINT  #3,A$;" : 

" ; LP (I);","? LA (1,1,1);","; LA (2,1,1) 

1476  GOTO  1493 

1477  REM  -  S-PARAMETER  LOADS 

1478  PRINT  A$; 

1479  INPUT  LP ( I ) , LS { I ) 

1480  IF  LP( I ) >N  THEN  PRINT  "PULSE  NUMBER  EXCEEDS  NUMBER  OF 
PULSES...":  GOTO  1478 

1481  IF  LS ( I ) >MA  THEN  PRINT  "MAXIMUM  DIMENSION  IS  10":GOTO  1479 

1482  IF  0$>"C"  THEN  PRINT  *3,A$r":  " ? LP ( I )  ;  "  ,  "  ; LS ( I ) 

1483  FOR  J*0  TO  LS(I) 

1484  A$=" NUMERATOR,  DENOMINATOR  COEFFICIENTS  OF  S~" 

1485  PRINT  A$;J; 

1486  INPUT  LA( 1 , I , J ) , LA( 2 , I , J ) 

1487  IF  0$>"C"  THEN  PRINT  #3 , A$ r J ; " : " ; LA( 1 , I , J ) ? " , " ? LA( 2 , I , J ) 

1488  NEXT  J 

1489  IF  LS ( I ) >0  THEN  1493 

1490  LS ( I ) = 1 

1491  LA( 1 , I , 1 ) *0 

1492  LA ( 2 , I , 1 ) *0 


1493  NEXT  I 

1494  FLG=0 

1495  RETURN 

1496  REM  **********  main  PROGRAM  ********** 

1497  REM  -  DATA  INITIALIZATION 

1498  REM - PI 

1499  P=4*ATN(1) 

1500  REM  -  CHANGES  DEGREES  TO  RADIANS 

1501  P0*P/180 

1502  B$="********************" 

1503  REM  -  INTRINSIC  IMPEDANCE  OF  FREE  SPACE  DIVIDED  BY  2  PI 

1504  G0*29. 979221# 

1505  REM  -  Q-VECTOR  FOR  GAUSSIAN  QUADRATURE 

1506  READ 

Q(l) ,Q{2) ,Q(3) ,Q(4) ,Q(5) ,Q(6) ,Q(7) ,0(8) ,Q(9) , Q( 10 ) , Q( 1 1 ) , Q( 12 ) 
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1507  READ  Q( 13) ,Q( 14) 

1508  DATA 

.288675135, .5, .430568156, .173927423, .169990522, .326072577 

1509  DATA  .480144928, .050614268, .398333239, .111190517 

1510  DATA  .262766205, .156853323, .091717321, .181341892 

1511  REM  -  E-VECTOR  POR  COEFFICIENTS  OF  ELLIPTIC 

INTEGRAL 

1512  READ  C0,C1,C2,C3,C4,C5,C6,C7,C8,C9 

1513  DATA 

1.38629436112, .09666344259, .03590092383, .03742563713, .01451196212 

1514  DATA  .5, . 12498593397, .06880248576, .0332835346, .00441787012 

1515  REM  -  IDENTIFY  OUTPUT  DEVICE 

1516  GOSUB  1582 

1517  PRINT  #3 , TAB( 20) ; B$ ; B$ 

1518  PRINT  #3, TAB( 22 ); "MINI-NUMERICAL  ELECTROMAGNETICS  CODE" 

1519  PRINT  #3,TAB( 36) ? “MININEC" 

1520  PRINT  #3 , TAB( 24) ; DATE$ ;TAB( 48 ) ; TIMES 

1521  PRINT  #3 , TAB( 20) ; B$ ; B$ 

1522  REM  -  FREQUENCY  INPUT 

1523  GOSUB  1135 

1524  REM  -  ENVIRONMENT  INPUT 

1525  GOSUB  1371 

1526  REM  -  CHECK  FOR  NEC-TYPE  GEOMETRY  INPUT 

1527  GOSUB  1552 

1528  REM  -  GEOMETRY  INPUT 

1529  GOSUB  1155 

1530  REM - MENU 

1531  PRINT 

1532  PRINT  B$;"  MININEC  MENU  ";B$ 

1533  PRINT  "  G  -  CHANGE  GEOMETRY  C  -  COMPUTE/DISPLAY 

CURRENTS" 

1534  PRINT  "  E  -  CHANGE  ENVIRONMENT  P  -  COMPUTE  FAR-FI ELD 
PATTERNS" 

1535  PRINT  "  X  -  CHANGE  EXCITATION  N  -  COMPUTE  NEAR-FIELDS" 

1536  PRINT  "  L  -  CHANGE  LOADS" 

1537  PRINT  "  F  -  CHANGE  FREQUENCY  Q  -  QUIT" 

1538  PRINT  B$ ; B$ ; B$ 

1539  INPUT  "  COMMAND  ";C$ 

1540  IF  C$-"F"  THEN  GOSUB  1135 

1541  IF  C$-"P"  THEN  GOSUB  621 

1542  IF  C$*"X"  THEN  GOSUB  1432 

1543  IF  C$«"E"  THEN  GOSUB  1366 

1544  IF  C$-"G"  THEN  GOSUB  1154 

1545  IF  C$-"C"  THEN  GOSUB  497 

1546  IF  C$-"L"  THEN  GOSUB  1457 

1547  IF  C$*"N"  THEN  GOSUB  874 

1548  IF  C$ < > "Q"  THEN  1531 

1549  IF  0$*"P"  THEN  PRINT  #3,  CHR$(12)  ELSE  IF  0$*"C"  THEN  PRINT 
♦  3,  "  " 

1550  CLOSE 

1551  GOTO  1619 

1552  REM  **********  NEC-TYPE  GEOMETRY  INPUT  ********** 

1553  OPEN  "MININEC. INP"  AS  #1  LEN»30 
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1554  FIELD  *1,2  AS  S$,4  AS  Xl$,4  AS  Yl$,4  AS  Zl$,4  AS  X2$,4  AS 
Y2$ ,  4  AS  Z2$ , 4  AS  RS 

1555  GET  1 

1556  NW-CVI(S$) 

1557  IF  MW  THEM  INFILE-1 

1558  RETURN 

1559  REM - GET  GEOMETRY  DATA  FROM  MIN  I  NEC .  INP 

1560  GET  1 

1561  S1-CVI(S$) 

1562  X1-CVS<X1$) 

1563  Y1-CVS(Y1$) 

1564  Z1-CVS(Z1$) 

1565  X2«CVS(X2$) 

1566  Y2»CVS(Y2$) 

1567  Z2*CVS(Z2$) 

1568  A( I ) -CVS ( R$ ) 

1569  IF  G<0  THEN  IF  Z1<0  OR  Z2<0  THEN  GOSUB  1574 

1570  PRINT  #3 , M  "sPRINT  *3, "WIRE  NO."?I 

1571  IF  X1-X2  AND  Y1-Y2  AND  Z1-Z2  THEN  PRINT"WIRE  LENGTH  IS 
ZERO.": GOTO  1549 

1572  GOSUB  1301 

1573  RETURN 

1574  IF  IZNEG  THEN  1578 

1575  PRINT "NEGATIVE  Z  VALUE  ENCOUNTERED  FOR  GROUND  PLANE." 

1576  INPUT  "ABORT  OR  CONVERT  NEGATIVE  Z  VALUE  TO  ZERO  (A/C)?  ";A$ 

1577  IF  A$»"A"  THEN  1549  ELSE  IF  A$-"C"  THEN  IZNEG-1  ELSE  1576 

1578  IF  Z1<0  THEN  Zl— Z1 

1579  IF  Z2<0  THEN  Z2—Z2 

1580  RETURN 

1581  REM  **********  IDENTIFY  OUTPUT  DEVICE  ********** 

1582  INPUT  "OUTPUT  TO  CONSOLE,  PRINTER,  OR  DISK  (C/P/D)" rO$ 

1583  IF  0$»"C“  THEN  F$-"SCRN: " sGOTO  1588 

1584  IF  0$-"P"  THEN  F$«"LPT1 sGOTO  1588 

1585  IF  0$  < > "D"  THEN  1582 

1586  INPUT  "FILENAME  ( NAME. OUT )" ;F$ 

1587  IF  LEFTS ( RIGHTS ( F$ , 4 ) , 1 ) ■ " . "  THEN  1588  ELSE  FS«F$+".OUT" 

1588  OPEN  F$  FOR  OUTPUT  AS  *3 

1589  CLS 

1590  RETURN 

1591  REM  **********  CALCULATE  ELAPSED  TIME  ********** 

1592  IH-VAL(MID$(TS,1,2) )-VAL(MIDS(OT$, 1,2) ) 

1593  IM-VAL ( MI D$ ( T$ ,4,2) ) -VAL ( MI D$ ( OT$ ,4,2) ) 

1594  IS-VAL(MID$(T$,7,2) )-VAL(MIDS(OTS,7,2) ) 

1595  IF  IS<0  THEN  IS-IS+60 : IM-IM-1 

1596  IF  IM<0  THEN  IM-IM+60 : IH-IH-1 

1597  IF  IH<0  THEN  IH-IH+24 

1598  T$-" : "+MID$(STR$( IS+100) ,3) 

1599  IF  IH  THEN  T$«MID$ ( STRS ( IH) , 2 ) +" s "-t-MIDS ( STR$ ( IM+100) , 3 ) +T$ 
ELSE  TS-MIDS(STR$(IM),2)-»-T$ 

1600  RETURN 

1601  REM  **********  CALCULATE  APPROXIMATE  TIME  REMAINING 

1602  I PCT-100*PCT 
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1603  T$«TIME$ 

1604  IH-VAL(MID$(T$,1,2) )-VAL(MID$ (OT$ , 1 . 2 ) ) 

1605  IF  IH<0  THEM  IH-IH+24 

1606  IM*VAL ( MI D$ ( T$ ,4,2)) -VAL ( MI D$ ( OT$ ,4,2)) 

1607  IS-VAL(MID$(T$,7,2) )-VAL(MID$(OT$, 7, 2 ) ) 

1608  IS«IS+60*(IM+60*IH) 

1609  IS* IS* ( 1/PCT-l ) 

1610  IM*INT( IS/60) 

1611  IS*IS  HOD  60 

1612  IH-INT(IM/60) 

1613  IM-IM  MOO  60 

1614  T$*'' :  "+MID$  ( STR$  ( IS+100 )  ,  3 ) 

1615  IF  IH  THEN  T$»MID$ ( STR$ ( IH) , 2 ) +M : "+MID$ ( STR$ ( IM+100 ) , 3 ) +T$ 
ELSE  T$ “MI D$ ( STR$ ( IM ) ,2) +T$ 

1616  LOCATE  CSRLIN.l 

1617  PRINT  Q$ ;  I PCT ;  " I  COMPLETE  -  APPROX  TIME  REMAINING  ,,T$" 

1618  RETURN 

1619  END 
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